January 1961 - Vol. 15, No. 73 
| 
(atiGAN 


ye" sony 


Mathematiea™™ 
of 
Computation 





A journal devoted to advances in numerical analysis, 
he application of computational methods, mathematical tables, 
high-speed calculators and other aids to computation 





Formerly: Mathematical Tables and other Aids to Computation 


Published Quarterly by the 
National Academy of Sciences—National Research Council 








Editorial Committee 
Division of Mathematics 
National Academy of Sciences—National Research Council 
Washington, D. C. 


H. rr, ene, 7 Mathematics Laboratory, David Taylor Model 


Washington 7 , D.C. 

ALAN cen Reina University of Liverpool, Liverpool 3, England 

P. C Hammer, University of Wisconsin, Madison 6, Wisconsin 

Evuaene Isaacson, New York University, New York 3, New York 

Y. L. Luks, Midwest Research Institute, Kansas City ‘10, Missouri 

Dantet SHanxs, Applied Mathematics Laboratory, David Taylor Model Basin, 
Washington 7, D. C. 

A. H. Tavs, University of Illinois, Urbana, Illinois 

R. 8. Varea, Case Institute of Technology, Cleveland 6, Ohio 

J. W. Wrencu, Jr., Applied Mathematics Laboratory, David Taylor Model Basin, 
Washington 7, D. C. 

D. M. Youne, Computing Center, University of Texas, Austin 12, Texas 


Information to Subscribers 


The journal is published quarterly in one volume per year with issues numbered 
serially since Volume I, Number 1. Starting with January, 1959 subscriptions are 
$8.00 per year, single copies $2.50. Other issues are available as follows: 
Volume I (1943-1945), Nos. 10 and 12 only are available; $1.00 per issue. 
Volume Il (1946-1947), Nos. 13, 14, 17, 18, 19, and 20 only available; $1.00 


per issue. 

Volume III (1948-1949), Nos. 21-28 available. $4.00 per year (four issues), 
$1.25 per issue. 

Volume IV—XII (1950 through 1958), all issues available; $5.00 per year, $1.50 
per issue. 


Microcard Edition 


Volumes I—X (1943-1956), Nos. 1-56 are now available on Microcards and may be 
purchased from the Microcard Foundation, Box 2145, Madison 5, Wisconsin, at a 
cost of $20.00 for the complete set. Succeeding volumes are available on request. 


Information to Contributors 


All contributions intended for publication in Mathematics of Computation and all 
books for review should be addressed to H. Polachek, Technical Director, Applied 
Mathematics Laboratory, David Taylor Model Basin, Washington 7, D. C. The 
author may suggest an appropriate editor for his paper. Manuscripts should be 
typewritten double-spaced in the format used by the journal. For journal abbrevia- 
tions, see Mathematical Reviews, v. 19, Index, 1958, p. 1417-1430, with the excep- 
tion of MTAC and Math. Comp. Authors should submit the original and one copy, 
and should retain one copy. 


Subscriptions, address changes, business communications and payments should be 
sent to: 


THE PRINTING AND PUBLISHING OFFICE 
Tue NatTionaL ACADEMY OF SCIENCES 
2101 Constitution Avenue 

Washington 25, D. C 


PUBLISHED BY THE 
NATIONAL ACADEMY OF SCIENCES—NATIONAL RESEARCH COUNCIL 
Mt. Royal and Guilford Avenues, Baltimore 2, Md. 


Second-class postage paid at Baltimore, Md. 











Applications of the Complex Exponential 
Integral 


By Murlan S. Corrington 


1. Introduction. The recent publication of an extensive table of the exponential 
integral for complex arguments [1] makes it possible to evaluate a large number of 
indefinite integrals not in existing tables, and to obtain values for the sine and cosine 
integrals for complex arguments. 


2. Definition of Exponential Integral. The definition used by the National 
Bureau of Standards will be used throughout, 


(1) E,(z) = / «du = RE,(z) + iIE,(z) 


where z = x + iy. R and J denote the real and imaginary components respectively. 
The integral converges if the upper limit is «e**, and is independent of a, so long as 
—3n S a S 32 [2]. To make £,(z) a single-valued function a branch cut is made 
just below the negative real axis, including the origin, such that z = x + iy = 


pe”, —x < 0 < x. This means that when the integral is evaluated for a point on the 
| negative real axis the contour must be indented above the pole at the origin. 
The values of £,(z) are given in the tables for the region 0 < @ S =; those for 
the region —x < @ < 0 can be obtained from the relations 


E,(2) = E,(z), z= 2 — ty. 
Here, as usual, f(z) means the complex conjugate of f(z). 


3. Relation to Exponential Integral for Real Arguments. The earlier tables of 
the exponential integral [3] for real arguments use the following definitions: 


(2) Ei(x) = ¢ < du = -4 ay 
—«o U —z U 


and 


(3) —Ei(—2) - | — du, 


where f is the Cauchy principal value. This gives the relations 

(4) E\(y) = —Ei(-y), y>0 
and 

F (5) —E,(—y) = Ei(y) + ix, y>0o 


4. Sine and Cosine Integrals for Complex Arguments. Let u = it, then equation 
(1) becomes 
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00 —it 


E,(z + iy) = / ie 


y—iz t 


/ cos t dt — i| sin ty 
y—iz t y—iz t 


— Ci(y — ix) — ile — Sil(y — iz)). 


(6) 


Interchanging letters gives the relations 


(7) E,\(y — tx) = —Ci(ax + ty) + alga — Si(x + ty)] 
and 
(8) E\(-—y + ix) = —Ci(x + ty) — if[gae — Si(a + ty)). 


Subtracting equations (7) and (8) gives 
Si(x + iy) = 44 [E\(y + ix) — Ei\(—y + izx)| + 4x 
(9) = 3(7Ei(y + tx) + TE\(—y + iz) + 2] 
+ #i[RE,(y + ix) — RE\(—y + iz)]. 
Adding equations (7) and (8) gives 
Ci(x + ty) = —3[RE\(y + ix) + RE, (—y + iz)] 
vies + HUE (y + tx) — TE\(—y + ix)). 


In the other quadrants the following relations can be used: 


(11a) Si(az — ty) = Si(x 4 iy), 
(11b) Si(—x + iy) = —Si(x + iy), 
(11c) Si(—x — ty) = Si(r + ty), 
(12a) Ci(x — iy) = Ci(x + iy), 
(12b) Ci(—x + iy) = Ci(a + ty) + im, 
(12c) Ci(—x — ty) = Ci(x + ty) — ix. 
When x = 0 in equation (6) 

E,\(ty) = —Ci(y) — ila — Si(y)], 
so 
(13) Si(y) = 3a + IE,(ty) 
and 
(14) Ci(y) = —RE,(ty). 


5. Relation to Previously Published Tables. Tables of sine and cosine integrals 
for complex arguments in rectangular coordinates have been published by Bleick 
[4]. His definition of the sine integral is not the one given here, and his function is 
not an analytic function. Since he defined the sine integral as 
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z+iy - 
(15) Sis(x + iy) = [ et at, 


ty 
the relation to the one used here is 


Si(x + iy) = Si(iy) + Sis(x + iy) 
Sis(x + ty) + HI[EA(y) — Ei(—y)] + 3. 


His definition of the cosine integral leads to the same values as the one here. 

A table of #,(z) for arguments in polar coordinates has been published by 
Mashiko for angles in the first quadrant [5]. Computer routines for the complex 
exponential integral are also available (6). 


(16) 


6. Table of Integrals. The following integrals, in which a and b are real, can be 
evaluated in terms of the real and imaginary components of the exponential integral: 
b 


sin t e ' , . 
| rap dt = 5 IRE(a) — RE,(a + 1b)! 





(17) 





in(—a) — RE,(—a + 1b)}, 





b —a a 
(18) [ ona dt = [x + IEx(—a + )] — F 1E(a +4), 
0 


b ‘ 
(19) [ seee-3 ie + IE(- —a + ib)] + 5 cTEx(a + ib), 


bar 


(20) | a = 5 e“[RE\(a) — RE,( ib)] 
0 a+ 
+ 2¢°(RE,(—a) — RE, ib)], 











i [ = fat = 29 [REA b + ia) + REs(—b + ia) + 2Ci(a)] 

21 

; pe > [E,(b + ia) + IE,(—b + ia)— 2Si(a) + 2], 
f ane a= mS —— [REx(b + ia) —RE; (—b + ia)] 

22 

ane 4+ 8° ITE,(b + ia) — IE,(—b + ia)], 

hens [ ae dt = = 5 608 a [RE,(b + ia) — RE,(—b + ia)] 


me 2 sin a [TE,(b + ia) -_ IE,(-— b+ ia)|, 
tcosht, 1. ’ * i 
(24) l e+e dt = 5 Sin a [[E,(b + ia) + IE\(—b + ia) — 2Si(a) + x] 
— 4 cosa[RE,(b + ia) +RE,(—b + ia) + 2Ci(a)], 





b t 


[ appu=- = — 924 pe (—b + ia) + Ci(a)] 
0 





on 
(25) cos a 


_ 008.4 | IE,(—b + ia) — Si(a) + nae 
a 2 
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b 


e' sin a 
[ ———, dt = —— [RE,(b + ia) + Ci(a)] 
2 2 
(9) 72? 


COS @ 





+ | ra(o + ia) — Si(a) + rl, 


a “ ¥ ; 1 
(27) [ Pie dt = sina | 1-0 + ia) — Si(a) + 5 r| 
— cos a [RE,(—b + ia) + Ci(a)], 


| 


i icgg? ' 
(28) l e+e dt = sina Ex + ia) — Si(a) + 5 | 
— cos a(RE,(b + ia) + Ci(a)], 





(29) [ oe a a 5 (RE(b) — RE,(—b)], 


b ~ 
(go) a = log yb + 5 [RE\(b) + RE;(—6)], 
0 


——__———————— @ = hg7 + = ~~ + b’) 


[ 1 — cos at cosh bt 
(31) o t 
4 (RE,(b + ia) + RE\(—b + ia)| 


where log y = 0.5772156649 --- = Euler’s constant, 


(32) [ sin af sinh at = tan” — 5 (rEi(b + ta) — 1B,(—b + ia)], 


(33) [ “— dt = = 5 (REC + ia) — RE,(—b + ia)), 
sin at cosh bt : ; : 
(34) [ ee dt = 5 EAC + ts) +3RA=b + ic). +o, 
0 
b 
(35) | sin ¢ log (a’ + #) dt = log a’ — cos b log (a’ + b’) 
‘ 0 
+ e* [RE,(a) — RE,(a + ib)] + e * [RE,(—a) —RE,(—a + 1b)], 
b 
(36) [ cos t log (a’ + ¢’) dt = sin b log (a’ + b’) 
0 


—e”([E,(—a + ib) + x] — e’ [E\(a + tb), 
b 


[ sint t log (a’ + ¢’) dt = cosh b log (a’ + b’) — loga’ 
0 
(37) — sin a [1E,(b + ia) + IE,(—b + ia) — 28i(a) + x] 
+ cosa [RE,\(b + ta) + RE,(—b + ia) + 2Ci(a)], 


b 


[ cosh ¢ log (a’ + ¢’) dt = sinh b log (a’ + b’) 
: 0 
(38) + sina [JE,(b + ia) — IE,(—b + ia)] 


— cos a|[RE,(b + ia) — RE\(—b + ia)), 





de 


(4 
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P t b &€ 
[ sin t tan’ — dt = —cosb tan’ — — — IE,(a + ib) 
(39) 0 a a 2 r 
| + Sx + [Esa + iby), 
4 t b, é 
| cos ttan™* — dt = sinb tan — + © [RE,(a + ib) — RE,(a)] 
)I, — 5 [RE\(—a + ib) — RE\(—a)], 
b : 
I sinh tan! dt = cosh b tan’? — “9 (RE,(b + ia) — RE: —b + ia)] 
a)}, (41) 7 a a 2 
~ > [IE,(b + ia) — 1E,(—b + ia), 
b 
[ cosh ¢ tan” “ dt = sinh b tan” : + sin a Ci(a) + cos al 3 r— sia) | 
0 — 
(42) + 5" [RE\(b + ia) + REs(—b + ia)] 
+ 5+ UEs(b + ia) + 1E(—b + éa)], 
ia)| a b l 
(43) [ : = tt = —tan~ = + IE,\(—a + ib) + =, a> 0, 
0 
1 —set > 
(44) [ o_ dt = tan” : + IE\(a + ib), a>0, 
0 
1 -—et 2 
(45) [ a ——S dt = | log ( + 5) ~ RE,(a) + RE,(a + ib), 
0 
P 9 9 9 
| e' log (a + t°) dt = e’ log (a + b*) — loga’ 
() 
(46) — 2sin a [JE,(—b + ia) — Si(a) + 31] 
»)| + 2 cos a[RE,(—b + ia) + Ci(a)], 
ib), “ 
[ e‘ tan” s dt = e’ tan “8 + sin a[RE,(—b + ia) + Ci(a)] 
(47) Jo a a 
+ ib), + cos allE,(—b + ia) — Si(a) + 4x}. 
The integrals of equations (19) and (20) have been tabulated by Bleick [4]. 
7. Extension of the Table of Integrals. Other more complex integrals can be 
derived from the above table by differentiation or integration with respect to 
Yi(a)|, aor b, using the Cauchy-Riemann relations, 
ORE,\(a + ib) — AE, (a + tb) 
0a “ 0b 
(48) at 
e 
= —- |b sin b — a cos b] 
+ ia)), a + b 
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and 
ORE,(a + ib) BS _d1E,(a + ib) 
ab be — | 
(49) it 
é . 
= eypl-4 sin b — b cos b]. 
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Evaluation of the Incomplete Gamma Function 
of Imaginary Argument by 
Chebyshev Polynomials 


By Richard Barakat 


During the course of some work on the diffraction theory of aberrations it was 
necessary to evaluate numerically the incomplete gamma function of imaginary 
argument y(», ix) for certain values of the parameter v. These integrals are special 
cases of the confluent hypergeometric function, and in the standard notation [1] 


a) y(v, ix) = [ e t’ dt Re(v) > 0 


= (ix)’v Fil», 1 + »v, —izx) 


where ;F; is the confluent hypergeometric function. Only the case of v real is con- 
sidered. 

When » isan integer, ,/, issimply a polynomial in «x. It is also possible to evaluate 
(1) in terms of the Fresnel integrals when (v — 1) is a half-integer. For vy = 4, (1) 
is proportional to the Fresnel integral. For other half-integer values, we can use the 
recurrence relation 


(2) y(1 + v, ix) = vy(», ix) — (ix)’e * 


to generate the necessary formulas. An alternate procedure utilizes the Lommel 
functions of two variables [2]; unfortunately, the Lommel functions have not been 
extensively tabulated. 

The remaining values of v can be treated by a Taylor series expansion for 
small values of x and by asymptotic developments for large values of x. The main 
difficulty is the intermediate range where x is neither large nor small. A possible 
method is to use Nielsen’s [3] representation 


a 
(3) yy, ix + iy) — y(v, ix) = & (ix) (1 — v),(—ix) “[1l — e “e,(ty)] 


n=0 


where e, (zy) is the truncated exponential series 


n . _ 
(4) en(iy) = Cy) : 
m=0 Mm. 
This method seems unnecessarily complicated. 

A more powerful method is to utilize the Chebyshev polynomials. Instead of 
expanding the exponential in (1) into a Taylor series, we expand into a series of 
Chebyshev polynomials. As Lanczos [4] has pointed out, ““While the expansion into 
powers on the basis of the Taylor’s series gives the slowest convergence, the expan- 
sion into the Chebyshev polynomials gives the fastest convergence.” 

It is convenient to transform the integral slightly by setting t = irq in (1) 


Received March 4, 1960. 
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1 
(5) y(v, iz) = (ia)’ [ e ‘%q"* dq. 
0 
Expanding the exponential into the series [5] 


(6) et = Di én(—2)"Jn(2)Ta(q) lq| 1 


where ¢, is Newmann’s factor (¢ = 2; ¢«, = 1 forn = 1) and T,,(q) isthe Chebyshev 
polynemial of the first kind 





To(q) =] 
(7) _ 71 (/2] (_1)*(n — 8 — 1)! we 
T.(9) = 5 2 ain ay. 


Substituting (6) into (5) and interchanging sum and integral (this is permissible 
because of uniform convergence of the series) leads to the integrals 


Bor) = . 
v 
(8) B,(») =| T.(q)q’ dq 


= n [n/2] (—1)*(n ey 1)! gn 
- 2 & sl(n — 2s)!(n — 28 + pv) © 





Finally (5) becomes 
(9) x(¥, ix) = (ig) én(—t)"Jn(2)B,(r). 


To obtain an error estimate on B, for large n, we use the relation 
(10) T,n(q) = cosné, q = cosé 
so that (8) becomes 





r/2 
(1) B,(v) = l cos n6(cos 0)” ' sin 6 dé 
~ 0(1). 
The asymptotic expression for J,(z) when n is large and z is fixed is given by [6] 
n 2 
(12) Ina) ~ GRY (j « n). 
n!\ 4 


Hence the dominating factor is J,(x). 


As an example of the power of the method, we consider the case of the Fresnel 


integrals (vy = 4), using as a standard the recent table by Pearcey [7]. The first 


thirteen values of B,(.5) are: 


By = 2.0 

B, = 0.66666667 B, = —0.49696970 
B, = —1.20 Bs = 0.42976370 
B; = —0.85714286 By = 0.40710869 
By, =  0.57777778 By = —0.40752711 
Bs = 0.52813853 By = —0.38709818 
B, = —0.54358974 By = 0.35566473 
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Omitting the superfluous factor (iq)* in (5) yields 


v(.5, iz) = y(5, ix)(ig)~* 
(13) = 2 «(—8)"Ja(2)Ba(5) 


a (7) (C(x) — iS(z)). 


The results for various values of x are tabulated below. The numbers in parentheses 
indicate the number of terms used in the series expansion. 


Relyi(.5, iz)] 
z=1 z=2 z=3 z=4 
Pearcey 1.8090490 1.335193 0.811910 0.461462 


series 1.80904840 (4) 1.3351931 (5) 0.8119099 (6) 0.4614622 (7) 
1.80902564 (3) 1.3351744 (4) 0.8118993 (5) 0.4614576 (6). 


Similar results hold for the imaginary part of y, . The small number of terms 
required for six-place accuracy is remarkable compared with the usual Taylor 
series expansion. 


Using this method, integrals have been computed corresponding to v = .2, .4, 
.6, 8 in the interval 0(.1)5. The integrals in question are 


1 
(14) vi(v, iz) = [ eq’ dq 


= (iq) ’y(», ix) 


(see equation 5). 
The first fifteen values of B,(v) are listed in Table 1, while in Table 2 we give 
the values of the real and imaginary values of 7;(», ix) to six decimals. 























TABLE 1 

B,/» 2 | 4 6 | 8 

~one sae! aes Binns Both. H. 
Bo 500000000 2.50000000 | 1.66666667 | 1. 25000000 
B, . 83333333 71428571 62500000 55555556 
B, —4.09090909 —1.66666667 — .89743590 | —.53571429 
B; —1.25000000 — .96638655 | — .76388889 | — 61403509 
B, 3.26839827 98484848 | 32887403 | 05952381 
B; .99358973 65203859 | 42658730 27324056 
Be ~3.08546292 — 90909091 | —.32988751 | —.11554623 
B; — .98290598 — .62150387 | — .39839182 — .25726340 
Bs 286921706 76839826 | 23742520 | 05315763 
By .90301003 .53065850 | 31171679 | 18010639 
By —2.76714376 — .72594072 | —.23077198 | —.06744618 
By | —.88701757 — 50877478 | —.29504438 | —.17180411 
By | 2.65095170 65829331 | 19072285 | 04254734 
By | 8439206 .4620702 | 2532754 | . 1368987 
Bu | —2.581791 — .630248 —.185156 | —.048825 
Bis | — .624846 — .335736 | —.187154 | —.108511 
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TABLE 2 
Refyi(r, éx)] 

x/» “Types Bee 6 8 x/» 2 

0 — |5.000000'2.500000|1.666667| 1.250000! 0 0 
.1  |4.997728'2.497918|1.664744| 1.248215) .1 | —.083281) 
.2 |4.990925|2.491682|1.658989| 1.242871) .2 | —.166251 
.3 |4.979626|2.481327|1.649432| 1.233999] .3 | —.248598| 
4 |4.963889|2.466908|1.636128| 1.221650 .4 | —.330016 
5 |4.943798|2.448505|1.619153| 1.205896 .5 | —.410206 
.6 |4.919457|2.426217|1.598600| 1.186830} .6 | —.488874) 
7 |4.890092|2. 400165] 1.574586) 1.164560, .7 | —.565736| 
-8 |4.858551|2.370489|1 .547245| 1.139217) .8 | —.640519) 
.9 |4.822300/2.337349|1.516730| 1.110945) .9 | —.712964| 
1.6 |4.782427|2.300922|1.483209| 1.079908] 1.0 | —.782825| 
1.1 |4.739134|2.261403|1.446869| 1.046282) 1.1 | —.849872 
1.2 |4.692643|2.219001|1.407910| 1.010259, 1.2 | —.913890) 
1.3  |4.643186|2. 1739391366545, .972044| 1.3 | —.974687| 
1.4 |4.591014/2. 126454/1.322998] .931850} 1.4 |—1.032085| 
1.5 |4.536384|2.076793|1.277506| .889904 1.5 |—1.085929) 
1.6 |4.479566|2.025210|1.230312| .846438| 1.6 |—1.136085| 
1.7 |4.420838|1.971971)1.181667| .801691| 1.7  |—1.182439| 
1.8 |4.360484/1.917344|1.131828| .755908| 1.8 |—1.224902| 
1.9 |4.298791|1.861603|1 .081054| .709337| 1.9 |—1.263406| 
2.0 |4.236050|1 .805022|1.029608| .662227/ 2.0 |—1.297904' 
2.1 |4.172552|1.747879| .977750| 614827, 2.1 |—1.328374| 
2.2 |4.108588|1.690447| .925742| .567385| 2.2 |—1.354815| 
2.3 |4.044445|1.632998| .873840| .520143| 2.3 |—1.377251 
2.4 |3.980406|1.575798| .822296| .473342| 2.4 |—1.395725 
2.5 |3.916747|1.519107| .771356| .427214| 2.5 |—1.410303 
2.6 |3.853735|1.463177| .721257| .381983) 2.6 |—1.421070 
2.7 |3.791628|1.408251| .672227| .337863| 2.7 |—1.428133| 
2.8 |3.730674|1.354558| .624483| .295059| 2.8 |—1.431615| 
2.9 |3.671107|1.302318| .578228| .253763| 2.9 |—1.431659) 
3.0 |3.613147|1.251735| .533655| .214153} 3.0 |—1.428423) 
3.1 |3.556998)1 203000} .490939| .176394| 3.1 |—1.422082 
3.2 |3.502850/1.156285| .450241| .140635| 3.2 |—1.412821 
3.3 |3.450873/1.111749| .411707| .107010| 3.3 |—1.400841 
3.4 |3.401222|1.069530| .375462| .075636| 3.4 |—1.386353! 
3.5 |3.3540301.029750| .341618| .046611| 3.5 —1.369574) 
3.6 |3.309413| .992512} .310265| .020018, 3.6 |—1.350732| 
3.7 |3.267467| .957899| .281477|—.004080| 3.7 |—1.330060 
3.8 |3.228268) 925076) .255308|—.025637| 3.8 |—1.307793 
3.9 |3.191871| .896787| .231794|—.044625| 3.9 |—1.284170 
4.0 |3.158312| .870358| .210951|—.061033, 4.0 |—1.259431| 
4.1 |3.127608| .846696| .192780|—.074871| 4.1 |—1.233815| 
4.2 |3.099754| .825790| .177259|—.086163| 4.2 |—1.207556) 
4.3 |3.074730| .807608| .164354|— .094951)_ 4.3 |—1.180887) 
4.4 |3.052494| .792104| .154010|—.101294 4.4 |—1.154033| 
4.5 |3.032988| .779212| .146158/—.105264| 4.5 |—1.127213) 
4.6 |3.016138| .768853| .140713|—.106951| 4.6 |—1.100636) 
4.7 |3.001852| .760931| .137576|—.106457| 4.7 |—1.074503| 
4.8 |2.990025| .755335| .136634|—.103895| 4.8 |—1.049002 
4.9 |2.980539| .751945| .137763|—.099392| 4.9 |—1.024311 
5.0 |2.973261| .750625 -140826|— .093084| 5.0 |—1 








.000593 


Im[yi(v, ix)] 
A 6 
0 | 0 | 0 
— .071380|— .062454)|—. 
— .142465|— .124630|—. 
— .212966|— .186254|—. 
— .282593] — .247052|—. 
— .351063|— .306759|—. 
— .418102|— .365115|—. 
— .483443| — .421868) —. 
— .546831|— .476778|—. 
— .608020| — 529616) = 
— .666783|— .580166} —. 
— .722903| — .628226| — 
— .776182| — .673611|—. 
— .826440|— .716151|—. 
— .873513| — .755697| — 
— .917259|— .792115|—. 
— .957555|— .825293}—. 
— .994298|— .855137|—. 
—1.027407|— .881576|—. 
—1.056824| — |904557|—. 
—1.082511|— .924050|—. 
—1.104452|— .940044|—. 
— 1. 122653|— .952550)—. 
—1.137140|— .961600| — . 
—1.147963| — .967244|—. 
—1.155190|— .969553}—. 
—1.158907|— .968616) —. 
—1.159223|— .964540)—. 
—1.156262|— 957449| —. 
—1.150165|— .947483| — 
—1.141088| — .934797|—. 
— 1.129203) — .919558| — . 
—1.114693|— .901946| — . 
—1.097754| — .882152| —. 
—1.078591|— .860375|—. 
— 1.057416] — .836824|—. 
—1.034451)— .811711|—. 
—1.009920| — .785256}—. 
— .984051|— .757678) — . 
— .957077| — .729200)—. 
— .929226| — .700045| — 
— .900728| — .670431|—. 
— .871809|— 640578) — . 
— .842691|— .610695|—. 
— .813590| — .580990| — 
— .784712|— .551660|—. 
— .756258| — .522895) — . 
— .728417|— .494873|—. 
— .701367|— .467764|—. 
— .675274|— .441722|—. 





— .650291)|— . 416892) — 





055512 
110761 
165486 
219430 
272340 
323971 
374084 
422454 


- 468863 


513108 


.554999 


594363 
631041 


-664892 


695794 
723642 
748351 
769857 
788113 
803095 
814797 
823233 
828439 
830468 
829390 
825297 
818295 
808507 


. 796071 


781139 
763875 
744456 
723067 
699904 
675167 
649064 
621806 
593607 
564681 


.535242 


505501 
475667 
445944 


- 416526 


387604 
359357 
331955 
305558 
280310 


. 256347 
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New Tables of Howland’s and Related Integrals 


By C. W. Nelson 


1. Introduction. In an earlier paper by C. B. Ling and the present author [1], 
values of the four integrals, 





(1) [ w* dw a (k 2 1) 
is "an smhwtw ki ics k a + 22’ (k 2 3) 
(2) II; fs “dw _ 2 xe dx (k = 1) 
II,* = EH EEGs 3 wtw =—<k! sinh 2x + 2x’ (k 2 3) 


were tabulated to 6D. These four integrals are called Howland’s integrals because 
they first appeared in Howland’s papers dealing with the stresses in a perforated 
strip (2, 3]. In another paper, C. B. Ling has reproduced the 6D tables of Howland’s 
integrals from [1] and added tables of other integrals derived from them [4]. 


2. Values of Howland’s Integrals to 9D. For most purposes, 6D tables of 
Howland’s integrals are adequate. This is usually true in plane-stress problems of 
elasticity involving straight bars of rectangular cross-section, for example. How- 
ever, the author has found that even in such problems the stresses are more easily 
evaluated in all regions of interest in the bar without resorting to contour integra- 
tion if 9D rather than 6D tables of Howland’s integrals are available. The 9D 
values of Howland’s integrals given in Table 1 were computed by the author and 
have been checked by using the following checking formulas [1]: 





> (1 — Ineus) = I, — 5, bm 2k(1 — In) = 212 — 3, 
k=0 k=l 
> (Iteus — 1) = 4, Di 2k( Tae —1)=¥, 
> Ile = 4 — Ih, > 2kIIy = 4 — 2M, 
> Ian = 3, > 2kII = 3. 
\ k=l k=2 


The method of computing the integrals is explained in [i] and this explanation will 
not be repeated here. 


3. Values of Jo, Is, and |? | to 18D. Recently the author has encountered a 


need for still more precise values of some of the Howland’s integrals. This occurs in 
the Hankel transform or Fourier-Bessel integral solutions [5, 6, 7) for problems of 
elasticity involving axially symmetric loading of a thick plate of infinite radius, 
Received June 6, 1960. 
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NEW TABLES OF HOWLAND’S AND RELATED INTEGRALS 13 
TABLE 1 
5 Values of Howland’s Integrals to 9D 
k | h he | Ih | Ihe 
1 0.768 574 538 | 0.220 119 581 | 
1}, 2 | 0.767 847 439 | 0.087 927 235 
3 | 0.827 710 296 | 2.038 710 667 | 0.043 347 862 | 0.460 713 719 
4 | 0.883 506 807 | 1.353 294 115 | 0.022 583 004 | 0.099 315 532 
1) 5 0.925 475 998 | 1.156 864 366 | 0.011 923 473 | 0.032 412 690 
3) 6 0.954 191 562 | 1.076 729 764 | 0.006 287 972 | 0.012 616 908 
7 0.972 698 993 | 1.039 251 312 | 0.003 295 013 | 0.005 391 117 
1) 8 0.984 124 180 | 1.020 537 600 | 0.001 713 298 | 0.002 433 000 
3) 9 0.990 949 179 | 1.010 870 147 | 0.000 884 147 | 0.001 135 997 
10 0.994 922 440 | 1.005 784 842 | 0.000 453 215 | 0.000 542 202 
nant 11 0.997 188 575 | 1.003 084 774 | 0.000 231 010 | 0.000 262 590 
ted 12 0.998 459 958 | 1.001 644 976 | 0.000 117 209 | 0.000 128 431 
1d’s 13 0.999 163 823 | 1.000 876 180 | 0.000 059 254 | 0.000 063 239 
14 0.999 549 306 | 1.000 465 841 | 0.000 029 870 | 0.000 031 283 
15 0.999 758 559 | 1.000 247 139 | 0.000 015 025 | 0.000 015 525 
” 16 0.999 871 321 | 1.000 130 809 | 0.000 007 546 | 0.000 007 722 
dss 17 | 0.999 931 719 | 1.000 069 073 | 0.000 003 785 | 0.000 003 847 
s of 18 | 0.999 963 903 | 1.000 036 390 | 0.000 001 897 | 0.000 001 919 
ow- 19 | 0.999 980 979 | 1.000 019 128 | 0.000 000 950 | 0.000 000 958 
sily 20 | 0.999 990 006 | 1.000 010 034 | 0.000 000 476 | 0.000 000 478 
nnd 21 | 0.999 994 762 | 1.000 005 252 | 0.000 000 238 | 0.000 000 239 
9D 22 0.999 997 261 | 1.000 002 744 | 0.000 000 119 | 0.000 000 119 
23 0.999 998 570 | 1.000 001 431 | 0.000 000 060 | 0.000 000 060 
and 24 | 0.999 999 255 | 1.000 000 745 | 0.000 000 030 | 0.000 000 030 
25 0.999 999 613 | 1.000 000 388 | 0.000 000 015 | 0.000 000 015 
26 0.999 999 799 | 1.000 000 201 | 0.000 000 007 | 0.000 000 007 
27 0.999 999 896 | 1.000 000 104 | 0.000 000 004 | 0.000 000 004 
28 | 0.999 999 946 | 1.000 000 054 | 0.000 000 002 | 0.000 000 002 
29 | 0.999 999 972 | 1.000 000 028 | 0.000 000 001 | 0.000 000 001 
30 | 0.999 999 986 | 1.000 000 014 | 0.000 000 000 | 0.000 000 000 
31 | 0.999 999 993 | 1.000 000 007 
32 | 0.999 999 996 | 1.000 000 004 
33 | 0.999 999 998 | 1.000 000 002 
34 | 0.999 999 999 | 1.000 000 001 
35 | 0.999 999 999 | 1.000 000 001 
36 | 1.000 000 000 | 1.000 000 000 
will where integrals such as 
- («int 7 + = cosh *) Jo(px)Ji(axr) dx 
d (4) : a [ »  s 2 
— ie a(sinh x + 2) 
Ss ° 
ve . must be evaluated. In equation (4), Jo(px) and J;(ax) are Bessel functions of 
' dius order zero and unity, respectively. Sneddon [5, 6] describes an approximate method 
244 of evaluating such integrals. He and his associates appear to have treated only cases 
where the double sign in the denominator has the plus value, but, presumably, their 
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method could be extended to cases where the double sign takes the minus value. 
Sneddon’s approximate method does not involve the use of Howland’s integrals. 

For a limiting case of the integral in equation (4), namely a = 0 and double sign 
taken as plus, Sneddon found that his approximate method gave errors of the order 
of 1% for p = 0 to 3. However, the author found that Sneddon’s approximate 
method applied to the integral in equation (4) gave an error of about 7% for 
p = 0, a = 1, and double sign taken plus. Comparable errors in certain other inte- 
grals were obtained by the approximate method so that, in a certain thermal stress 
problem (radius of heated region equal to thickness of plate), the errors in normal 
stresses on the axis of the plate ranged up to 21.8% of the greatest normal stress 
occurring anywhere on the axis or 9.3% of the greatest normal stress occurring 
anywhere in the plate. 

The author has felt that it is desirable to be able to obtain the values of integrals 
such as those considered in the foregoing two paragraphs with an accuracy better 
than that obtainable by Sneddon’s approximate method. The more accurate values 
may at least be used to check a few values obtained by the approximate method in 
cases of doubt, and the availability of accurate values of the integrals occurring in 
thick-plate problems may even be found essential in the extension of the Hankel 
transform method to problems not yet considered. The author first attempted to 
evaluate the stresses in a thick-plate problem [7] with the aid of the 9D values of 
Howland’s given in Table 1 and found that only a very limited range of values of 
the parameters p and a in integrals such as that of equation (4) could be success- 
fully dealt with. For this reason, the 18D values of two of the four Howland’s 
integrals, J, and J',*, for even integral k, were computed. These are given in Table 2 
and discussed further in the following. Green and Willmore [9] encountered integrals 
similar to those discussed in the present paper except that theirs contained only 
one Bessel function. They used Howland’s integrals for the simpler integrals they 
had to evaluate and Sneddon’s approximate method for the rest. Apparently the 
Fourier-Bessel integral method was first applied to thick plate problems by Lamb 
[10], who made no attempt to evaluate the integrals. Dougall [11] proposed that 
integrals such as the example in equation (4), but containing only one Bessel 
function, be evaluated by contour integration, but apparently no one has cared to 
perform the task involved which appears to be considerable even though the roots 
of sinh z + z = 0 are known. 

If the integral taken as an example in equation (4) is expressed as a series of 
Howland integrals by expanding the numerator of the integrand in powers of z, 
only even powers of x will occur in this expansion. Thus the integral can be evalu- 
ated if adequate values of J, or [,* are available for even integral k, and if the series 
converges or can somehow be summed. From a consideration of various loading 
conditions of a thick plate, it appears to the author that as long as the axially- 
symmetric loading of the plate consists only of boundary loads on the plane faces 
(i.e. no body forces), all other integrals needed can also be evaluated from tables 
of J; and J;,* for even integral k. Accordingly, the values of J;, and J;,* given to 18D 
in Table 2 are believed to be the only Howland integrals needed in a fairly broad 
class of thick-plate problems. Again the reader is referred to [1] for the general 
method of computing the integrals in Table 2. However, it should be stated here 
that the method depends on the basic equation 
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ie. TABLE 2 
Values of the Howland Integrals In and I, to 18D 
zn k Ik i,* 
ler - ae ed 
ite 2 0.767 847 439 133 919 047 
‘or + 0.883 506 806 508 692 590 1.353 294 115 170 484 009 
“ 6 0.954 191 561 826 139 064 1.076 729 763 674 217 127 
r 8 0.984 124 180 148 424 614 1.020 537 600 064 959 340 
= 10 0.994 922 439 853 445 374 1.005 784 842 224 523 505 
al 12 0.998 459 957 947 832 383 1.001 644 976 253 520 390 
28S 14 0.999 549 305 562 626 177 1.000 465 841 012 174 418 
ng 16 0.999 871 321 426 371 778 1.000 130 808 809 455 335 
18 0.999 963 903 069 165 323 1.000 036 389 774 380 024 
| 20 0.999 990 005 851 207 354 1.000 010 033 628 030 387 
aed 22 0.999 997 260 778 826 414 1.000 002 744 455 935 026 
ter 24 0.999 999 255 282 148 058 1.000 000 745 402 213 094 
1es 26 0.999 999 798 878 365 743 1.000 000 201 210 025 403 
in 28 0.999 999 945 988 927 609 1.000 000 054 022 369 865 
in 30 0.999 999 985 565 214 614 1.000 000 014 436 216 216 
kel 32 0.999 999 996 158 384 196 1.000 000 003 841 795 571 
. 34 0.999 999 998 981 377 142 1.000 000 001 018 645 284 
to 36 0.999 999 999 730 790 958 1.000 000 000 269 211 822 
of 38 0.999 999 999 929 059 585 1.000 000 000 070 940 758 
of 40 0.999 999 999 981 355 380 1.000 000 000 018 644 662 
= 42 | 0.999 999 999 995 111 469 1.000 000 000 004 888 537 
d’s 44 0.999 999 999 998 721 023 1.000 000 000 001 278 977 
a 46 0.999 999 999 999 666 045 1.000 000 000 000 333 955 
e2 48 0.999 999 999 999 912 959 1.000 000 000 000 087 041 
‘als 50 0.999 999 999 999 977 351 1.000 000 000 000 022 649 
niv 52 0.999 999 999 999 994 116 1.000 000 000 000 005 884 
hear 54 0.999 999 999 999 998 473 1.000 000 000 000 001 527 
h 56 0.999 999 999 999 999 604 1.000 000 000 000 000 396 
the 58 0.999 999 999 999 999 898 1.000 000 000 000 000 102 
mb 60 | 9.999 999 999 999 999 974 1.000 000 000 000 000 026 
hat 62 | 0.999 999 999 999 999 993 1.000 000 000 000 000 007 
ssel 64 | 0.999 999 999 999 999 998 1.000 000 000 000 000 002 
| to 66 1.000 000 000 000 000 000 1.000 000 000 000 000 000 
ts 
- n+l 
s of (5) ->, “= * - 
An where 
alu- © ntk-t 
ries (6) : a 2 [ zx dx 
ling n 2(k!) Jo sinh" x 
lly- and that the values in Table 2 have been checked by using the second and fourth 
Aces of equations (3). When these checks were first applied, errors were indicated. To 
bles eliminate these errors, they had to be located and this was facilitated by applying 
i8D another checking formulat 
oad + For » = 1, use 
eral " * . 
here (7a) 2 ([>"| - :) = 7 — In2 
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= 2m | 2m 1 
(7) ee el Oe cB ho al 


Equation (7) is obtained by writing the series of integrals represented by the left 
member of the equation, interchanging the order of summation and integration, 
performing the summation which requires merely recognition of the Maclaurin 
expansion of cosh x — =—s * and then performing the evaluation of the resulting 
elementary integral as follows: 


1 f° 2” (x cosh x — sinh x) dr a" \* 1 
() 5/ ~ 








sinh” x e ~ 2(n — 1) sinh™™ z |o 2(n — 1)° 


TABLE 3 
Veluce of|* | to 18D 





























El - | el 
1 1.051 799 790 264 645 000 21 | 0.036 396 239 976 639 766 
2 0.450 771 338 684 847 857 22 | 0.034 712 213 240 700 944 
3 0.283 656 164 817 032 302 23 =| : 0.033 177 094 523 417 174 
4 0.206 411 435 020 280 449 24 0.031 771 973 852 978 737 
5 0.162 091 502 769 567 468 25 0.030 481 012 454 106 163 
6 0.133 389 997 047 394 398 26 0.029 290 843 691 027 776 
7 0.113 302 756 992 243 931 27 0.028 190 108 942 707 765 
8 0.098 463 515 074 533 508 28 0.027 169 094 240 899 480 
9 0.087 055 830 293 630 464 29 0.026 219 443 043 962 235 
10 0.078 014 012 373 164 903 30 0.025 333 927 165 280 539 
11 0.070 671 880 431 666 952 31 0.024 506 262 569 448 302 
12 0.064 591 720 775 032 384 32 0.023 730 960 108 696 094 
13 0.059 474 138 718 298 180 33 0.023 003 203 705 225 661 
14 0.055 107 451 486 079 519 34 0.022 318 750 267 406 497 
15 0.051 337 774 971 088 869 35 0.021 673 846 947 047 716 
16 0.048 050 565 748 631 057 36 0.021 065 162 331 082 628 
17 0.045 158 809 284 589 149 37 0.020 489 728 904 903 199 
18 0.042 595 221 342 283 185 38 0.019 944 894 690 604 315 
19 0.040 306 957 521 379 715 39 0.019 428 282 397 575 541 
20 0.038 251 937 970 666 201 | 40 0.018 937 754 758 470 797 
TABLE 4 
7 : © xJi(ox) dx 
Values of | ich ots 
© xJi(px) dx 

,; | rwee. 

0 0 

1 0.471 282 

2 0.250 403 

3 0.166 664 

4 0.125 000 

5 0.100 000 








ol 
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left 
ion, 
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166 
44 
i74 
137 
163 
176 
765 
480 
235 
D39 
302 
094 
661 
497 
716 
628 
199 
315 
541 
797 
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The values of [=| are recorded in Table 3, since they are the starting point for 
obtaining J, and [,*, and because they may be needed as in the example leading to 
Table 4 in the following. It was not felt necessary to include values of [| for 


m = 4, 6, 8, --- because, at worst, only a few such values might be needed by a 
reader, and these can be computed from Table 3 by using equation (16) given in 
the following. 


4. Application of Howland’s Integrals in Evaluating Related Integrals. Pro- 
cedures for evaluating integrals such as those in equation (4) for moderate values 
of p and a, say for p and a each less than unity, have been amply discussed in con- 
nection with similar integrals treated in references [2], [3], [4] and [7]. If either p or 
a (or both) should be considerably greater than unity in integrals such as those 
of equation (4), some special study of the integral in question is usually required 
and an effort should be made to determine some simple expression which the 
value of the integral approaches asymptotically for large values of the parameter 
or parameters. Usually this can be done by considering the physical problem in 
which the integral arose and by examining the known approximate solution for 
* rJi(px) dx. 

sinh x + x 
to be evaluated for a series of values of p covering the range 0 S p < ~. The inte- 
gral arises in a three-dimensional elasticity problem involving axially symmetric 
loads on a thick plate. By considering a limiting case of the physical problem or by 
other methods, it can be shown that 


“ tJ i(px) dx 1 
(9) A sinh = + 2 ~ 2p ; for large p. 


a limiting case of the problem. For example, suppose the integral [ 


Thus the range of values of p which must be considered is reduced to the range from 
zero to the lowest value at which equation (9) gives the result with sufficient 


accuracy. The integral in question may be expressed as a series of integrals in the 
form 


(10) f Hite 5% <n Sate 


sinha+2 421 sinh" x 


where the integrals in the right member may be evaluated by contour integration 
following the method described for a similar integral in [8]. For example, the first 
two integrals required are 

(11) [ aJi(px) dx _ 1 


- — Qalx,(xp) — 2xy(2rp) + 3x,(3xp) — ---] 
~ sinhxr p 


x2 
wo | a SS = = + Qelea(xp) +. 2a (2xp) + Bu (Bx0) + -- 


— Qn’ plko(xp) + 4xo(2ep) + Yxo(3xp) + ---] 


where xo(2) and x,(x) are modified Bessel functions of the second kind of order zero 
and one respectively. 


By evaluating the first four integrals and applying an Euler transformation to 
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the resulting alternating series in equation (10), the values of the left member of 
equation (10) for p = 4 and p = 5 were found to be as given in Table 4. The entries 
in Table 4 were completed as follows. For p = 0 and p = 1, the integral was evalu- 
ated from 


J ,(px) dx 3 3 3°55, 
ans (RSS = oe — 2 Feet 2 Teh 


For p = 2 and p = 3, a Kummer transformation was applied to equation (13) so 
as to obtain 


C 4 oo n 
(14) [ tJi(px) de = > (—1)"" [ a'Ji(on) de + Remainder Series 


o snhz+2 aa1 sinh* x 
where 


n=l 


“safe f cof 


n=l 


Remainder Series = 2{ Is - >» Lok tt Bi 


(15) 


The remainder series converges for p < 5 but the convergence is, of course, very 
slow as p approaches 5 and, even for p = 2 and p = 3, it is best to speed the con- 
vergence by applying an Euler transformation. To evaluate the series in equation 
(15), it is necessary to use Table 2, Table 3, and the equation 


|"|- 1 {fta— set mV [m 2) 
n\| (m—1)m n? n 


_ (n+1) |} 
n wae f 


which is equivalent to the first of equations (4) in reference [1]. 





(16) 
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Two Theorem Tables of Matrix Algebra 


By Gilbert C. Best 


Abstract. It is shown in this paper how it is possible, by a slight modification 
of the mode of presentation, to state many of the theorems of matrix algebra in 
table form. Due to the large number of theorems presented the tables give only the 


statement of each theorem. It is not considered feasible to also attempt to indicate 
sources and proofs. 


1. Introduction. Many theorems are of such a nature that they can be stated 
in table form. Take, for example, the theorem which reads, “If A isa non-singular 
symmetric matrix then so is A™*.” As a first step we restate this in “operator” form. 
We define the “inverse-forming” operator 0 to be such that OA = A“, where 
A is assumed non-singular. Then the theorem reads, ““The property of symmetry 
is invariant for the operator O"'.” Finally, defining “‘v” to indicate invariance, the 
theorem may be presented in an “operator vs. property” table by a v in row 0" 
and column “symmetry.” This theorem appears in Table I in row 3 column 6. This 
method of presentation has the obvious advantage of considerable compactness. 

Again, consider the theorem 


(1) (A*)* = (A™)". 


Defining a “transposing” operator 0” such that 07A = A’ this theorem can be 
rewritten in “operator” form as : 


(2) 070"°A = 0'0'A. 


The theorem can then be stated as: “The operators 0” and O™ are commuta- 
tive.’”’ In this last form the theorem in question is readily presented in an “‘operator 


vs. operator” table by defining k to indicate commutativeness. Then (1) is given 
as a k in row O’ and column O", or, equally well, by a k in row O“ and column 0”. 


2. Operators. In general, any transformation of a matrix may be considered to 


define an operator. Also, for implied given fixed coefficients a, --- a, and an 
ordered set (A;, Az, ---, Am) of matrices, the linear combination > int a;A; 


may be considered an operator. Similarly, any function of a matrix or matrices 
may be considered an operator. For example, the forming of the determinant, the 
extraction of the characteristic roots or characteristic vectors, or the product of 
two matrices, may all be considered operators. A numbered list of the operators 
used in this paper is shown below. The operators in the rows of Table 1 and the 
rows and columns of Table 2 preserve the same numbering. 


3. Properties. The idea of a property also may be generalized beyond the well- 
known and clearly defined properties such as symmetry, skew-symmetry, singu- 
larity, positive definiteness, etc. For example, the fact that a given matrix A has 
specific roots, e.g., 2, —5, 7 say, is a property, and this property is unaffected by a 


Received April 22, 1960. 
19 











GILBERT C. BEST 


List oF OPERATORS 





* Symbol Explanation 
0 I Identity Operator 
1 Or Transposer 
2 0 Conjugator, i.e., sends A = (a,;;) into A = (@;;) 
3 O- Inverse Former (if inverse exists) 
4 O04 Adjoint Former (or adjugate, see [1]) 
5 Ow kth Compound Former, see [1] 
6 PQ Equivalent Matrix Former* 
7 PPr Congruent Matrix Former* 
8 PPr Hermitely Congruent Matrix Former* 
9 PP Similar Matrix Former* 
10 a Orthogonally Similar Matrix Former* 
11 UU~ Unitarily Similar Matrix Former* 
12 > Sum Former 
13 O« Former of Product with Scalar a ~ 0 
14 LNCO Former of Linear Combination (coefficients assumed fixed)* 
15 Il Former of Product of Matrices 
16 OP Former of Matrix Raised to Power P 
17 POL Former of a Polynomial from a Matrix (coefficients as- 
sumed fixed) * 
18 FUNC Function Taker* 
19 DET Determinant Taker 
20 tr Trace Taker 
21 R Characteristic Root Taker 
22 V Characteristic Vector Taker 


* Operator matrices or coefficients are over the same field as the matrix A to 
which they are applied unless otherwise indicated. 


List OF PROPERTIES 





* 


& Property 

1 Diagonal 

2 Block Diagonal (also called “direct sum’’) 

3 Triangular. 2 types* : upper and lower 

+ Block Triangular. 2 types* : upper and lower 
5 Non-Singular 

6 Symmetric 

7 Hermitian 

8 Skew Symmetric 

9 Skew Hermitian 

10 Orthogonal 

11 Unitary 
12 Normal 

13 Diagonable, i.e., similar to a diagonal matrix 
14 Symmetric Positive Definite (real) 

15 Hermitian Positive Definite 

16 Definite. 4 types* positive definite, positive semi-definite, negative 

semi-definite, negative definite 

17 Specific Rank 

18 Specific Index 

19 Specific Invariant Factors 
20 Specific Characteristic Equation or Characteristic Roots 
21 Specific Characteristic Vectors 


22 Specific Determinant 
v in Table 1 implies that the matrix remains the same type. 
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similarity transformation. If a matrix has a specific determinant, rank, or index, 


that is also a property. A list of properties used is shown below, and the columns 
of Table 1 preserve the same numbering. 


4. Rules of Application. If an operator is applied to more than one matrix at a 
time the following rules hold: 

(1) If the operator O is applicable to a single matrix, i.e., is a matrix-to-matrix 
operator, and O is applied to the ordered set (A; , Az, --- , Am), then one obtains 
(OA, , OAz2, --- , OAm). 

(2) If the operator O is applicable to two matrices and yields one matrix and O 
is applied to the ordered set (A; , Az, --- , Am), then O is first applied to A; , A, 
yielding A,.’, then O is applied to A.’A; yielding A,’, etc. 

If a set of matrices (A; , A2, --- , Am) is said to have a property P, then each 
matrix of the set has property P. 
















































































TaBLe 1 
Operators vs. Properties 
| | | | 4 ie eye 
| 4 2/3 | | | lei#!] | 32 | § | | | | 
§ & 3 2/\z a . 3 = | Fi | 3 | Z = | 
~|2\3/2/3l2\ale)6/2 Farad i leiszis|¢ 
S19/S/G 12/8 /Slale|SlelalslaiS\s] os) | 3/181 3\¢ 
ela flalZlelelelelsia|elela\eigaialsl=l=|§ 
2\S\/2/8/8 SISSIES SiS EL SIS ISEI S| Elalals 
AlRlealalz)aie)alaiojsizjala) sam | 4/4/60 oa 
Op. 1)2) 3) 4] 5) 6) 7) 8 | 9 | 10) 11) 12) 13 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 
"| we Ber 7 We Fey a he RG Se He Ok PO he OL 
1| OT vivi ti tlol|oi|vjiolololoio)\o olelololololo| | 1 
219 vivivivjvjviviovjololololoviolololojo| -* 2 
3; 0 v olelelolelelelolelelole vijviviole lv 3 
4| OA vijvivivlvijviole viviv|olvivie | |e 4 
5) O® v v viv» v } | 5 
6| PQ ae | v | v 6 
7 PPT viv} v| | v }vjvio} 7 
8| PPT vi je v viviv v | 8 
9| PP“ | | » v Baba |v|19 
_re o|o| v |v | jvie vivijov|v| |» fl0 
11) VUU+ v| le v | v |v | lolele lo|o| |o jl 
12) > vivivl|o }vj}ol|olo| | jololb 12 
13| O* vjvlololololriolr| | jojo) ++i +o) +) 0 v 13 
14] LNCO|F|F|F|F| |Flol|olo| | | | 14 
15 vivjvjoljo| kik) vio) | | kik 15 
16| OF viviviolvio|o} lvl|v|vi|olole | v 16 
17) POL |F\F\|F\F| |Fi\ov| | |Fi lIrir v 17 
18| FUNC | F| F| F) FY | Fie} | | F| lrir lo} 18 
| | | } 
— —|-—| — -— — — 
1| 2 3| 4] 5 6 7| 8} 9) 10] 11/12] 13] 14 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 | 
Symbol Explanation 
a maintains type for even n but for n odd: neg—pos 
b maintains type, though a change from semi into non-semi is possible 
F invariant for LNCO, POL, FUNC over a field F for matrices over F 
k invariant if and only if the matrices commute 
r invariant for a, LVNCO, POL, FUNC over the real field 
t changes type 
v invariant 


al fora >0O 
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TABLE 2 
Operators vs. Operators 























& 7 |< |2 olR IRIE IA 2 8 S| a |S = 5 
ISPOHSTSISIRINIVIN IMIS |IALS IS |B o|/S |e [Sl aiele 
Op. 1 | 3} 4/5 6 | 7 s|9 10 | 11 12 | 13 | 14| 15 | 16 | 17 | 18 | 19 | 20| 21 22 | 
1| OF lr |k kik} ke] |e) | |k| [klk o|k\|k\|k|k\|ki/k}] {1 
2| 0 |e) I) k| klk | | | | [ki kl[ ri kl kil rirl kl kl ke) kl 2 
3.02 |kl ki) kik} | klk} eI | o|k k ki 13 
404 |kikik} | | | |k\ k|k o|k 4 
5| O® | k| k k | | | She k\k 5 
6| PQ | | | | | kik k 6 
dee |") 1 | | | | dy | | BIEL 3 
9} PP- | |kl ki | | | | | [el el el el el el el klk 9 
10} PP+ | k| || k| | | | | | | | Rl e ele) RL Re) Rl elk 10 
juu+ | | |kik} | | | | | [el el| kl e| kl Rl Rl Rik 11 
12 kik) | | lk) kl k| k| klk) |k | k 12 
13 kik} | | |e] ele] el kl el ae] [k k\|k| {13 
4) LNCO|k\r| | | a k) | | | k k\k 14 
15) [] ok) olo|k| | [kl k| ke} k 15 
16| OF kl k\| ki k|k| | | k| k| k| k k| {16 
17;}POL |k{r| | | | | | k| k| k | k| 17 
18} FUNC|k|r| | | | | ki k| k| k\ |18 
19} DET | k| k| k| } | taratae | k| k l19 
20) tr ki k| | * kK) k) ke) ke) | k 20 
21| R k\k|k] | | | | | |kl& k\ k\ k (21 
22) V wry a | | 22 
nail 2 EE SSR a Ss ee Se ee oe: ee ee ee rn pe eee Oe 
| | | 
1a,2)] #1 4] s 6|7 8 9 | 10} 11 | 12 | 13 14 | 15 | 16 | 17 | 18 | 19 | 20 | 21 | 22 
| ' | | | 
Symbol Explanation 
I two applications of operator equivalent to the identity operator 
0 order of product reverses, e.g. if P = operator for row in table PI][(AB) = 
[[P(BA) = P(B)*P(A) 
k operators are commutative 
r operators are commutative for LNCO, POL, FUNC over the real field 
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Improved Formulas for Complete and Partial 
Summation of Certain Series 


By Herbert E. Salzer and Genevieve M. Kimbro 


Abstract. In two previous articles one of the authors gave formulas, with nu- 
merous examples, for summing a series either to infinity (complete) or up to a cer- 
tain number n of terms (partial) by considering the sum of the first 7 terms S; , or 
some suitable modification S; , closely related to S; , asa polynomial in1/j. Either 
S.. or S, was found by m-point Lagrangian extrapolation from S;,, Sj1, ---, 
Sjo-m4i1 to 1/7 = 0 or 1/j = 1/n respectively. This present paper furnishes more 
accurate m-point formulas for sums (or sequences) S; which behave as even fune- 
tions of 1/j. Tables of Lagrangian extrapolation coefficients in the variable 1/7 
are given for: complete summation, m = 2(1)7, jo = 10, exactly and 20D, m = 11, 
jo = 20, 30D; partial summation, m = 7, jo = 10, n = 11(1)25(5)100, 200, 500, 
1000, exactly. Applications are made to calculating z or the semi-perimeters of 
many-sided regular polygons, Euler’s constant, 


pie. 1 “era . 
1+ Sgatw@a i eteo for j= @ 





(Catalan’s constant), calculation of a definite integral as the limit of a suitably 
chosen sequence, determining later zeros of J,(2) from earlier zeros for suitable », 
etc. A useful device in many cases involving sums of odd functions, is to replace S; 
by a trapezoidal-type S; which is seen, from the Euler-Maclaurin formula, to be 
formally a series in 1/7’. In almost every example, comparison with the earlier 
(1/j)-extrapolation showed these present formulas, for 7 points, to improve re- 
sults by anywhere from around 4 to 9 places. 


1. Introduction. In two earlier papers, [1, 2], one of the authors gave tables for 
both complete summation (all terms, to infinity) and partial summation (up to a 
certain number of terms) of certain kinds of slowly convergent series. In the case 
of partial summation, divergent series were also included, provided that a suitable 
auxiliary series could be found of the desired slowly convergent type and simply re- 
lated to the original divergent series. The essential idea in both cases is to regard 
the sequence S; , the sum of the first 7 terms of the series, as the values for x = 1/7 
of an interpolable function S(z) to which the slight extrapolation from specified 
S;,toj = © (x = 0) ortoj = k (x = 1/k), k > jo where S,j, is the last speci- 
fied S; , yields good accuracy. The approximating formula for S(2) was an m-point 
Lagrange polynomial of the (m — 1)th degree in x which at x = 1/7 assumes the 
prescribed value S; , for the last m values of j ending at jo = 5, 10, 15 or 20, from 
which we extrapolated to either j = © (x = 0) orj = k > jo(x = 1/k). Numerous 
examples which yielded surprisingly high accuracy for a variety of sequences S; in 
both complete and incomplete cases, attested to the wide applicability of consider- 
ing S; a smooth function of 1/j, even when we were in complete ignorance as to the 
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actual analytic expression for S; = S(1/j) or of a theoretical justification for con- 
sidering S(1/j) as an approximate polynomial in 1/). 

However, a still further improvement in m-point formulas for both complete 
and partial summation is applicable to a wide class of sequences where S; = S(1/7) 
behaves like an even function of 1/j. Thus by taking x” as argument instead of z, in 
an m-point Lagrangian extrapolation formula for x = 0 (complete summation) or 
a value close to 0 (partial summation) based upon those same final m values of S; , 
we should get accuracy equivalent to (2m — 1)th degree instead of (m — 1)th 
degree. As will be seen from the illustrations below, the resulting improvement is 
often quite impressive. 

There is no hard and fast classification of all the varied problems to which these 
newer formulas are applicable. The reason is that even if a problem does not seem 
offhand to involve a sequence of that even-function type, often with a very slight 
transformation, regrouping, or alteration, one sees that it really is amenable to this 
more accurate treatment. 

Of course, every sequence to which these improved extrapolation formulas for 
arguments 1/j’ are particularly applicable can also be handled by the earlier formulas 
employing arguments 1/j, because any polynomial in z’ is also a polynomial in z, 
but with considerably less accuracy for the same number m of points and the same 
last 7 = jo. But the converse is not true—we cannot in general expect these newer 
summation formulas to work well when applied indiscriminately to sequences where 
the earlier method may give very high accuracy. One way of realizing this is to 
think of the non-constant part of a well-behaved function of z near t = 0 being 
approximated by Cz. Extrapolation employing x° = y as the variable, near x = 0, is 
like extrapolation for ~/y based upon a polynomial approximation in the variable y. 
But, as anybody who has attempted to interpolate in a table of square roots near 
zero has found out, +/y, although continuous at y = 0, possesses a singularity due 
to an infinite derivative. 


2. Other Related Articles. The idea of the extrapolation to x = 0 for argument 
y = x has been employed for just the linear case in the well-known “h’-extrapola- 
tion process”, or “deferred approach to the limit’’, which has been extensively 
treated in the literature on the numerical solution of differential equations (first 
introduced by L. F. Richardson [3, 4]). The argument x or h corresponds to two 
conveniently small values of a mesh-length, say h; and he . Richardson’s process has 
been generalized to higher powers beyond h’ by several writers, notably G. Blanch, 
[5] and H. C. Bolton and H. I. Scoins [6]. However, the only reference that was en- 
countered by the writer which was concerned with problems where the approxima- 
tion might be considered as a purely even function of h having more than a single 
term, has been M. G. Salvadori [7]. Besides some sets of 2-point coefficients for h’- 
and h*-extrapolation, Salvadori tabulates 3-point coefficients for (h’, h*)- and (h', 
h*)- extrapolation, and 4-point coefficients for (h’, h*, h®)- and (h‘, h°, h®)-extrapola- 
tion. The values of A are in the form 1/n; , where n; are sets of small integers rang- 
ing from 2 to 8. Salvadori gives applications to numerical differentiation and inte- 
gration, as well as to some boundary value problems and characteristic value 
problems. 
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3. Formulas for Complete Summation. In choosing a jp suitable for most com- 
plete summation purposes, we wish to obtain a substantial increase in accuracy 
over the use of the earlier formulas in [1], which has already been proved to be very 
accurate, without having coefficients that might be too cumbersome. It is also de- 
sirable to give exact values rather than decimal values, because in highly accurate 
formulas the theoretical or truncation error might be considerably smaller than the 
computing error arising from the use of rounded decimal entries. But we must also 
take account of the fact that the fixed points 1/j° in place of the older 1/7 makes 
the exact fractional form of the extrapolation coefficients have around twice as 
many digits in both numerator and denominator, which adds considerably to the 
amount of time to do an example. 

In the present paper it seems that a very convenient choice is jo = 10, for all 
cases ranging from the 2-point through the 7-point. In other words we give formulas 
for linear through sextic Lagrangian extrapolation formulas for functions of the 
variable y = 2° taken at x = 1/j, or arguments y = 1/7 at j = 10,9, --- ,10— 
m + 1 for m = 2(1)7. This is equivalent to quadratic through twelfth degree 
accuracy for even functions in z = 1/j. The extrapolation formula to obtain the 


complete sum S from the partial sums Sy, Sy, --- , Sy—m4: is the very simple 
m—1 
(1) S~ > A{hosSr0-s : 


The coefficients A{f}o_; are given in Table 1 in exact fractional form BST}o_;/D{>’, 
so that (1) may be most conveniently employed as 


m—1 
(2) S ~ (1/D{) > Bir te :Sie - 


In no case through m = 7, does D{f? have more than ten digits exclusive of final 
zeros, which is convenient in the division. The values of A{f}o_; are given also to 20 
decimals in Table 2. 

Although the 7-point formulas for jo = 10 are very accurate, as will be apparent 
from the examples below, we give also in Table 3 for possible use in some kind of 
isolated key calculation where extreme accuracy is sought, even at the expense of 
considerable computing labor, the coefficients in the 11-point formula, ending at 
So , given exactly, to be employed in 


10 
(3) S~ (1/Die”) 2 Biiis0—<Sa0—: ‘ 


Formula (3) is exact for any even polynomial in x = 1/j up to the 20th degree. To 
avoid too much non-essential numerical work, no illustrations were given of the 
use of Table 3, since the resulting accuracy is so high by comparison with the results 
of using Table 1 or 2, that an excessively large number of significant digits is needed 
to reveal its full extent. But Table 3 should be kept in reserve for a summation 
problem requiring unusual precision. 

The formula for A$”,_; is obtained rather simply from the well-known defini- 
tion of the m-point Lagrangian interpolation coefficients where we have fixed points 
1/jc, 1/(jo — 1)*, --+ , 1/(jo — m + 1)’ and set the variable y = 2° = 1/)° equal 
to 0 to correspond toj = ~. 
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TABLE 1 





(m) (m) (m) 
Aio10~i = Bioi0-:/Dio 




















m= 2 | m= 5 
(2) . 
Bis = —81 Bis = 2034 43488 
Paw = 160 Bi, = —23001 55599 
= 2 | Bis = 82879 44704 
| BS, = —1 17517 54833 
: ams | Bes, = 56875 00000 
Bi, = 19456 o — 1269 77760 
Bi, = —59049 | 
Df = 2907 A 
| {3 =  —75703 12500 
m=A4 BS, = 17 57751 73632 
os | Bioy = —123 97838 67861 
| TB mon Bi, = 359 05926 50456 
Biss = 308 45088 | Bi, = —448 74915 24087 
= = | Bi? 10= — 200 20000 00000 
ae - ‘a ror | D® = 3 35221 28640 
a 0 Jo 
m= 7 
BY, = 54190 40768 
Bi; = —31 54296 87500 
BS. = 474 59296 88064 
Bi?, = —2761 33679 65995 
BS, = 7181 18531 89120 
Bi, = —8388 15723 34857 
BSo10 = 3575 00000 00000 
DD? = 50 28319 29600 
” ( 4)?" Ba i 
(4) AS-i _ =—1 » eee 





’ [jo — k)® — (jo — i)’ 


k= 


o 


where in |’, k = 7 is omitted. 


4. Illustrations of Complete Summation. 

A. Example 1. Considering the circle as the limiting case of inscribed regular 
polygons of j sides, as j — ©, the quantity z is the limit of the semi-perimeter, 
j sin a, where a = 180°/j = x/j, as j — «.* Now the approximation S; = j sina = 

* Although this example affords a splendid illustration of the improvement of (1/j*)- 
extrapolation over (1/j)-extrapolation, it suffers from the aesthetic defect of having the value 


of x occurring implicitly in every S; in the various powers of a needed to compute sin a. In 
other words, there is definitely something ‘‘circular’’ in this example. 
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TABLE 2 
A{pi0-« in Decimal Form 
(2) ‘ 
Ajoy = —4.26315 78947 36842 10526) , 6 emit v 
Aftie = 5.26315 78047 36842 10526410) — 0. See Soeee SENSS Gnee 
Ai, = 5.24355 64435 56443 55644 
awe ‘Alor = —36.98404 36201 18987 76605 
(6) 
As = 6.69281 04575 16339 seg2gi4%9s = 107. 11111 74954 31220 92142 
A{, = —20.31269 34984 52012 38390\4ie9 — —133.86654 44631 80008 84564 
Aro = 14.6198 30409 35672 5146241010 = 99. 72174 44482 66636 08913 
nee m=4 m= 7 
A a 
Aig; = — 4.80592 32026 14379 084971 4°, = 0.01077 70418 88152 99926 
Afts =  28.55599 12854 03050 10893) 40) _ 
. ‘Aios = —0.62730 63998 75844 32029 
Af}, = —51.41650 54179 56656 34675| 40) _ : : “ 
a 5 Tene eens 3oo78\4 i = 9.43840 15984 01598 40160 
yes ont an e24S/Afor = —54.91570 11329 03951 53140 
wienit fos =  142.81482 33272 41627 89522 
Ajos = —166.81830 92541 1 5 
Aives = 1.60219 78021 97802 19780) 47" _ o pn a prepa paired pte 
A{p, = —18.11463 36098 54198 08949; —" 
- Aft; =  65.27083 72237 78400 24899) 
AS), = —92.54970 97523 21981 42415 
Afpio =  44.79130 83361 99977 06685| 
TABLE 3 
Asozo-i = Bio2o-:/D5o” 
Bye = 74096 20000 00000 00000 00000 
Byun = —35 37615 48335 31708 54782 90644 
S012 = 649 45974 08685 61313 24915 22048 
BS: = —6200 60319 26092 91850 74192 35023 
wee BY. = 34801 60376 52150 35629 23772 47744 
Sots = —1 21941 46052 37160 60638 42773 43750 
B32, = 2 73659 70208 28851 47761 53823 64160 
Soir = —3 92511 27655 98026 11495 97941 97770 
Bos = 3 47343 22454 05086 94470 03616 05120 
Boy = —1 72481 59320 99496 29170 21217 51885 
BSi2o = 36718 51008 00000 00000 00000 00000 
DS” = (32124 40751)(38698 35264)(23 58125) 
suler a 2 93153 05663 14310 15219 20000 _ 
eter, ; 
a= S(1/j) is seen to be an even function of 1/7 which equals x for 1/7 = 0. Therefore 
1/52)- we expect an m-point Lagrange polynomial approximation for variable 1/ '? to be 
ys Sl considerably more accurate than a polynomial in 1/j. Following are the values of 
a. In the semi-perimeters j sin a to 25D, which were obtained from a table of sin a to 
30D originally published by Herrmann [8]. For 7 = 4(1)6, 9, 10, sin a was copied 
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from Herrmann’s table, and for 7 = 7, 8, sin a was computed by Taylor’s theorem 
employing Herrmann’s entries as key values: 





». 


Sj: Semi-perimeter = j sin a 





2.82842 71247 46190 09760 33774 
2.93892 62614 62365 64584 35298 
3.00000 00000 00000 00000 00000 
3.03718 61738 22906 84333 03783 
3.06146 74589 20718 17382 76799 
3.07818 12899 31018 59739 68965 
3.09016 99437 49474 24102 29342 





SOON 


— 





In the above values of S; , as well as S; given in the other examples, the accuracy 
of the last few places, although highly probable, is still not absolutely guaranteed. 
However, in every example the values of S; are certainly correct up to the number 
of places needed to guarantee that the “computational error” in the final answer 
(which is due to the error in the S; multiplied by the extrapolation coefficients 
A{w) is appreciably less than the deviation of the answer from the true value. This 
latter “truncating error” is thus made to stand out clearly, and it indicates the 
theoretical accuracy of the extrapolation formula, regardless of the number of places 
carried in the work. In practice we do not often know at the outset of an example 
just how many places are needed in the S; to assure us that the computing error will 
be dominated by the truncating error. Sometimes when the theoretical accuracy 
turns out to be unexpectedly fine, the example must be done again, carrying more 
places, to prevent the computing error from obscuring the truncating error. 

The results of the extrapolations employing (1) or (2), for m = 7, gave for z, 
(whose true value to 20D is 3.14159 26535 89793 23846), the answer 3.14159 
26535 89793 179 ... which is correct to within a unit in the 16th decimal. The 
extent of the improvement over the earlier (1/j)-extrapolation formulas is appar- 
ent from the result of 3.14159 280 ... obtained by the corresponding 7-point (1/7)- 
extrapolation coefficients, which deviates from + by 14 units in the 7th decimal. In 
other words, the error in the use of this newer formula is only around 0.4-10~° of 
that in the older one. The greater power of this newer method in this present exam- 
ple may be further illustrated even for m = 4, where (1/j’)-extrapolation yields 
3.14159 2650 ..., or accuracy to around 3 of a unit in the 8th decimal, whereas 
the corresponding 4-point (1/7)-extrapolation formula gives no better than 3.1411 

. , which is off by $ of a unit in the 3rd decimal. In fact, the answer even by 
2-point (1/7’)-extrapolation, namely 3.1413... , is still better than the above 
3.1411.... 

It is interesting to note that the use of (1/j’)-extrapolation on the semi-perim- 
eters gives this great improvement only for the inscribed polygons, and it will not 
work well for the circumscribed polygons, upon which it was also tried. A reason 
that would lead us to expect poor extrapolation results, even though the corre- 
sponding semi-perimeter j tan a is still an even function of 1/j, is that the series 
for tan a converges poorly by comparison with sin a. Thus for a = 2/4, occurring 
in S; = S,, the remainder after the term involving the sixth power of 1/)’, is con- 
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siderably greater for j tan a, so that the use of (1) or (2) for m = 7 is not nearly 
so good as for j sin a. 

B. Example 2. The sequence for Euler’s constant 


j 
— ime (1/r) — log i = 0.57721 56649 01532 86061 to 20D 
r=1 ) 


has been treated earlier by (1/j)-extrapolation ({1], p. 358). Applying (1) or (2), 
for seven points, directly to S; = }°?.1 (1/r) — log j yields the very inaccurate 
0.593, the reason being that S; does not behave like an even function of 1/j. The 
older (1/j)-extrapolation formulas, employing jo = 10, gave 0.57721 41 ... and 
0.57721 56695 ... by the 4- and 7-point formulas with respective errors of around 
13-10°° and 4-10~*. To improve upon these results we must modify our S; se- 
quence into an even function of 1/j having the same limit y. This is easily accom- 
plished by replacing the last 1/r in the summation, namely 1/r = 1/j, by half its 
value, or 1/27. At first sight there is an apparent motivation in that the new finite 
summation is suggestive (at one end anyhow) of the more accurate trapezoidal 


i 
rather than rectangular approximation to the integral i (1/r) dr. This trapezoidal 
1 


motivation happens to lead to the correct choice in this present example, but in 
general it does not yield a sequence that is even in (1/7). The true motivation lies 
in the Euler-Maclaurin summation formula applied to log j. The general formula 
is expressible as 


a+jw 
a f(x)dz = (Bh thtfet---+fiat +i) > (fi — fr’) 


3 5 
) yt yer w 7 v 
sop Si — fe") — aera Ui? — fe") + --- 19). 
Now (5) does not denote a complete equality, since the Euler-Maclaurin formula 
is an asymptotic expression that is given with a remainder term. Employing (5) 
heuristically for w = 1,a = 1 and f(x) = 1/z, the right member of (5), exclusive 
of the (4fo + fi + --- + fj. + 4f;) and an undisclosed remainder term, is an even 
function of 1/(j7 + 1), from which, replacing 7 by 7 — 1, 

. We FE 1 1 


is an even function of 1/j, so that the same is true of the sequence 


(5) 
yy 


, = 1 J eas i i _ j 
s/a(1+h+it + +2) log j 


whose limit, as 7 — ©, is also equal to +.* 
Since the older m-point (1/7)-extrapolation formula is linear in S; (or S,/) and 





* The reader is cautioned that the above heuristic demonstration is not to be understood 
as a proof that we have a convergent infinite series in (1/j*) from which we can “‘prove’”’ that 
the ‘“‘constant’”’ term in S;’ is y by taking the limit as 7 — ©. The fallacy there would be in 
that there is no ‘“‘constant” term because the fo , fo’, fo, --- terms in (5) yield for f(z) = 1/z 
a divergent sequence. Actually S;’ is defined only up to any fixed order derivative, say f $ a. 
and it then consists of terms in 1/j?, constant terms and an integral formula for the remainder. 
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yields exactly zero for any polynomial in 1/j having no constant term, up to the 
(m — 1)th degree, the above-mentioned 4- and 7-point results will not be changed 
by use of S,’ instead of S;. But the improvement is very noticeable when S,’ is 
employed with (1/j’)-extrapolation. Following are the terms in the modified 
sequence S,’ to 20D: 





ee 
j S;'= 2 > + 3; — ogi 


ru * 





0.57203 89722 13442 71450 
0.57389 54208 99232 95873 
0.57490 71974 38611 66585 
0.57551 84223 73258 12347 
0.57591 56011 77306 92889 
0.57618 81210 76479 02991 
0.57638 31609 74208 28424 





SOON OS 


_ 





The use of the 7-point formula in (1) or (2), where jo = 10, upon S,’, gave an an- 
swer of 0.57721 56649 0143 ... which is correct to a unit in the 13th decimal 
(i.e., 5 places more than (1/j)-extrapolation). Use of just the 4-point formula in 
(1) or (2) gave an answer as good as 0.57721 56647 5 ... which is correct to 
within 13 units in the 10th decimal (i.e., 4 places more than (1/7)-extrapolation). 

C. Example 3. A different type of sequence is encountered in the evaluation of 

1 
oa 1 

the definite integral [ ee 
59945 30942. One obvious sequence to consider is S; which is formed by dividing 
the interval (0, 1) into 7 equally spaced intervals and letting S; be the sum of the 
rectangles of height 1/[1 + (r — 1)/j] and width 1/j, for r = 1(1)j, but that fails 
to behave as an even function of 1/7. However, the trapezoidal rule, or 


1/1 1 1 1 1 
sy =1 (5+ . tee tii +3), 
"GN 14g 14+2/ 1+(@j-)/j 4 
according to the Euler-Maclaurin formula (5), where now w = 1/j, a = 0, and 
both f;” is fixed as well as fo”, being at the endpoints 1 and 0, is seen to have a 
truncating error that is formally a series in 1/;°. The values of S,’, in either exact 
form, or to 20D, are as follows: 





dx = log 2, whose value to 20D is 0.69314 71805 














‘ 7 ae 1 

: sr=3(i+ 2 ain +t) 

4 1171/1680 = 0.69702 38095 23809 52381 
5 1753/2520 = 0.69563 49206 34920 63492 
6 9631/13860 = 0.69487 73448 77344 87734 
7 2 50241/3 60360 = 0.69441 94694 19469 41947 
8 2 00107/2 88288 = 0.69412 18503 71850 37185 
9 5 66803/8 16816 = 0.69391 76020 05837 29995 
10 1615 04821/2327 92560 = 0.69377 14031 75427 94323 





The 4- and 7-point (1/7)-extrapolation, jo = 10, gave values of 0.69314 86 ... and 
0.69314 7176 ... , correct to 1} units in the 6th decimal and 4 unit in the 8th 
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decimal respectively. The (1/j*)-extrapolation was performed for every m-point 
formula from m = 2 through m = 7, with the following results: 








m | value of S | deviation | = | value of S deviation 
| — | 

2 | 0.69314 81... 10-* 5 | 0.69314 71805 67... 10-" 

3 0.69314 7188... 10° 6 | 0.69314 71805 6054... #-10-" 

4 0.69314 71807 1... 13-10-” | 7 | 0.69314 71805 60046... 10-* 





The improvement over (1/j)-extrapolation in the 4- and 7-point results is by four 
and five places respectively. 
D. Example 4. A somewhat more sophisticated application of (1/j)-extrapola- 
tion is in the summation of the series for Catalan’s constant, or 
1 1 1 


T: = 1 -gtep-azt:: 

H. T. Davis [10] gives a full discussion of Catalan’s constant, including an account 
of the earlier work of J. W. L. Glaisher, and he also reprints Glaisher’s 32-decimal 
value of T; = 0.91596 55941 77219 01505 46035 14932 38. Since the series for 
T:2 is absolutely convergent, it may be regrouped as 


1 1 1 1 1 1 
. 1+(-5+5)+(-p+5) + +(-@ yt &+ + - 


The general term u, , r > 0, of T2 is equal to - ae > _ 7? which is an odd function 


of r or 1/r. Thus, as in the preceding example, employing (5) with w = a = 1, the 
modified sum 


j-1 ¢ >° 
S,’ = S; = 4 u; = a eg pe l6r = ¢ . ww .) 


= (16r? — 1) (1672 — 1)? 





is again seen to be formally an even function of 1/j, having the same limit S which 
is approached by S,.* The values of S,’ to 20D are as follows: 


: “ys i 10 De 16) 
laa ray (16r?— 1)? = 2 \ (167? — 12 











a 0.91798 69831 73330 85103 
5 0.91724 36100 54163 02747 
6 0.91684 71757 66868 06945 
7 0.91661 06554 47552 03321 
8 0.91645 81601 71966 79489 
9 0.91635 40724 61230 31205 
10 0.91627 98501 91732 37910 








* Although in Example 4 we know the explicit formula for / f(z) dx, we may expect this 
1 


principle to be applicable also in cases where i f(x) dx, f(z) odd, or for that matter also f;<”? 
1 


for odd p, is not known in closed form, and where S;’ may still be regarded formally as a series 
in 1/}?. 
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Use of the 7-point (1/j)-extrapolation, j, = 10, upon either S; or S,’, while not 
identical in accuracy, because now the difference of $(16j/(16j’ — 1)”) isno longer an 
exact polynomial in 1/j, gave results very close to each other, namely 0.91596 
55973 ... and 0.91596 55980 ... with respective deviations of 4-10~* and 3-10~°. 
The use of (1/7’)-extrapolation, i.e., (1) or (2), for m = 7, while giving the poorer 
answer of 0.91596 74 ... with a deviation of 2-10°° in working with the S; se- 
quence (as was to be expected), gave upon working with the S,’ sequence the highly 
accurate 0.91596 55941 7714 ... , which is correct to $-10-", showing a gain in 
accuracy of around 5 places. 


5. Formulas for Partial Summation. Given the first ten terms of a sequence 
S; which behaves as an even function of 1/j, we might wish to find by (1/’)-extrap- 
olation S,,” > 10, instead of going to the limit as 7 — «. The purpose of this 
section is to improve what was accomplished in [2] where just (1/j)-extrapolation 
was employed. The m-point formula for S, which occurs usually as a sum of the 
form > 7-0 u, , is obtained by setting x = 1/n’ in the Lagrange interpolation coeffi- 
cients whose fixed points are 1/jo, 1/(jo — 1), --- , 1/(jo — m + 1)’. In the 
present instance, in order to avoid too much tabulation, since now besides jo and 
m, n is also a variable, being no longer just ©, we consider a choice of j) and m 
which shall be suitable for most problems and which shall give a substantial in- 
crease in accuracy over the (1/7)-extrapolation formulas previously given which 
were based upon jo = 10 and m = 7 [2]. Thus it is natural to take 7) = 10 and 
m = 7 for these present formulas also. The argument n = 11(1)25(5)100, 200, 500, 
1000, and all coefficients are given exactly. This range of n is not quite so extensive 
as in the previous paper because the arguments 1/)’ in place of 1/j, 7 = 4, 5, ---, 
10, n, increase the labor in computing the exact forms, which also have consider- 
ably greater bulk in figures. To find S, = S(n), we employ the extrapolation 
formula in the following form: 


10 
(6) S(n) = 2) Aj(n)S;. 
Every set of coefficients A ;(”) is given in the exact fractional form of C;(n)/D(n) 


where D(n) is the least common denominator for each n. Thus it may help the 
computer to have 


(7) S(n) = (1/D(n)) 2, Ci(n)S;. 


In (6) and (7) the jo = 10 is understood as well as m = 7. When also n is under- 
stood, we may employ for (7) the somewhat more concise 


10 
(7’) Sn = (1/D) 2, CiS;. 
jJ= 
In (7), or (7’), the D(n), or D, is given also in the form of factors having no more 


than 10 digits, exclusive of terminal 0’s, to facilitate the divisions on a ten-bank 
desk calculator. The C;(n) and D(n) are shown in Table 4. 





TABLE 4 A;(n) = C,(n)/D(n) 
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The coefficients A ;(n) = C;(n)/D(n) were calculated directly from the formula 
10 
a []’ (n’ — F’) 
(8) Aj(n) = 23; 
IT’ @ - #) 
k=4 
where k = j is absent from ||’. Both the calculation of A ;(n) and the determina- 
tion of D(n) was facilitated by expressing each of the factors in the right member of 
(8) in terms of powers of primes. 
To facilitate the use of (8) for desired values of m other than in this present 
table, we notice that we may express A ;(n) as 


10 


, Ss - -» 
(9) A;(n) = B; - ut (n K) , where 
n2 
- 
a 
(10) B; Il’ (7? =f k*) 


k=4 


is independent of n. The exact, as well as 30 decimal, values of the fundamental 
quantities B; are given in the following Schedule 1. 


ScHEDULE 1 





10 
j B; = j# I’ UG? — k) 
k=4 











+ “anes = 0.01077 70418 88152 $9926 41103 75221 
5 | = ee = -—0.62730 63998 75844 32028 87647 33209 
6 | 2 ie = 9.43840 15984 01598 40159 84015 98402 
7 | = "2530 0080 = —54.91570 11329 03951 oF 25117 94669 
8 | 78 79605 ~ 142.81482 33272 41627 89522 26664 64497 
9 EE a On = —166.81830 92541 16626 40764 80794 74705 
10 | ieee = 71.09731 48193 65042 96325 41775 64463 





6. Illustrations of Partial Summation. 

A. Example 5. Suppose that in Example 1 above, instead of passing to the limit 
as j7 — © to obtain zx, we wished to calculate Sx , or the semi-perimeter of a 20- 
sided regular polygon from the semi-perimeters of the 4- through 10-sided regular 
polygons. We have Sx = 20 sin 9°, whose value to 20D is 3.12868 93008 04617 
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38020. Using the same values of S; as in Example 1, we find by the earlier method 
of (1/j)-extrapolation [2] Sx. = 3.12868 93076 ... which is correct to around a 
unit in the 8th decimal. But use of the present tables for (1/7*)-extrapolation in (6) 
or (7), for n = 20, yields the highly accurate Sx = 3.12868 93008 04617 359... , 
correct to about 2 units in the 17th decimal, showing a gain of around 9 places. 

B. Example 6. As an illustration of a different type of problem that does not 
correspond to one in complete summation, consider the case where from the first 
few known zeros of some higher mathematical function, we wish to obtain the value 
of some later zero, say the nth. As will be seen below, there are circumstances when 
it is preferable to choose as the sequence S;, 7 < jo, from which to extrapolate, 
some suitable even function of 1/7 which may not be a function of the jth root, and 
yet from S;,7 > jo, the jth root, is readily obtainable. 

Consider the problem of finding the later zeros of the spherical Bessel functions 
J2m44(z) from either tabulated earlier zeros or some other suitable function of m. 
In the general asymptotic formula for z,“”, the nth zero of J,(z) cos a — Y,(z) sin a, 
namely, 


(a a “Se 4y* — 1 

ng =(n +5 i) “~ Bin +  — Dr — a] 
_ _(4e* — 1)(280* = 31) 
384{(n + 4v — })x — a}? 





(11) 





ox --+ [11], 


set a = 0 and v = 2m + 3. Then from (11) it is apparent that 
(12) Sram = (n + m)[zinty — (n + m)x] 


has a formal expansion in even powers of 1/(n + m), which could serve as the basis 
of an extrapolation formula. 

However, after searching for ready-made tables of Zn , hone were found capa- 
ble of testing the full potentialities of Table 4. To avoid extra labor, we shall first 
illustrate this principle of (1/»’)-extrapolation with a smaller example limited to 
the available published 6D values of z$7? as far as n = 6 [12]. The problem is to 
calculate z§7? for n = 6, whose published value is 24.727566, from the four preced- 
ing values of z$7> = 11.704907, 2{72 = 15.039665, 2$/2 = 18.301256 and 2{)) = 
21.525418. In other words, since m = 2, the problem is to find S, from S,, S; , S¢ 
and S; , from which 2$?} is found from (12). From (8), with Ils. replaced by 
[]‘i-:, we find A4(8) = —a¥te, As(8) = rB-44$H, Ac(8) = —H4$$4¢ and A,(8) = 
$23543 from which S; = > j-4 A;(8)S; = —3.241393. Finally, from (12), z$/> is 
found to be 24.727567, which deviates by only 10° from the published value.* 
Comparing with (1/v)-extrapolation based upon those same values of S, — S;, 
and where A,(8) = —4, As(8) = 7$%, Ac(8) = — $y, A7(8) = $54, we find S,; = 
— 3.241225, from which z\)> is found to be 24.727588, which deviates by 0.000022 
from the published value. 





* Since we started with 6D values, it is not possible to estimate from this example the 
possibly higher theoretical accuracy in (1/»*)-extrapolation, which is just the truncation error 
when the example is done with a sufficiently large number of places both initially and in the 
course of the work. 
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For a similar example employing Table 4, and revealing the full accuracy of (6) 
or (7), we choose a modification of S,4m, Say Snim, Where 


(13) Sram = (n + m)[Zimiy — (nm + m)z], 


and where now 25,+4 , instead of being the nth zero of J2m4;(x), is defined as a pre- 
assigned number of terms of the right member of (11) (for a = 0, v = 2m + 3) 
which is the same for every n. For the lowest values of n, there will be considerable 
deviation between the true value of the root z}n; and the function 2522.4 which is 
(n + m)x + an exact odd polynomial in 1/(n + m), making 8S,,,.,, an exact even 
polynomial in 1/(n + m). But at the inconvenience of having to compute S,,,, for 
the initial values of n, we may employ (6) or (7) to extrapolate for S,,,, for some 
larger n to get 252.4 which will agree with the true value of the root 24%, to very 
high accuracy. Taking (11) as far out as 1/{(m + 34» — 4)x — a}’, we have for 
a = 0,» = 2m + Sand p= 4” = (4m + 1)’, 


5 u—1_ (wu —1)(7 — 31) _ (u— 1)(83y" — 982u + 3779) 











Snim = — Be = 3: Bat(n + m)? “15 + 2%5(n + m)* 
_(u- 1)(6949u* — 153855” + 15857434 — 6277237) 
(14) 105 - 2% - x7 -(n + m)s 


(wu — 1)(70197u* — 24 79316u* + 480 10494,” 
___ — 5120 625484 + 20921 63573) * 
40320 - 2" - x°(n + m)8 
4) 


Suppose that the problem is to calculate the 14th zero of Js/2(z) or z$)3’. Then 
m = 1, and we should want to find 8. using Table 4 upon S, — Sy, after which 
we obtain 2{)2 from (13).** From (14) and then (13), 2{)2’ which is equal to 2$)2 
to around 14D, is found to be 47.06014 16127 6054. A quick examination of the 
ratios of successive terms in (14) indicates without having to compute the 
1/(n + m)” term that, to 14D, z{}? is actually 47.06014 16127 6053. Following 


are the calculated values of 8; , for 7 = n + 1 = 4(1)10, to 16D (last figure ap- 
proximate) : 














fa 5; 

4 —0.97371 85140 72535 8 
5 | —0.96680 12788 75286 8 
6 | —0.96311 73960 26803 8 
7 | —0.96092 04667 12625 8 
8 | —0.95950 42113 72512 2 
9 | —0.95853 75688 13022 8 
10 | —0.95784 82845 01448 5 

Employing the older (1/j)-extrapolation, we find S;,; = —0.95622 28677 507 ... 


and from (13), 2)2 = 47.06014 16126 6... which agrees with the true value of 


* The coefficients through 1/(n + m)*® are from Watson [11], and the coefficient of 1/(n + m)® 
is from Bickley and Miller [13]. 

** This particular problem could, of course, be set up equally efficiently computationwise 
by writing S,4: = ao + ai/(n + 1)? + --- + a,/(n + 1)”, where a; is independent of n. But 
this present method works as long as we know somehow the values of S; . 
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zsj2. to a unit in the 10th decimal (12th significant figure). But the (1/7*)-extrapo- 
lation yields $5 = —0.95622 28662 9517 ... and from (13), 2}2) = 47.06014 
16127 6055 ... which almost agrees with the true value of z;}}’ to 14 decimals (16 
significant figures). 
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Approximation of Curves by Line Segments 


By Henry Stone 


1. Introduction. Most methods of linear and nonlinear programming developed 
up to this time were designed to find the maximum or minimum value of a linear 
or nonlinear function inside a region bounded by hyperplanes. In applications of 
these methods it is often found that the formulation leads to constraints which are 
nonlinear. For example, in gasoline blending problems the relationship between 
lead concentration and octane number of a blend (the lead-susceptibility curve) 
has been found to be exponential in form. One method of dealing with such prob- 
lems has been to approximate the given curves by a series of broken line segments 
[1]. The usual method of finding such an approximation has been, by visual examina- 
tion of the graph, first to decide on the number of line segments, second, to select 
the intervals over which each line segment is to apply, and third, to draw in what 
appears to be good linear approximations to the curves over the selected intervals. 
Since the problem of obtaining a best fit of broken line segments to a curve does not 
seem to have been previously investigated, it is the purpose of this paper to formu- 
late the problem, give a closed form solution when the given function is quadratic, 
show a general numerical method of solution, and apply this numerical procedure 
to the lead susceptibility curve. 

In the case when the function is quadratic, an interesting and simple result was 
obtained. It was found that in fitting N lines over some interval (a, b), that the 
points at which one line segment was discontinued and the next line segment started 
were equally spaced over the interval. This result allowed the equations for the 
lines to be expressed in a very simple form. 


2. Formulation of the Problem. Given a known nonlinear relationship 
(1) y = f(z), U Sr Suy. 


It is desired to obtain an approximation of the form 


a, + biz, mSztSiu 
(2) - jn ee ee 
ie + bya, Un ST S Uy 
which is best in the least squares sense. If the points wu , v2, --- , Uw are specified 


in advance, the problem reduces to the simple case where the best fit over each 
interval is obtained separately by the usual least squares procedure. The results 
developed below will give the solution to this problem as a special case of the 
general solution. The more general case, treated in this paper, occurs when the end 
points of the line segments uw , U2, --- , Uv-1 are not specified in advance, and when 
only the number of line segments, N, is specified in addition to the relationship 
given by Equation (1). 


Received December 28, 1959. 
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Following the least squares formulation, we wish to obtain the values of the 
variables which minimize 


F(a, a, +++, dy, bi, be, -++, by, i, Ue, ***, Uys) 
(2) N uj 
=D We) - a — 62) ae. 

Jmk OUj-1 
The method to be used is the usual normal equation technique of finding the point 
at which the partial derivatives of the function (2) are zero. This procedure leads 


to the following system of 3N — 1 equations to be solved for the 3N — 1 variables 
M&,°** , Unw-i: 


a =0=-2f" 


(3a) da, I”) dz + 2a,(uy — up-1) + by(uy — uy.) 
p =1,2,---,N 
(3b) s, =0= -2 ‘ af(x) dx + a,(u,’ — uss) + $b,(u,’ — ud_s) 
p =1,2,---,N 
OF = 0 = 2f tp) (apis — ay) + Qty fl tty) (bps — by) 
(3e) 


— 2up(dpiibpi1 — Gpby) — Up (bp41 =~ b, ) i a5 41 + ay 
p = 1,2,---,N —1. 
For the moment let us assume that the values of w, --- , uw, are fixed. In 
this case the set of Equations (3a) and (3b) are seen to be nonhomogeneous and 
linear in terms of the a’s and b’s, and, further, each a, and b, pair may be obtained 


from the solutions of N pairs of linear simultaneous equations. Carrying out this 
computation we obtain 


3 4 
yo art} (uy: + Upp + u5—1)1 (uy » Ups) 
(4a) \ 
— 2(uUp + Upa)d (up, Ups) ’ 
6 
(4b) Bp = rap (2H (tn ta) — (Up + tps) (5 Ups) 
where 


I (ty , Up) = & f(x) dz, 
(4c) eer 


J (up, Upa) = hy af (x) dex. 


Up—1 


Equations (4a) and (4b) represent the general solution to the problem when the 
points uw , --: , Uv; have been selected in advance. When this is not the case then 
to proceed further we substitute the a, and b, from Equations (4) into Equations 
(3c), giving a system of N — 1 nonlinear equations in ™ , --- , uy. to be solved. 
In general no closed form solution for this set of equations can be found. In a later 
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section of this report a numerical method of solution will be described which is prac- 
tical for large-scale computer application. In the case when f(x) is quadratic, some 
simplification of the problem occurs which allows a simple closed form solution to 
be obtained. This result is derived in the next section. 
3. Solution for f(x) Quadratic. In the case that f(x) is quadratic, that is, 
(5) f(z) = rx? + sx +1, 
the function to be minimized is 
N uj 
(6) F=>, / (ra” + sx +t — a; — bx)’ de. 
j=l “uj-1 


The integrand of this function can be simplified by consolidating the linear portion 
of f(z) with the a; and b; . So that, if we let 


a; = (a; — t)/r 


B; = (b; — s)/r, 


(7) 

Equation (6) becomes 
N uj 

(8) F=fr yg (x? — a; — Bx)’ dz, 
j=l Yujii 


and it is only necessary to find the minimum of Equation (8) with respect to the 
a; , 8; and u;. The constant multiplier, r’, in (8) drops out in the derivation of the 
normal equations. 

Following the procedure described in the last section with now f(x) = 2’, the 
a, and 8, are obtained from Equations (4) and have the following simple form: 


ap = —3l(Up + Up)” + 2UpUy—| 


Bp = Up + Upi, p=1,2,---,N 


(9) 


Before proceeding to the solution of the general problem, let us first examine the 
case of N = 2. In this case Equations (3c) reduce to a single equation in u. Sub- 
stituting the results of Equations (9) into (3c) we obtain after some straightforward 
but tedious algebra the following equation in 1% . 


(10) 4( ut. — w)u — 6( U2” — wo )uy” + 4(us — wo )m — (us — w') = 0 
which has for its only real root 

(11) Us = 3(U2 + Ww). 

With w given by (11) and a, B; , a2, Bs by (9), it is easy to verify that 
(12) a + Bit, = a2 + Boy 

and in terms of the original variables 

(13) a, + dy, = ae + dom ; 


that is, the lines of best fit intersect at their common end point. 
To proceed now with the case of fitting N lines to the quadratic, we now make 
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the additional continuity assumption, verified in the case N = 2, that each pair of 
adjacent lines intersect at their end points, uw , uz, --- Uy_; , so that instead of the 
system of nonlinear equations (3c) to solve, we have instead the system 


(14) Ap - BpUp = Apt + Bpiitlp ’ i 1, 2, a a N — i, 
and substituting for a, and 8, from (9) we have 
(15) 2U gly =, Us oy 2UpUp+i <s Up41 eae 1, 2, Ties? N-1 


which further simplifies to 
(16) Up — 2p + Ups = 0, p=1,2,---,N —1, 


Equation (16) may be regarded as a system of linear homogeneous difference 


equations, with two known boundary points, u% and uy . The general solution is of 
the form 


(17) Up = ky > kop, 


and evaluating the constants, k, and k, , from the conditions at p = 0 and p = N 
we have finally 


, 


(18) tp = to + © (uw — we), p =1,2,---,N-—1 


as the values for the common end points. This equation gives the interesting result 
that for the continuous case, the end points are spaced out at equal intervals over 
the whole interval (1% , uw). In applications of the method it has been found most 
convenient first to determine the u’s from (18), next calculate the a’s and 6’s from 
(9) and finally the a’s and b’s from 


a, = Ta,t+t 
bp = rBp + 8, p=1,2,---,N 


where r, s, ¢ are the coefficients of the original quadratic (5). 


(19) 


4. A Numerical Method of Solution. The basic problem is to find a solution of 
the normal equations (3); that is, to find the a’s, b’s and w’s satisfying 


oF a aE : 
(3a) da, wo 0, Pp TT , 1, 2, »N 
q oF os = ¢ “*-. y 
(3b) a, p = 1,2, »N 

oF 

——_ == = 2 eee N _ - 
(3e) ia 0, p = 1,2,---,N —1 


It has been shown that these equations are usually so complicated that the solution 
cannot be obtained directly. A general method in such cases is to assume a trial 
solution and derive linear equations for small additive corrections. The solution of 
this set of linear equations leads to a new set of trial values and the process can be 
repeated until the corrections become negligible. 
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Description of the method becomes much simpler if we adopt a new notation. 
To this end let 


4 =a, 6. = de, ee 5 Oy = Gy, 
Ovi = Dy, On+2 = be, oe Bow = by, 
Beni = Mh, bow42 = Ue, ne O6sn-1 = Un-i. 


The system of equations to be solved is then 
oF 


3" —— - bad ae 
(3’) 20, 0, p =1,2,---,3N—-1 
Let the trial values be denoted by 0,", --- , 03-1 ; expanding (3’) in Taylor’s series 
about the point 6°, and dropping higher-order terms, we have 
oF oF , 0, OF 
20 —_ = 0 ~~ — 6; ~_ 6; a = eee —_ 
(20) 90, a9 + > (0; — 05) Soragp? P= Le BN 1 


where the notation (@F/0@,’), etc. means that the indicated derivative is to be 
evaluated at the point @. Let AQ = @ — @° be the vector with components, 
6, — 6, -*-, Osv-1 — O3v-1, let G be the vector with components, (dF /30,), 
--» , (OF /d0$n-1), and let H be the (3N — 1) by (3N — 1) matrix whose (4, 7) 
component is (d°F/90,00;’). Then in matrix notation the right-hand side of Equa- 
_ tion (20) becomes 


(21) G + AO-H = 0 


and if H is non-singular, then 0", the second approximation to the solution of (20), 
is given by 


(22) Q' = 0° — HG. 


A FORTRAN computer program for the IBM 704 has been written and tested 
which will solve this problem for any function, f(x), which is differentiable. A 
FORTRAN subroutine must be supplied for each f(z), which will compute f(k), 
I(ky , ke), J (ki , ke) and (af(k)/dk). Further details can be supplied by the author. 


5. Application to Lead Susceptibility Curves. As described in the Introduction, 
part of the purpose of the research described in this paper was to find the best 
method of linearizing the lead susceptibility curve. It has been found experimentally 
that for every blend examined, the relationship between octane number of the blend 
and quantity of tetraethyl lead added is adequately given by 


(23) O-N- =k + ke = f(z), 0 


IIA 


zs 3.0 


where O-N- is the octane number of the blend, x is the quantity of lead in c.c., 
and k, , k2 , c are known constants representative of the blend. It has been indicated 
above that this problem may be treated as a single-parameter problem in c alone, 
since k, and ke may be absorbed into the fitted lines. Thus, instead of hav ng the 
integrand of Equation (2) represented by terms of the form 


(ky + kee? — a; — bya)? 
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let us set 


(24) 


(25) 


aj = (a; — ky)/ke 
B; = b;/ke, 
and the integrand becomes proportional to 
(e* — a; — Bz)’. 


The method and computer program described in the last section were used to 
obtain values of the line and end point parameters for fitting two, three and four 
line segments to the exponential function, exp(—czr). The accompanying Tables, 1 
through 3, give the normalized line slopes, 8, the normalized intercepts, a, and the 
values of the end points, u, for the number of lines, N, equal to 2, 3, 4, as a function 


TABLE 1 


Line Parameters, N = 2 


APPROXIMATION OF CURVES BY LINE SEGMENTS 
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} Intercepts Slopes 
c Common end point, #: 
a a2 Bi B 
0.1 | 0.998 0.980 —0.093 | —0.080 1.385 
0.2 | 0.994 0.931 —0.174 —0.129 1.400 
0.3 | 0.988 0.867 —0.246 —0.157 | 1.360 
0.4 | 0.981 0.798 —0.309 —0.170 | 1.316 
0.5 | 0.973 0.728 —0.367 —0.175 | 1.276 
0.6 | 0.966 | 0.661 —0.420 —0.173 | 1.235 
0.7 0.958 0.598 —0.469 —0.168 | 1.196 
0.8 | 0.950 0.540 —0.515 —0.160 1.155 
0.9 0.943 0.488 —0.558 —0.152 | 1.116 
1.0 0.935 0.442 —0.600 —0.142 | 1.080 
1.2 | 0.923 0.363 —0.680 —0.124 1.008 
1.5 | 0.906 0.275 —0.793 —0.101 | 0.912 
TABLE 2 
Line Parameters, N = 3 
Intercepts Slopes End Points 
a a as Ar | Bs Bs “1 “2 
0.1 | 0.999 | 0.993 | 0.976 | —0.096 | —0.088 | —0.079 | 0.785 | 1.783 
0.2 | 0.998 | 0.975 | 0.919 | —0.186 | —0.156 | —0.125 | 0.782 | 1.775 
0.3 | 0.995 | 0.943 | 0.840 | —0.266 | —0.202 | —0.146 | 0.835 | 1.820 
0.4 | 0.992 | 0.907 | 0.756 | —0.340 | —0.237 | —0.154 | 0.830 | 1.810 
0.5 | 0.988 | 0.870 | 0.673 | —0.411 | —0.265 | —0.153 | 0.801 | 1.773 
0.6 | 0.985 | 0.835 | 0.595 | —0.478 | —0.285 | —0.147 | 0.774 | 1.738 
0.7 | 0.981 | 0.798 | 0.523 | —0.541 | —0.299 | —0.137 | 0.756 | 1.708 
0.8 | 0.978 | 0.767 | 0.460 | —0.604 | —0.312 | —0.127 | 0.724 | 1.663 
0.9 | 0.974 | 0.737 | 0.404 | —0.665 | —0.323 | —0.117 | 0.694 | 1.618 
1.0 | 0.971 | 0.711 | 0.355 | —0.725 | —0.333 | —0.107 | 0.666 | 1.573 
1.2 | 0.966 | 0.664 | 0.275 | —0.841 | —0.349 | —0.088 | 0.615 | 1.487 
1.5 | 0.959 | 0.608 | 0.192 | —1.012 | —0.371 | —0.066 | 0.548 | 1.365 
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TABLE 3 
Line Parameters, N = 4 

































































Intercepts Slopes End Points 

a a aw | a Bi B2 Ba Ba m | we “3 
0.1 |1.000/0.996)0.988/0.973/—0.097| —0.091|—0.085 —0.078/0.571/1 .268)2 .077 
0.2 |0.999/0.985)/0.955/0. 908) —0.188)—0.166|—0.143)—0.120/0.607/1 .296)/2.096 
0.3 |0.997/0.968/0.910/0.824|—0.274|—0.226|—0.181|/—0.140)0.607/1 .293)2.087 
0.4 |0.996)/0.948)/0.859|0.734| —0.356| —0.276|—0.206| —0.146/0.595)1 .275/2.068 
0.5 |0.994'0.927\0.805\0.646| —0.434/—0.318 —0.221)—0.143)0.577/1 .246)2 .040 
0.6 |0.992/0.904/0.752)0.562)—0.508)—0.352| —0.228| —0.134/0.563)1 .222/2.016 
0.7 |0.989/0.880)0.700/0.486)—0.580)—0.381 —0.231|—0.124/0.549/1 .199|1 .991 
0.8 |0.987/0.859)0.654/0.420| —0. 650) —0.409| —0.232|—0.112/0.530)1 .165)1 .953 
0.9 |0.985/0.841/0.614/0.363| —0.720) —0.435)-—0.233|—0.101/0.507/|1.124/1.910 
1.0 |0.984/0.824/0.577/0.313| —0.789| —0.460|—0.232| —0.091/0.485/1 .084/1 . 866 
1.2 |0.981/0.795|/0.516)0.234/ —0.925| 0.508 —0.230| —0.072|0.445)1 .008)1.778 
1.8 0.977/0.760(0.445)0.154)—1.126/—0.574 —(.227|—0.051 0.393/0.907)1 .649 

TABLE 4 
Maximum Fitting Errors 
max(e~°* — aj — Bjx) 
N=2 | N=3 N=4 

0.1 0.0016 | 0.0010 0.00055 

0.2 0.0059 0.0032 0.0016 

0.3 0.0119 | 0.0050 0.0028 

0.4 0.0189 0.0080 0.0043 

0.5 0.0265 | 0.0114 | 0.0062 

0.6 0.0344 0.0150 0.0083 

0.7 0.0423 | 0.0190 0.0106 

0.8 0.0500 0.0223 0.0127 

0.9 0.0574 0.0255 0.0145 

1.0 0.0645 0.0286 0.0162 

1.2 0.0774 0.0341 | 0.0193 

1.5 0.0936 | 0.0410 0.0231 





of the exponential parameter, c. The actual intercepts and slopes, a and b, for any 
particular problem are obtained from 


@; = ky + koa; 
bj = keB;. 


It has been found preferable in practice, for reasons of accuracy, to determine 
the common end points of the lines from the relation 


(26) 


a; + bu; = Ojai + Oj410; 
giving 
(27) PR as NE 2 t 
* Dyn — By 
In order to provide the user of this method with a guide to the number of line 
segments that should be used for any particular application, an error analysis was 
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carried out. The maximum difference between the exponential function, e, and 


the lines a; — b;z were computed for the two, three, and four line cases. The results 
are shown in Table 4 where 


max(e * — a; — 6;2), 0s75 30 


is given as a function of c for these cases. To transform the error given in the nor- 


malized form to error in terms of octane number for any particular case, the trans- 
formation 


Error (O-N-) = ke-Error (Table) 
is used where k, is the constant shown in Equation (23); that is, in the equation 
O-N- =k + ke™. 


The tables given for the lead susceptibility curves can be used for curves of a 
more general character, namely those of the form 


(28) giz) =nt+ret+re™, Wo 


In this case the least squares function is 


lA 
Nn 
IA 
s 
= 


N vj 
(29) F=>, / (ry + rez + ree" — a; — b,z)* dz. 
j=l Swj_1 
Let a new variable x be defined by 
xr=3- s-* q dx = __ 3dz é 
Wy — Wo Wyn — Wo 
The range of x isO S x S 3.0, and F becomes 
N uj 
(30) F« > | (e* — a; — Ba) dx 
j=l uj—1 


which is identical in form, except for a proportionality constant, to the previous 
case. The transformation equations are 


a; = J la; — ri + (b; — re) wo) 
r3 


(31) B; = (>) (b; — rs) 


q { Wj — Wo 
u; = 3([ ———__}. 
Wn — Wo 


The inverse transformation equations to the original units are 
38; 
Wn — Wo 


(32) a; = 71 + rsa; — (b; — re) wo 


Wo + (5) U;. 


1. W. W. Garvin, H. W. Cranpatu, J. B. Jonn, R. A. Speviman, ‘Applications of 
linear programming in the oil industry,’’ Management Sci. 3, 1957, p. 407. 


bj =1r2+ 
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TECHNICAL NOTES AND SHORT PAPERS 


Numerical Integration Using Sums of 
Exponential Functions 


By F. C. Ledsham 


Let y be a function of the independent variable z, and let 


dy 

(1) . es a 

Suppose that the numerical values of g be known for the n + 1 discrete values 
of x given by x=2» + rh, where r runs through the integers from zero to n and h is 
some fixed increment. Let f(x) be some function of z chosen to coincide with q at 
these particular z-values. For convenience we shall denote g(x» + rh) by q,, with 
a similar convention for other functions. In order to obtain a numerical estimate 
of the increment of y over any chosen range of xz, say y, — ym, one may then use 
the approximation 


zo+kh zo+kh 
(2) m-m= { qa] far, 
zo+mh zo+mh 

the accuracy of which will depend upon the choice of the function f(x). Note that 
k and m may have any real values (positive, negative or zero) and, in particular, 
are not confined to the integers nor to the region from zero to n. 

It is customary to take f(x) to be the lowest order polynomial satisfying the 
conditions specified above. As is well known, it is not then necessary to determine 
this polynomial, and it is possible to write 


zotkh n nm 
(3) [fae =hY ACh, myn;r) f, = hD Alk,m,n;7) a, 
z r=(0 r=0 


ot+mh 


where the A’s are constants satisfying the equation 
(4) D A(k, m,n;r) = k — m, 
r=0 


and which may be tabulated once and for all as functions of k, m, n and r. 

If m = 0 and k = n, then this last procedure leads to the Newton-Cotes series 
of formulae, of which Simpson’s rule (for n = 2) is probably the best known. If we 
put k = m + 1, and also let this quantity equal n + 1 or n, then we get, respec- 
tively, extrapolation and check formulae which may be used for step-by-step 
numerical integration of first order differential equations. Other combinations of 
m and k also lead to useful formulae, and a selection of these (for various values of 
n) is included in a paper by Bickley [1]. 

It has been pointed out by Greenwood [2] and by Brock and Murray [3], that a 


Received April 5, 1960. 
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polynomial form of f(z) does not always provide the best approximation to q(z). 
In particular, these authors give examples in which f(z) is better taken as a sum 
of exponential functions, in the form 


(5) f(z) = ya exp(aiz), 


where the a; are chosen to suit the problem concerned, and may be complex or 
zero. 

Greenwood considers two particular forms of equation (5), both involving real 
values only of the a; . For his first type he puts a; = i; while, for the second, he 
puts a; = i — 4n where i, in both cases, runs through the integers from zero to n. 
The second (symmetric) case is only applied by Greenwood to even values of n. 
He also takes m = 0 and k = n, to produce formulae analogous to those of the 
Newton-Cotes series. 

Brock and Murray are concerned with the step-by-step numerical integration 
of first order differential equations. They consider more general cases in which the 
a; are complex, and tailored to fit the particular problem to hand—with the help of 
earlier (and less accurate) solutions of the differential equations concerned. 

The papers of Greenwood and Brock and Murray both contain practical ex- 
amples of uses of the ideas expressed above, and discuss-the magnitudes of the 
errors involved. 

Using equation (5), with particular values given to the coefficients a; , we may 
obtain an equation of the form 


ezotkh n n 
(6) fdx = > Blk, m,n;h,r)f, = > B(k, m,n; h,r)q. 
Zo+mh r=0 r= 


Unlike the corresponding equation (3), h does not occur naturally as a factor on 
the right hand side of this equation, and the coefficients B have to be recalculated 
for every change in this quantity. Incidentally, Brock and Murray, in their work, 
do take out a factor h—and consequently calculate coefficients equivalent to B/h 
in our notation. 

If one of the a; be zero, then we have the relationship 


(7) >: B(k, m, n;h, r) = h(k — m), 
r=( 


corresponding to equation (4). If none of the a; be zero, then this last equation 
does not hold exactly—though it remains approximately true, and serves as a 
useful check during the calculation of the B’s. 

If one of the a; be zero then let this be a, . We then have 


zo+kh 
i f dx = a,h(k — m) 
(8) zo+mh p 
+> al lexp {a;(2% + kh)} — exp {a(x + mh)}). 

i=? Qj 


(is) 


If none of the a; be zero then the first term on the right hand side of this equation 
drops out, and the restriction i ¥ s is omitted from the summation. 
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From the n + 1 equations 
(9) f=a= 2 ai exp{a;(xo + rh)}, 


the a; , or rather a; exp(ai2o), may be found in terms of the g, , enabling us to put 
equation (8) into the form of equation (6). In previous applications of these ideas 
[2, 3] this step, or its equivalent, has been carried out numerically. It is therefore 
considered that there may be some interest in analytical formulae into which one 
might substitute directly in order to obtain the required coefficients B(k, m,n; h, r), 
such as those given below. 
Write 
(10) & 


a; exp(aizo) | 
exp (ah) _,f 


so that equations (9) may be written in the form 


(11) fe =a = D biti, (r = 0,1, ---- n). 
i=0 


For any value of n, it may be verified that the solution of these n + 1 equations 
leads to the following symbolic equations for the 6; : 


(12) —_— Tix: (t, — q) 
LLixi (t; — ts) 
where, after expansion, the powers of the q’s are lowered to represent suffixes, and 
the term of the numerator originally independent of g is taken as the coefficient of 
qo . For example, if n = 2 we would have 
‘a (t, — g)(e — q) 
(i — &)(te — b) 
which would be interpreted as giving 
a titego — (ti + teqr + m 
eas ? 
(4 — to)(te — tb) 





(symbolically ), 


0 





and similarly for b; and b . 
If, as in equation (8), a, be zero then a, = b, and that equation becomes 


zotkh n b; 
(13) [fac = bak —m) + & Bas 0%. 
zo+mh nA ay 


Substituting into equation (13) the values of b; calculated from equations (12) 
leads immediately to the required form (6). It will be noted that the b; do not 
involve k and m, so that they do not have to be recalculated if integrations are 
required over ranges corresponding to more than one pair of values of these param- 
eters. 

If two of the a; be complex conjugates, then so also will be the corresponding 
terms of the summation on the right hand side of equation (13). Instead of cal- 
culating both those terms completely, it is therefore only necessary to find the real 
part of one of them and then to double it to obtain the sum of the two. 

Equations (12) and (13) have been used to recalculate some of the coefficients 
quoted in the papers mentioned above. The first example taken was Greenwood’s 
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symmetrical case with n = 6. The second was taken from the paper by Brock and 
Murray, with n = 3 and the a; consisting of two pairs of complex conjugates. 
Using the full 10-figure capacity of a desk machine, cancellation reduced the 
accuracy of the results obtained to some four significant figures in the first example 
and six in the second. The procedure given here should, however, be readily adapta- 
ble to electronic digital computers, and the increased capacity of those machines 
should enable the coefficients B(k, m, n; h, r) to be calculated to any accuracy 
likely to be required in practice. 

Sevenoaks 

Kent, England 


pt — G. Bicker, ‘‘Formulae for numerical integration,’’ Math. Gaz., v. 23, 1939, p. 


2. R.E. GREENWOOD, “‘Numerical integration for linear sums of exponential functions,” 
Ann. Math. Statist. v. 20, 1949, p. 11 


3. P. Brocx & F. J. Murray, ‘The use of exponential sums in step by step integration,”’ 
MTAC, v. 6, 1952, p. 63-78 and 138-150. 


New Factors of Mersenne Numbers 
By Edgar Karst 


I have tested for prime factors all the Mersenne numbers 2?—1 corresponding 
to prime exponents p in the interval 3000 < p < 3500. The limit of the search for 
factors was 9p? when no factor was previously known; otherwise the limit was 3p’. 

The nineteen new prime factors of Mersenne numbers found by this search are 
displayed in the following table. Factors corresponding to smaller values of p have 
been listed in a paper by Brillhart and Johnson [1]. 


p New factors of 2? — 1 
3037 145 777 
3041 5 565 031 
3067 22 063 999 
3083 15 914 447 
3119 230 807-14 222 641 
3121 31 509 617 
3167 12 237 289 
3181 127 241 
3191 40 895 857 
3253 46 452 841 
3257 4 032 167 
3299 19 873 177 
3329 665 801-1 005 359-26 225 863 
3391 1 519 169 
3433 5 952 823-12 688 369 


An extensive table by Riesel [2] includes smaller prime factors of Mersenne 
numbers corresponding to p = 3037, 3041, 3119, 3121, 3181, 3257, 3299, 3329, 
3391, and 3433 in the preceding table. 

Brigham Young University 
Provo, Utah 

1. Jonn Brritwart & G. D. Jonnson, “On the factors of certain Mersenne numbers,” 
Math. Comp., v. 14, 1960, p. 365-369. 

2. H. Rresex, “Mersenne numbers,” MTAC, v. 12, 1958, p. 207-213. 


Received September 14, 1959; in revised form July 26, 1960. 











SURVEYS AND EXPOSITORY PAPERS 


A Bibliography on Approximate Integration 


By A. H. Stroud 


This bibliography is an attempt to provide a reference to the knowledge con- 
cerning approximate integration methods, and to make better known the recent 
contributions to this field. The two main subjects of this bibliography are quadra- 
ture (or numerical integration) formulas—weighted linear sums of values of the 
integrand (and possibly derivatives of the integrand)—and finite difference ap- 
proximations. Papers dealing with error theory or other properties of these methods 
are included. 

Some topics which are not included are: applications of integration methods 
to the solution of differential equations; graphical methods; and tables or rational 
function approximations for functions expressed as integrals. 

There are several reference works which give good treatments of parts of this 
subject. Some of the best of these—Hi_pEeBranp [1], BuckrneHam [1], and 
MineEur [3}—give good summaries of the well-known integration methods for 
functions of one variable. 

The recent book by V. I. Krytov [7] is probably the best all around book to 
date concerned with approximate integration for functions of one variable. 

The small book by Nikou’skm [3] treats “best”? quadrature formulas and ex- 
tremal problems related to quadrature formulas. 

The two papers by Hammer [2] and Stroup [3] together give a good survey of 
what is known concerning quadrature formulas for functions of more than one 
variable. 

KantorovicH & Kry.ov [1] give a good discussion concerning methods for the 
approximate solution of integral equations. 

For references to works giving quadrature formulas, interpolation tables or 
tables of integrals, see FLercHer, Mituter & RosENHEAD [1]. 

A few references of historical importance have been included here. A more 
complete list of papers prior to 1864 is given by BrerEens DE Haan [1], who gives 
60 papers concerned with approximate evaluation of integrals. 

The titles of the foreign language papers, other than French, Spanish and 
German, usually have been translated into English. Included are the volume, 
page references to reviews in Mathematical Reviews (MR) and Zentralblatt fiir 
Mathematik (Zb!1.). 

I am indebted to Prof. P. C. Hammer who suggested and encouraged this work, 
and to Prof. J. C. P. Miller for his corrections and suggestions. 

The journal abbreviations are those given in Mathematical Reviews, Index for 
Volume 19, pages 1417-1430. 


Received December 15, 1959; in revised form June 13, 1960. The preparation of this bibliog- 
raphy was supported in part by the Office of Ordnance Research, and in part by the Graduate 
Research Committee of the University of Wisconsin. 
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REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


1[A-F, G-I, K, M}—Nia M. Burunova, Spravochnik po matematicheskim tab- 
litsam, Dopolnenie N. 1 (Handbook on Mathematical Tables, Supplement No. 1), 
Izdatel’stvo Akad. Nauk SSSR, Moscow, 1959, 184 p., 26 cm. Price 9 rubles. 


This is a first supplement to the important work by Lebedev and Fedorova [I], 
which appeared in 1956. The arrangement is similar to that of the original, with a 
Part I in fifteen chapters listing the contents of the various tables, which are re- 
ferred to by number, and a Part II giving the numbered references, separately for 
each chapter. 

With two small exceptions, the titles of the chapters are unchanged. The trans- 
lated heading of Chapter IV, formerly “Decimal and natural logarithms,” is now 
simply “Logarithms’’; this allows the inclusion of a few tables of logarithms to base 
2. Chapter XV, besides the former “Prime numbers, factors, products, quotients 
and fractions,’ now covers also “Conversion from one system of numeration to 
another’’; this accommodates a few tables which deal with binary-decimal and simi- 
lar conversions. A three-page list of functions, giving references to both Handbook 
and Supplement, also shows a few additional categories not involving alteration of 
chapter headings. 

As in the original Part II, references to works in Russian are naturally given in 
Cyrillic characters, but the far more numerous references to other works are given in 
roman characters. A mathematician knowing no Russian could probably contrive 
to use the greater part of the book. 

Like the original work, the supplement is very welcome. 

A. F. 


1. A. V. Lesepev & R. M. Feporova, Spravochnik po matematicheskim tablitsam, Moscow, 
1956. See RMT 49, MTAC, v. 11, 1957, p. 104-106. 


2{|A-E, G, J, M, Q]—Frrepricu O. Rineies, Mathematische Formelsammlung, (7th 
ed.) Sammlung Géschen Band 51/5la, Walter de Gruyter & Co., Berlin, 1960, 
320 p., 16 em. Price DM 3.60. 


This little German volume, from the Sammlung Géschen, is an enlarged version 
(42 extra pages, 3 extra figures) of the 1956 edition written by O. Th. Biirkeln. It is 
a classified collection of formulae and standard theorems from all branches of 
(essentially) undergraduate mathematics. As such, there is some, but not too much, 
overlap with similar collections, such as the well-known Smithsonian Mathematical 
Formulae. It is entirely different, however, from the numerical analysis handbook 
in the Sammlung Géschen, Formelsammlung zur praktischen Mathematik, by Giinther 
Schulz. 

One (probably not anticipated) use to students of this little boox is as a pleasant 
means of picking up a basic mathematical German vocabulary. The sentences are 
short, factual, and not unduly encumbered by grammar or philosophy. 

The sixteen chapter headings are: Arithmetik und Kombinatorik, Algebra, Zah- 
lentheorie, Elementare Reihen, Ebene Geometrie, Stereometrie, Ebene Trigono- 
metrie, Spharische Trigonometrie, Mathematische Geographie und Astronomie, 
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Analytische Geometrie der Ebene, Analytische Geometrie des Raumes und Vektor- 
rechnung, Differentialrechnung, Integralrechnung, Funktionentheorie und kon- 
forme Abbildung, Differentialgeometrie, and Differentialgleichungen. 

D.S. 


3[C, L}—L. K. Frevet & J. W. Turtey, “Seven-Place Table of Iterated Log. 
(1 + x),” The Dow Chemical Company, Midland, Michigan, 1960. Deposited in 
UMT File. 


The n-fold iteration of log, (1 + 2x) is explicitly given by 


nm < 1f1,n-19 ‘ 
oe (=) n[ b+ ri Con" |, 

2 i=3 t 2 «a 
where the C..’s are numerical rational fractions. Using nine terms of this expansion 
the authors have compiled a 7D table containing 4000 different entries for n = 
0(.05)10 and z = 0(.05)1. All computations were programmed on the DATA- 
TRON 220, and the output in tabular format was printed directly by a Teletype 
printer. The recorded values are accurate to +4-10~’. 





Autuors’ SUMMARY 


4[F|—C. L. Baxer & F. J. Gruenspercer. The First Six Million Prime Numbers, 
The RAND Corporation, Santa Monica, published by The Microcard Foundation, 
Madison, Wisconsin, 1959, 8 p., 16 x 23 em. + 62 cards, 7.5 x 12.6 em. Price $35.00. 


This unusual table is arranged on microcards which present 124 photographs. 
Each photograph (except the first and last) displays 39 pages of tabulation. Each 
page lists 1250 prime numbers. The primes range from 1, which is counted as a 
prime, to 104395289. Each line of a page contains 25 consecutive primes. The first 
prime in the line is given completely; only the last three digits of the other 24 
primes are given. The rank of a prime, once it is located in the table, is given by 
an cbvious formula in terms of its page number, line number, and position in that 
line. It is just as easy to find isolated values of r(x), the number of primes Sz. 

It is obvious that the very high condensation of information achieved in this 
list of primes is not won without some difficulty, namely, the fact that it is invisible 
to the naked eye. A quite strong pocket magnifying glass or a microcard reading 
machine is required to read the data. Any use of the table other than finding 
whether or not a given number is prime or evaluating (x) for isolated values of 
zx is really impractical. For example, to determine the number of twin primes in the 
56th million or to calculate a sum involving consecutive primes—problems in 
which one must carefully keep one’s place in the table—would be difficult indeed. 
Of course, such problems should be done by an electronic digital computer any- 
way. 

There is an interesting description of the IBM 704 program used to generate 
the list of primes on punched cards, which also contain the differences between 
consecutive primes. 

D. H. LeaMer 


University of California 
Berkeley, California 


Eprror1at Nots.—If 2 is counted as the first prime, which is the current practice, then 
the six millionth prime is 104395301, which happens to be the first member of a prime pair. 
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5(F].—H. Gupra, M. 8. Cueema, A. Menta & O. P. Gupta. Edited by J. C. P. 
Mier, Representation of Primes by Quadratic Forms, Royal Society Mathe- 
matical Tables No. 5, Cambridge University Press, New York, 1960, xxiv + 
135 p., 29 em. Price $8.50. 


The eight tables in this work give solutions a, b and k, n to two related Dio- 
phantine equations 


sp=a+D and ckp =n’ +D, 


with «x = 1 and 2, D = 5, 6, 10, and 13, and for all primes p < 100,000 for which 
solutions exist. In gach table there are approximately 2400 primes, that is, about 
one-fourth of all primes < 100,000. 

The subtitle indicates that this volume is Part I of a larger work and implies 
that other values of D will be forthcoming. The significance of these particular 
values, 5, 6, 10, and 13, is that in all quadratic number fields, R(«/—D), these 
are the smallest positive D’s for which unique factorization of the algebraic integers 
is lacking. In all four cases the class number is 2 (there are two classes of ideals) 
and this is associated with the two values of «. For about one-half of the p’s for 
which —D is a quadratic residue a solution exists for x = 1, and for the remaining 
one-half, for x = 2. It may be noted that the x = 1 table is always somewhat shorter 
than its companion x = 2 table. This is as expected, since there are generally more 
primes of the form 4Dm + N than of the form 4Dm + R if R is a quadratic residue 
of 4D and N is not (see MTAC, v. 13, 1959, p. 272-284). 

There is an interesting, twelve-page introduction to the background ideal and 
class number theory. It is boldly stated there, p. xii, that unique factorization 
exists for square-free D > 0 only if D = 1, 2, 3, 7, 11, 19, 43, 67, and 163. How- 
ever, this has never been fully proven. This background theory culminates in 
several theorems due to H. P. F. Swinnerton-Dyer. 

The tables were done by hand, “with the help of two long strips of paper—A 
and B,” and were carefully checked in a variety of ways. Why the solutions to 
kp = n°’ + 6 for p = 7 and p = 31 are listed, on p. 36, as k, n = 6, 6 and 22, 26 
respectively rather than the obvious 1, 1 and 1, 5 is not clear to the reviewer. 
But this slip seems to be exceptional. 

The tables may be used as lists of prime ideals in the four fields, R(4/—D). 
The dust jacket suggests a second application, that they “may contribute to the 
understanding of such unsolved questions as ‘Is the number of primes of form 
n’ + 5, orn’ + 6, etc., infinite or finite?’ ” In this respect, however, it may be re- 
marked that the data here are rather meager since the primes are <10°, while 
in 10 minutes an IBM 704 can count such primes up to 3.24-10" (see MTAC, v. 
13, 1959, p. 78-86). In fact the 16th prime of the form n’ + 5 is 86441 and is 
the largest listed in this volume, while the 4368th prime of that form is 32,371, 
926, 089. Nonetheless, the data are sufficient to indicate at least rough agreement 
with the Hardy-Littlewood conjecture. Since 


Ls(1) = 1.405 > L6(1) = 1.283 > Lyo(1) = 0.993 > Ly3(1) = 0.871 


we should expect the relative density of primes tc increase as we progress from 
n +5ton’ + 6 ton’ + 10 ton’ + 13 (Math. Comp., v. 14, 1960, p. 324-326). 
This is indeed the case. 
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An apparent anomaly concerned n’® + 13, where there were numerous primes 
from n = 0 to 68 and from n = 264 to 298, but none in between. This striking 
maldistribution was most alarming, and threatened dire consequences to the Hardy- 
Littlewood conjecture, until the real explanation was found—pages 105 to 120 
were missing in the reviewer’s copy. Aside from this gross lapse, the volume has 
the usual elegance of the Royal Society Mathematical Tables. 

D. S. 


6[F, L|—C. B. Hase.erove in collaboration with J. C. P. Mituer, Tables of the 
Riemann Zeta Function, Royal Society Mathematical Tables No. 6, Cambridge 
University Press, New York, 1960, xxiii + 80 p., 29 cm. Price $9.50. 


These important and fascinating tables are concerned primarily with ¢(} + it), 
the zeta function for real part 4, and with its zeros. This complex function is ex- 
pressed both in cartesian and polar forms: 


¢ (3 + it) = we (4 + it) + ise (4 + it) 
ine foe. 


In the latter, @(¢) is continuous, with @(0) = 0, and the signed modulus Z(t), 
given by 


_ we T+ 48). ,; 
Z(t)=-r TG + 31) | o(3 + it), 
changes sign at every zero. The nth zero, y, , is the nth solution of Z(y) = 0 and 
the nth Gram point, g, , is the solution of @(g,) = n. 

Table I gives Rf(4 + it), 9€(4 + it), Z(t), and @(t) to 6D for ¢ = 0(0.1)100. 
Ge(1 + it) and gf(1 + it) are also listed. 

Table II gives Z(t) to 6D for ¢ = 100(0.1)1000. 

Table III has two parts. Part 1, for n = 1(1)650, gives yn, gnt, and @, = 
(1/x)pht’(} + ty.) to 6D and | ¢’(4 + 7y,)| to 5D. Part 2, for n = 651(1)1600, 
gives y, to 6D and | ¢’(4 + 77,)| to 5D. 

Table IV gives Z(t) to 6D for four other ranges of ¢, 


¢ = 7000(0.1)7025, 17120(0.1)17145, 100000(0.1)100025, 250000(0.1)250025. 


Also for these four ranges are given the zeros to 6D and derivatives to 5D. There 
are 28, 32, 38, and 42 zeros in the four ranges respectively. 

Table V gives (1/r)phI'(} + it) to 6D for ¢ = 0(0.1)50(1)600(2) 1000. 

The inclusion of g, and ¢, in part 1 of Table III allows the reader to study 
Gram’s “Law” which states that the zeros and Gram points are interlaced: 


Yn-1 < Jn—2 < Yn < Jn-1 < Yn+1 - 


A violation occurs if | ¢, | > 4. The first violation is for n = 127. The first double 
violation is for n = 379 and 380—.e., there are three Gram points between 737 and 
zs0 - In all, there are 22 violations in these 650 zeros. Gram’s ‘‘Law’’ may also be 
expressed by saying that the complex ¢(3 + it) approaches its zero via the 4th 
or 3rd quadrant. Thus the following statistics for these 650 zeros are of interest: 
4th quadrant, 320 cases; 3rd quadrant, 308; 2nd quadrant, 13, and Ist quadrant, 9. 
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oU% +iff 
Fic. 1.—t (4 + it) in the complex plane, for 0 < ¢ S 30, 


This suggests the conjecture* that Limy.. (1/N) >? ¢, = 0, that is, that 
ia + tn) 


is positive real in the mean. Since violations against the “Law” are associated with 
exceptionally close zeros or exceptionally large values of Z(t), some extremes dis- 
covered by D. H. Lehmer are of interest, and are listed in Table IV: 


vm = 7005.062 866 and ymsi = 7005. 100 565 (m = ?) 
Z(17123.1) = 18.955 257. 


The tables of ¢(4 + it) and of Z(t) were both computed on two different ma- 
chines, EDSAC and the Manchester Mark I, from two different asymptotic formu- 
lae, one due to Gram, the second a modification of the Riemann-Siegel formula 
due to Lehmer. A great deal of difficulty was experienced in eliminating discrepan- 
cies between the two methods, and three or four erroneous tables had to be dis- 
carded. 


* Note added November 10, 1960. For the first 650 zeros this mean value, ¢., shows such a 
marked trend toward zero that it is even probable that |Z @,|/ ~/N also tends toward zero. 
Taking only the ‘‘worst” points—that is, where |2@,| reaches a new maximum—sample 
values of |z on| / WN are .0513, .0437, 0380, .0328, .0279, 0251, and .0217 for n = 8, 33, 64, 126, 
256, 379, and 606, respectively. It should be noted that ¢, oscillates as it tends towards zero 
—up ton = 650 it changes sign 265 times. In particular, all 22 of the Gram violations, | be | > 4, 
are associated with such sign changes. But this latter is almost inevitable up to n = 650, 
since the extreme values of = @, up to this limit are only +.520578 for n = 567 and 
— 533885 for n = 606. 
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As to uses of the tables, Haselgrove himself has used the values of the first 
600 zeros in his disproof of Pélya’s conjecture. For later results on this theme, see 
R. 5S. Lehman, “On Liouville’s function,’ Math. Comp., v. 14, 1960, p. 311-320. 

It is not inconceivable that study of these tables of ¢(4 + it) may inspire 
some investigator to a new approach to the Riemann Hypothesis. Similarly, the 
table of ¢(1 + 7) can be studied in connection with proofs of the prime number 
theorem. For both of these “uses,” however, a graphical presentation is highly 
desirable, and it is regretted that a good collection of graphs was not included in 
this volume. For example, a graph of ¢(3 + it) in the complex plane versus the 
parameter {—say from 0 to 30—is particularly interesting. See Fig. 1. Problem for 
the reader: If, in Fig. 1, the variable ¢ is thought of as time, explain the initial 
counterclockwise motion in the orbit and the subsequent clockwise motion with a 
shorter and shorter mean period. Hint: Consider the formula 


ro) = (1-2) 


2 n=l n*® 


for a fixed value of s = 4 + it. For what mean value of n does a block of consecu- 
tive terms in the series have all its terms in phase; in which block are alternate 
terms out of phase; and finally, which terms add in an essentially random manner? 

Since it is known that Haselgrove has also computed complex zeros of the 
closely related L(s) and other Dirichlet character series, and since he has not in- 
cluded these results, it would be desirable to issue a companion volume to make 
these related tables generally available. 

Two very minor errors were noted in the Introduction. 

Page xii, line 5: change {(p,) to ¢’(pn). 

Page xvi, line 6 from bottom: change (3.7) to (3.32). 
With respect to Liénard’s table of ¢(n) for n = 2(1)167, which is mentioned on 
page xx, a recent extension should be noted: J. W. Wrench, Jr., ‘Further evalua- 
tion of Khintchine’s constant,’’ Math. Comp., v. 14, 1960, p. 370-371. 





D. S. 


71G, K, X, Z]|—AntHony Ratston & HersBert S. Witr, Mathematical Methods 
for Digital Computers, John Wiley & Sons, Inc., New York, 1960, xi + 293 p., 
27 em. Price $9.00. 


This book contains contributions from twenty-four research workers in numerical 
analysis and related fields. Two of the contributors serve as editors, and, as the 
title implies, all contributors are definitely high-speed-computer oriented. It is 
interesting to note that ten of the twenty-four authors are from universities. 

In the introduction it is stated that ‘the major purpose of this book is to present 
many—but by no means all—of the more commonly used tools of the modern 
numerical analyst along with some of the more promising newly developed 
methods. The motivation behind this presentation is not only to gather together in 
one place a partial survey of modern numerical methods but also, in each case, to 
acquaint the reader with the interplay between computer capabilities and processes 
of analysis’’. 

The editors have done a good job in carrying out the purpose of the book. They 
have coordinated the efforts of their colleagues in a remarkable way. About twenty 
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important and fairly representative problems are discussed by experts in various 
fields and the presentations, for the most part, are organized in a uniform manner. 
All but one of the authors describe programs for solving their problems in terms 
of the following topics: 
. The function of the program 
. A mathematical discussion of the problem 
. A summary of the calculation procedure 
. A complete flow diagram 
. A box by box description of the flow diagram 
. Standard subroutines required by the program 
. A representative sample problem 
. Memory requirements 
. An estimation of the running time 
10. A list of references 
Actual coding of the programs is omitted, as is any reference to a particular com- 
puting machine. 

The book is divided into six parts: generation of elementary functions, matrices 
and linear equations, ordinary differential equations, partial differential equations, 
statistics, and miscellaneous methods. 

Part I contains an interesting discussion of the generation of elementary func- 
tions by means of polynomial and rational approximations. Ten of the more com- 
mon functions are examined in detail, and comparisons of the errors obtained using 
various approximation methods are included. Because of the nature of its subject 
matter this part has a format different from that discussed in Paragraph 3. Part II 
contains six chapters covering matrix inversion, systems of linear equations, and the 
matrix eigenvalue problem for symmetric matrices. Part III contains four chapters 
concerned with the numerical solution of ordinary differential equations, and Part 
IV contains five chapters discussing parabolic, elliptic and hyperbolic partial 
differential equations. There are four chapters in Part V covering multiple regression 
analysis, factor analysis, autocorrelation and spectral analysis, and the analysis of 
variance. The six chapters in Part VI discuss methods for numerical quadrature, 
Fourier analysis, linear programming, network analysis, and the solution of poly- 
nomial equations. 

This book should serve as an excellent reference work for those workers who need 
to solve problems using electronic digital computers. For this reason the format, 
described above, is especially good. The book is not intended to be a textbook for 
introducing the subject of numerical analysis, but most students of the subject will 
find useful information here. Lists at the end of the various chapters contain over 
two hundred references relating to the problems under consideration. 

The reviewer would like to see a supplementary volume appear containing many 
of the important methods not covered in the present volume. For instance, Part 
II does not cover the Givens method for finding the eigenvalues of symmetric 
matrices. Omitted entirely is any discussion of methods for finding the eigenvalues 
of arbitrary matrices. Other readers will find additional important methods they 
would like to see covered in a second volume. 


— 


Connor WW 


Rosert T. GREGORY 
The University of Texas 
Austin, Texas 
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8[G, Z].—Franz E. Houn, Applied Boolean Algebra, An Elementary Introduction, 
The Macmillan Company, New York, 1960, xx + 139 p. $2.50 (paperback, 
offset print). 


The author gives as much information as it seems possible to include conven- 
iently ina text of this length starting with the veriest elements concerning the use of 
Boolean algebra in the design of switching circuits—especially for digital computers. 

Examples of engineering applications are given, but no complete design of a 
very extensive arithmetic unit is undertaken. Hence, the ultimate dependence of 
the designers of computers on Boolean algebra is not completely illustrated, although 
it is strongly and correctly implied. 

The crux of any applications of Boolean algebra to the design of switching cir- 
cuits lies in computational schemes for writing fairly efficient statements of Boolean 
propositions. This problem is faced by the author, but only to an extent which per- 
mits the reader (and problem worker) to understand the nature of the difficulties 
which are encountered and some of the procedures which promise to be helpful. In 
this regard the pamphlet is no less informative than most of the other textbook 
material available, but additional reference to computational efforts would have 
been welcome. 

The electrical elements to be used are described abstractly in a reasonable way, 
and altogether the presentation is self-contained, lucid, and reasonably illustrated by 
problems. No sophistication is required in the reader except for motivation. 

Attention is not restricted to circuit design, and the standard applications of 
Boolean algebra are treated to an extent which is indicated in the chapter headings 
listed below: 

Introduction 

Boolean Algebra as a Model of Combinational Relay Circuitry 
Boolean Algebra as a Model of Propositional Logic 

The Boolean Algebra of the Subsets of a Set 

The Minimization Problem 

The Binary System of Numeration (Appendix I) 
Semiconductor Logic Elements (Appendix IT) 

It would have been helpful if the author had included various alternate nota- 
tions. He uses A for “or’’ and no symbol for ‘‘and.” A short table of notations would 
be helpful to the neophyte, for not ail authors have the thoughtfulness to describe 
their notation. The bibliography is not extensive. 

The printing is by photographic offset process from typed copy, and there was a 
considerable amount of smearing and a number of ghost images in the review copy. 
However, these defects did not make reading seriously difficult. 

This pamphlet should be a handy introduction to Boolean algebra for many 
users and a useful adjunct to the texts for several courses which might be offered in 
colleges. 


C. B. ToMpxKINs 
University of California 


Los Angeles 24, California 


9 [K).—B. M. Bennett & P. Hsu, Significance Tests in a 2 X 2 Contingency Table: 
Extension of Finney-Latscha Tables, July 1960. Deposited in UMT File. 


In testing the significance of deviations from proportionality in a 2 X 2 con- 
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tingency table 
With Attribute Without Attribute Total 
Series I a A-a A 
Series II b B-b B 
Total r N-r N 


with fixed marginal totals A = B, and a/A 2 b/B, Finney (1948, Biometrika 35, 
145) supplied the exact probabilities and the significant values of b at levels .05, 
.025, .01, .005, respectively for A = 3(1) 15. An extension of the same table for A = 
16(1) 20 was prepared by Latscha (1953, Biometrika, 40, 74). The present tables 
further extend the exact probabilities (4 decimal places) to the case: A = 21(1) 25 
and the significant b values at the four levels of the Finney-Latscha tables. 


AvutHors’ SUMMARY 


/10{K].—-Morpecai Ezexret & Kart A. Fox, Methods of Correlation and Regression 


Analysis, 3rd edition, John Wiley & Sons, Inc., New York, 1959, xv + 548 p., 
24 em. Price $10.95. 


It is very pleasing to see that Dr. Ezekiel’s well-known text, Methods of Correla- 
tion Analysis, has been modernized. The older excellent book, which for many years 
stood alone in the desert of statistical literature, has now been joined by a growing 
array of fine statistical text material. Professor Fox was clearly a wise choice for 
Dr. Ezekiel to make for the co-author of the present edition. 

This book is aimed mainly at the non-mathematical reader, in that algebraic and 
computational methods are stressed. There is also a wealth of practical material 
drawn from applied research. However, the mathematical statistician will also find 
much of great value; for example, the use of digital computers for data analysis is 
something which not all statisticians are yet well aware of. 

In general, the treatment of regression and correlation is quite comprehensive, 
and results from the recent theoretical literature have been utilized throughout the 
book. As an example of this point, mention may be made of the inclusion of material 
dealing with regression and the analysis of variance, time series and errors, and the 
fitting of simultaneous relations. 

The book is divided into seven main sections, with twenty-six chapters in all. 
The main sections are: 

. Introductory Concepts 

. Simple Regression, Linear and Curvilinear 

. Multiple Linear Regressions 

. Multiple Curvilinear Regressions 

. Significance of Correlation and Regression Results 

. Miscellaneous Special Regression Methods 

. Uses and Philosophy of Correlation’and Regression Analysis 

There are also three appendices, providing a glossary and important equations, 
methods of computation, and some technical notes. The inclusion of an author index 
is also a commendable feature. 

Those working in the field of economics, agriculture, and business statistics will 
find this text of much value, since a good deal of the material is slanted toward 


SIoo rr Wh 
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those areas. For teachers, statisticians in general, and members of operations re- 
search teams, this revised edition will be found to fill a real need. For applied work 
in fields like engineering, medicine, sociology, and psychology, the full treatment of 
the basic concepts of regression and correlation will be of immense value, and some 
slight reconsideration of the illustrative examples given will often provide insight 
into the real problems in these other fields. 

Thanks are due to both the authors and publishers for making this material 
available. It is a text that should be in the library of every technical organization. 


Harry WEINGARTEN 
Special Projects Office 
Department of the Navy 
Washington, D. C. 


11(K].—E. Frx, J. L. Hopes & E. L. Lenmann, “The restricted chi-square test,” 
included in Probability and Statistics, edited by U. GrenanpeR, Almqvist & 
Wiksell, Stockholm; John Wiley and Sons, New York, 1959, pages 92-108 [See 
the following review]. 


This paper contains some new tables of the power function of chi-square tests, 
i.e., of the non-central chi-square distribution, for small degrees of freedom. As 
is well-known, the power function 8 = (a, f, \) of a chi-square test depends on 
three parameters: a, the “level of significance” at which the test of the null hy- 
pothesis Hy is conducted, i.e., the probability of the test falsely rejecting Hy when 
it is true; f, the “degrees of freedom” of the test; and \, the ‘‘non-centrality param- 
eter,” which measures the “distance” of the alternative H = H(X) under con- 
sideration, from the null hypothesis Hy. The tables of this article give \ to 3D as 
a function of 8 = 0.5(0.1)0.9, 0.95, for f = 1(1)6 and a = 0.001, 0.005, 0.01, 
0.05(0.05)0.3, 0.4, 0.5. The quantity tabulated is that value of the parameter \ 
which satisfied the equation 





—(A 1 +2k-—1 — 2 
2 2k—1 1/2)z 
e am f é a dz 


fo ki! 20/42 (F72) J a) x z= 8 
where f = number of degrees of freedom and x;(a) is such that 


1 4 f-1, —(1/2)22 
a Sa ds = «. 
2ODAT(F/2) Ix jay” , —— 


These tables thus supplement those of E. Fix (1949), reviewed in MTAC, 
v. 4, 1950, p. 206-207. 
CHURCHILL EISENHART 


National Bureau of Standards 
Washington, D. C. 


12(K]—Utr Grenanper, Editor, Probability and Statistics, Almqvist & Wiksell, 
Stockholm; John Wiley & Sons, New York, 1959, 434 p., 24 em. Price, $12.50. 


“Once it had been suggested that a book of studies in probability and statistics 
should be presented to Harald Cramér in honor of his 65th birthday, the authors 
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needed little or no persuasion to contribute.” This volume, subtitled The Harald 
Cramér Volume, is the result. It constitutes a fitting tribute to Harald Cramér, 
Professor of Mathematical Statistics and Actuarial Mathematics at the University 
of Stockholm for over three decades, then President, and now Chancellor of the 
University of Stockholm—Sweden’s outstanding figure in the mathematical theory 
of probability, mathematical statistics, and actuarial mathematics. 

The twenty-one contributors to this volume are from six countries—England 
(2), Finland (1), France (1), India (1), Sweden (3), United States (13); five are for- 
mer students of Professor Cramér; and all have enviable international reputations 
in their respective fields. Their respective contributions are arranged within the vol- 
ume in alphabetical sequence, with page-lengths indicated in parentheses, and are 
as follows: T. W. Anderson, “Some scaling models and estimation procedures in the 
latent class model” (30); M.S. Bartlett, ‘The impact of stochastic process theory on 
statistics” (11); J. L. Doob, “A Markov chain theorem” (8); G. Elfving, “Design 
of linear experiments” (17); W. Feller, “(On combinatorial methods in fluctuation 
theory” (17); E. Fix, J. L. Hodges and E. L. Lehmann, “The restricted chi-square 
test” (16); U. Grenander, “Some non-linear problems in probability theory” (22); 
M. Kac, “Some remarks on stable processes with independent increments’ (9); 
D. G. Kendall, “Unitary dilations of Markov transition operators, and the corre- 
sponding integral representations for transition-probability matrices” (23); P. Levy, 
“Construction du processus de W. Feller et H. P. McKean en partant du mouvement 
Brownien” (13); P. Masani, “Cramér’s theorem on monotone matrix-valued func- 
tions and the Wold decomposition” (15); P. Masani and N. Wiener, “Non-linear 
prediction” (23); J. Neyman, “Optimal asymptotic tests of composite statistical 
hypotheses” (22); H. Robbins, “‘Sequential estimation of the mean of a normal pop- 
ulation” (11); M. Rosenblatt, ‘Statistical analysis of stochastic processes with sta- 
tionary residuals” (30) ; C. O. Segerdahl, “‘A survey of results in the collective theory 
of risk” (24); J. S. Tukey, “An introduction to the measurement of spectra’’ (31); 
8.8. Wilks, ‘“Non-parametric statistical inference’’ (24); H. Wold, “Ends and means 
in econometric model building” (80). Only one of these papers (E. Fix et al) contains 
a mathematical table of general interest, which is considered separately in the im- 
mediately preceding review. 

All in all, this volume provides a panoramic and stimulating view of the work at 
the frontiers of probability and statistical theory and some of their applications. The 
individual scientist, unless he is intimately concerned with the theory of stochastic 
processes, will very likely find only a small fraction of its contents of direct interest 
to him, and, therefore, may consider the volume much too expensive for personal 
acquisition. On the other hand, it is the type of volume that one expects to find in 
the library of a university where research in probability and mathematical statistics 
and their applications is conducted at the post-graduate level, and in the libraries 
of other organizations where research is carried out in the above and related fields. 

Finally, it is to be regretted that the volume does not contain either a photo- 
graph or a biography of Professor Cramér at this milestone in his career. 


CHURCHILL EISENHART 


National Bureau of Standards 
Washington, D. C. 
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13[K].—Masaak1 Srsuya, “Modal intervals for chi-square distributions,” Ann. 
Inst. Statist. Math., v. 9, 1958, p. 225-236. 


Let f(x), —2© <a < ©, bea continuous unimodal probability density function, 
and let L, U(L < U) satisfy the conditions: [ f(z) de = 1 — aand f(U) = 
L 


f(L). Then | L, U | is called the (1 — a)-content modal intervs! for f(z). In this 
paper f(x) is taken to be the x’ density function with ¢ degrees uf freedom. Appli- 
cations are discussed, and a method of computation of modal intervals with given 
content is developed for this case. Tables are presented of L to 4D, U to 3D, 
P.(2’ <= L) to 5D, and (U — L)/(U + L) to 4D, for ¢ = 3(1)30(10)100 and 
1— a= 8, 9, .95, .99. A previous table [1] gave values of 1 — a, effectively for 
(U — L)/U + L) = 01, .05, .10, .20, and @ = 4(4)80. 


C. C. Crate 
University of Michigan 


Ann Arbor, Michigan 


1. M. A. Grrscuicx, H. Rusin & R. SitGreaves, ‘Estimates of bounded relative error 
in particle counting,” Ann. Math. Statist., v. 26, 1955, p. 276-285. 


14[K].—Masaak1 Srsvya & Hipeo Topa, “Tables of the probability density func- 
tion of range in normal samples,” Ann. Inst. Statist. Math., v. 8, 1957, p. 155-165. 


This paper gives details of the calculation of the probability density function 
fn(w), for which there are included tables to 4D corresponding to w = 0(.05)7.65 
and n = 3(1)20, where n represents the size of the normal sample. Cadwell’s 
formula [1] is cited as the basis for this calculation. 


W. J. Drxon 


University of California 
Los Angeles, California 


1. J. H. Capwe t, “The distribution of quasi-ranges in samples from a normal popula- 
tion,’’ Ann. Math. Statist., v. 24, 1953, p. 603-613. 


15[K, 8, V]—O. G. Surron, Mathematics in Action, Harper & Bros., New York, 
1960, xvi + 236 p., 20 cm. Price $1.45 (Paperback). 


In the words of its author, “This book, written primarily for the layman, will 
prove . . . of interest also to students, especially those in the upper forms of schools 
or in the first years in the university. It is a view of the part played by mathematics 
in applied science, as seen by a mathematical physicist.” 

Chapter 1, ‘““The Mathematician and his Task,” begins by discussing the mean- 
ing of theories in physics and the role of mathematics in the development of these 
theories. Chapter 2, “The Tools of the Trade,” gives special attention to complex 
numbers and to the development of the calculus and the differential equations of 
mathematical physics. The approach is essentially that of the physicist (‘‘an infini- 
tesimai quantity [is] one which does not exceed the smallest change of which we can 
take cognizance in our calculations’). 

The remaining five chapters, entitled respectively, “Ballistics or Newtonian Dy- 
namies in War,” “An Essay on Waves,” ‘““‘The Mathematics of Flight,” “Statistics 
or the Weighing of Evidence,” and ‘“‘Mathematics and the Weather,” are essentially 
independent essays that not only provide illustrations of applied mathematics in 














inn. 


tion, 


this 
ppli- 
iven 
3D, 
and 
y for 


error 


func- 
-165. 


ction 
)7.65 
well’s 


ON 


opula- 


York, 


, will 
‘hools 
natics 


nean- 
these 
mplex 
ons of 
infini- 
ve can 


n Dy- 
tistics 
atially 
ties in 











REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 93 


action, but also serve to emphasize fundamentals in both classical and modern 
physics. An especially notable example is the approach to Heisenberg’s «xcertainty 
principle through the problem of the bandwidth required to resolve pulses of short 
duration. Perhaps the most stimulating chapter is the last, in which the author con- 
cludes that “Certain apparently sensible questions, such as the question of weather 
conditions . . . several days ahead, are in principle unanswerable and the most we 
can hope to do is to determine the relative probabilities of different outcomes.” 

There may be some question as to whether the non-mathematical layman will 
be able to follow all of the development of the last five chapters, and the author is 
occasionally guilty of extravagance. (Few aerodynamicists, even in the United King- 
dom, would be willing to admit that F. W. Lanchester’s esoteric volumes Aerody- 
namics and Aerodonetics “played a part in aerodynamics not unlike that exercised 
by Newton’s Principia in astronomy.’’) These things notwithstanding, the reviewer 
believes that the author has succeeded admirably in reaching the goal described in 
the opening quotation of this review. Indeed, he goes beyond this goal, and the book 
(especially the individual essays) is warmly recommended to practicing applied 
mathematicians, as well as to laymen and students. 


Joun W. Mites 
University of California 


Los Angeles, California 


16[K, 8S, W].—C. West Cuurcuman & Puitsurn Ratoosn, Editors, Measurement: 
Definitions and Theories, John Wiley & Sons, Inc., New York, 1959, vii + 274 
p., 24 em. Price $7.95. 


The December 1958 meetings of the American Association for the Advancement 
of Science included a two-day symposium on Measurement. The dozen papers of 
this symposium are included in this volume, together with a related paper sep- 
arately invited but not delivered. 

The book is divided into four parts. Part I, “Some Meanings of Measurement,’’ 
consists of four papers concerned primarily with defining and characterizing the con- 
cept of measurement. Part II, “Some Theories of Measurement,” contains three 
papers that approach the definition and characterization of measurement in more 
formal language—mathematical and logical symbols are more in evidence here—and 
thus warrant grouping together in a separate class. Part III, “Some Problems in the 
Physical Sciences,”’ contains three papers that deal with theoretical and practical 
aspects of measurement in classical and modern physics, plus a paper on “rare 
events” that is out of place in this volume. Part IV, “Some Problems in Social Sci- 
ence,” contains only two papers, one on inconsistency of judgments as a measure of 
psychological distance and one having to do with experimental tests of a probabilis- 
tic theory of economic behavior. 

There are in all fourteen contributors, one paper being co-authored. The disci- 
plines they represent are: philosophy (5), psychology (3), psychophysics (1), physics 
(2), mathematics (1), statistics (1), economics (1), and accounting (1). Philosophy 
is thus somewhat over-represented; astronomy, the biological and earth sciences, 
not at all. 

In the Preface it is stated that the Symposium “was designed to present con- 
trasts in approaches to the problems of measurement.” To this end, “the partici- 
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pants were chosen from different disciplines” and selected particularly “because it 
was known that they had different viewpoints on the meaning and significance of 
measurement.” Consequently, this volume “is a book of contrasts,” not a unified 
collection of interrelated essays on measurement, and definitely not a textbook on 
measurement. At best it presents a broad picture of current thinking on the defini- 
tion, nature, and functions of measurement, against a background of measurement 
needs and practices in various disciplines at the middle of the twentieth century. 
But the “picture” is not all in sharp focus. This is due not so much to differences in 
expository skill of the authors, as to differences in their objectives. Some seek sharp- 
ness of definition at the price of a narrow field of applicability; others demand a 
broad field of applicability at the sacrifice, if necessary, of sharpness of definition. 
Furthermore, much of the discussion is at a level of abstraction so far removed from 
the day-to-day practice of measurement in scientific and industrial laboratories that 
many who have devoted their lives to measurement of the properties of animate and 
inanimate things will find a large fraction of the volume very foreign to them, if not 
entirely unintelligible. Nevertheless, it is a volume that one will expect to find in the 
library of a university or college where research is conducted at the postgraduate 
level. 


CHURCHILL EISENHART 
National Bureau of Standards 


Washington 25, D. C. 


17[K, Z].—JapanesEe STANDARDS AssocrATION, Random Number Generating Icosa- 
hedral Dice (20-face Dice), 6-1 Ginza-higashi, Chuo-ku, Tokyo. Price $2.50 per 
set of 3 dice + postage $.70 (up to 9 sets). 


This device is a set of three icosahedral dice made of plastic material. The dice 
are different colors, red, yellow and blue, so that ordered triplets of digits may be 
generated. Each decimal digit appears on two faces of each die. 

The dice were presumably intended to measure 15 millimeters between parallel 
faces. However, the casting was not particularly good and the measurements listed 
below were recorded between the ten pairs of faces on the new dice tested. This re- 
view was unaccountably lost for several months, and in the intervening period there 
has been considerable flow of the plastic materiai so that the measurements are cur- 
rently considerably worse and actually meaningless, for the faces are clearly no 
longer plane. 


Red Yellow Blue 
14.74 14.90 15.09 
14.97 14.96 14.71 
14.82 14.95 14.93 
14.88 14.97 14.93 
14.94 15.05 14.93 
14.76 14.95 14.83 
14.85 15.07 15.09 
14.96 14.97 14.74 
14.83 14.93 14.87 
14.79 14.90 15.03 


When the dice were new they were tested by 800 rolls each on a level felt surface 
conforming to the specification of ordinary dice tables. Standard tests [1, 2] applied 
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to the individual dice and to the group of three detected no bias. Therefore, it seems 
reasonable to assume that the variations from a true regular icosahedron are minor 
in terms of the application intended. 

Less meticulous samplings have been made with the aged dice, but it seems un- 
likely that a person generating random decimal digits at a rate which can be met by 
these dice would notice any serious bias. 

The frequencies with which the digits appear can be changed slightly by the 
usual standard means. These would include weighting to displace the center of grav- 
ity (an awkward and cumbersome method at best, and, at worst, one which is diffi- 
cult to disguise) or applying wax to one or more faces to increase the probability that 
the waxed face will be on the bottom after the throw. 


C. B. Tompkins 
University of California 


Los Angeles 24, California 
1. THe Ranp Corp., One Million Random Digits and 100,000 Normal Deviates, The Free 


Press, Glencoe, Illinois, 1955. 
2. C. B. Tompxins, RMT 11, MTAC, v. 10, 1956, p. 39-43. 


18[L].—_F. M. Henverson, Elliptic Functions with Complex Arguments, The Univer- 
sity of Michigan Press, Ann Arbor, 1960, v + 38 p. + 160 unnumbered p., 29 
em. Price $8.00. 


The tables and charts on the 160 unnumbered pages are in four equal parts, 


relating in turn to the functions sn w, cn w, dn w, and E(w) = i dn’w dw, each 
0 


for 19 values of the modular angle sin™ k, namely, 1°, 5°(5°)85°, 89°. Each opening 
has a table on the left and a corresponding chart on the right. If w = u + w, 
the tables are all for u/K = 0(.1)1, v/K’ = 0(.1)1, and the charts cover the same 
unit square on a scale such that the side of the square is about 16.2 cm. In Parts 
I-III, the quantities tabulated, to 4 figures without differences, are the real and 
imaginary parts x, y of sn w, cn w, dn w, respectively. In Part IV, if Ey + if; = 
E(w), the relations E(K) = E, E(K + iK’) = E + i(K’ — E’) have led to the 
tabulation of the normalized quantities Ey’ = E,/E, E,’ = E,/(K’' — E’). The 
lines drawn on the charts are curves of constant x or y in Parts I-III, and curves 
of constant E,’ or E,’ in Part IV. Seven-figure values of the complete integrals 
are provided. The information given is sufficient to enable the four functions con- 
cerned to be evaluated for any point w in the complex plane. 

The well-known tables of Spenceley and Spenceley [1] were used as a source 
of values for real w, whence the imaginary transformation and the addition formu- 
las were used to compute the values of the functions of u + i. The computation 
was mostly done on an IBM 650. 

The Introduction contains enough information about elliptic integrals and 
functions to explain the tables and charts to anyone not previously acquainted 
with the subject. It also contains several applications to potential problems. It is 
pleasant to find such a valuable contribution to mathematical tabulation made 
by a civil engineer. The charts were constructed by a group of five Turkish naval 
officers at the University of Michigan. 

One could wish that italic type (available and used in other contexts) had been 
used for mathematical symbols in the Introduction, but no such minor matter 
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can obscure the importance of the volume, which cannot fail to be found very 
useful. Only in the case of Part IV has the reviewer heard of any similar or related 
table of comparable scope, namely the table of the Jacobian zeta function by Fox 
and McNamee [2]. 
A. F. 
1. G. W. Spencetzey & R. M. Srencetey, Smithsonian Elliptic Functions Tables, Wash- 
ington, 1947. See MT'AC, v. 3, 1948, p. 89-92. 


2. E.N. Fox & J. McNameg, “The two-dimensional problem of seepage into a cofferdam,”’ 
Phil. Mag., 8.7, v. 39, 1948, p. 165-203. See MT'AC, v. 3, 1948, p. 246, 252, also v. 7, 1953, p. 190. 


19[L].—L. Lewin, Dilogarithms and Associated Functions, MacDonald & Co. Ltd., 
London, 1958, xvi + 353 p., 21 em. Price 65 Shillings. 


The functions treated here are for the most part special cases of the Lerch zeta 
function, which can be defined by the series } 70 2”/(n + b)’, |z| <1, b nota 
negative integer or zero. To describe the text and tables, it is convenient to give 
some notation. Let z = x + iy, where x and y are real and i = »/—1. Then 


(1) li.(z) = — [ t* In (1 — t) dt, Li,(z) = [ t*Lins(t) dt; 


Tin(z,a) = / (t + a)~ arctan é dt, Ti(z) = Tin(x,0), 
0 


(2) 
Lin(iy) = 2-*Li,(— y’) + iTin(y); 
8 8 
Cl.(6) = - { In (2 sin} t) dt, Clea(0) = i Clons(t) dt, 
0 2 0 
(3) é 
Clon41(0) -_ Liens: (1) sa [ Clon(z) dt; 
C) ei”? C) ef 
(4) Glan(0) + iClon(0) = 2), Clansa (8) + iGlansa(@) = 20 
k=1 k=1 


The function Gl,(@) is a polynomial in @ of degree n. 

Chapter I deals with the dilogarithm function Ii,(z). The function T7i2(x) is 
considered in Chapter II; Ti2(z, a) in Chapter III. Ci.(@), @ real and positive, is 
Clausen’s integral, and is studied in Chapter IV. Chapters V and VI take up 1i,(z) 
for n = 2 and 3, respectively, and the analysis of this function for general values of 
n is the subject of Chapter VII. The general relations in (3) and (4) are also studied 
in this chapter. Chapter VIII deals with series expansions and integrals which can 
be expressed in terms of the basic functions in (1)—(4). Chapter [X is very useful. 
It is a compendium of results derived in the previous chapters. It also contains a 
survey of mathematical tables. 

A description of the tabular material in this volume follows. 


Table I. Lin(x), n = 2(1)5 
z=0(.01)1.0, 5D 
Table II. Tin(y), n = 2(1)5 


y = 0(.01)1.0, 5D 
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Table III. Cl,(3xa), mn = 2(1)5 
0(.01)2.0, 5D 
= 2(1)5 
a = 0(.01)20, 5D 
Table V. Li2(r, @) = real part of Li(z), z= re”. 
r=0(.01)1.0, 6 =0(5°)180°, 6D. 


2 
ll 


Table IV. Gl,(42a), 


3 


Throughout Table V the symbol z should be replaced by r. No information is 
supplied for interpolating in the tables. 

The volume is replete with striking and curious results, some of which have been 
rediscovered a number of times, and the book should prevent future duplication of 
effort. There is a well-detailed table of contents and index. An extensive bibliography 
is also given. 


Y.L.L. 


/ 20(L].—Vera I. Pacurova, Tablitsy integro-eksponentsial’not funkisii 


E,(z) = [ eu “du 


(rabies of the Exponential Integral Function E(x) = f eu” au), Akad. 
1 


Nauk SSSR, Vychislitel’nyy Tsentr, Moscow, 1959, xii + 152 p., 27 cm. Price 
9.60 rubles. 


This volume from the Computational Center of the Academy of Sciences of the 
USSR deals with well-known integrals which depend on the exponential integral 
when » is a positive integer n. There are three tables. Table I (pages 3-52) is re- 
produced, with acknowledgment, from the NBS table calculated for a report of 
1946 by G. Placzek and Gertrude Blanch (see MTAC, v. 2, 1947, p. 272) and more 
widely disseminated in 1954 in [1]. The table gives Z,(z) to 7 or more decimals 
for n = 0(1)20, x = 0(.01)2(.1)10, and also 7-decimal values of the auxiliary 
functions E,(z) — x ln x and E;(x) + 32° Inz for x = 0(.01).5 and x = 0(.01).1, 
respectively; these last two ranges need transposing in the sub-title on page 1. 

The other two tables are original. It is not stated what machines were used in 
computing them. Table II (pages 54-62) gives e*E, (x), n = 2(1)10 to 7 decimals 
(6 figures) and e* to 7 figures, all for z = 10(.1)20. Table III (pages 64-151) 
gives eE,(x), v = 0(.1)1 to 6 or 7 figures and e* to 7 figures, all for z = 
.01(.01)7(.05)12(.1)20. No table gives differences. 

A short introduction contains mathematical formulas and recommendations 
about interpolation. For integral n, the formula d'E,(x)/dx" = (—1)'E,_,(2) 
enables the tabulated function values themselves to be used for interpolation by 
means of Taylor’s series. A table is given showing the accuracy attainable in inter- 
polating various functions linearly or with 3 or 4 Taylor terms or with 3, 4, or 5 
Lagrange terms. ; 


A. F. 


1. Nat. Bur. Stanparps, Appl. Math. Ser. No. 37, Tables of Functions and of Zeroes of 
Functions, U. S. Government Printing Office, Washington, D. C., 1954, p. 57-111. See RMT 
104, MTAC, v. 10, 1956, p. 249. 
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21(L).—Eart D. Raryvituz, Special Functions, The Macmillan Co., New York, 
1960, xii + 365 p., 24 cm. Price $11.75. 


This aptly titled, interesting, extremely well written book is based upon the lec- 
tures on Special Functions given by the author at the University of Michigan since 
1946. The author’s aim in writing the book was to facilitate the teaching of courses 
on the subject elsewhere. As an instructor in such a course, the reviewer feels certain 
that this most welcome text is to be accorded a warm reception on many college 
campuses. 

More than fifty special functions receive varying degrees of attention; but, for 
the sake of usefulness, the subject is not approached on the encyclopedic level. 
Many of the standard concepts and methods which are useful in the detailed study 
of special functions are included. There is a great deal of emphasis on one of the 
author’s favorite subjects: generating functions. Two interesting innovations are 
I. M. Sheffer’s classification of polynomial sets and Sister M. Celine Fasenmyer’s 
technique for obtaining recurrence relations for sets of polynomials. Functions of 
the hypergeometric family hold the center of the stage throughout a major portion 
of the text. The book concludes with a short current bibliography which should 
enable the reader to begin a more detailed study of the fielc. 

There are twenty-one chapters in all, and the book may be roughly divided into 
four distinct parts: 

(1) Two short preliminary chapters, 1 and 3, deal separately with infinite prod- 
ucts and asymptotic series, respectively. 

(2) Chapter 2 treats the gamma and beta functions, and chapters 4, 5, 6, and 7 
are devoted to the hypergeometric family: the hypergeometric function, generalized 
hypergeometric functions, Bessel functions, and the confluent hypergeometric func- 
tion, respectively. 

(3) Chapter 8 is concerned with the generating function concept, as a prepara- 
tion to chapters 9, 10, 11, and 12, which consider orthogonal, Legendre, Hermite, 
and Laguerre polynomials, respectively. Chapter 13 contains I. M. Sheffer’s classi- 
fication of polynomial sets; chapter 14 contains Sister M. Celine Fasenmyer’s tech- 
nique for obtaining recurrence relations for polynomials; and chapter 15 contains 
symbolic relations among classical polynomials. There follow three polynomial 
chapters: 16, 17, and 18, on Jacobi, ultraspherical and Gegenbauer, and other poly- 
nomials, respectively. 

(4) The concluding three chapters, 19, 20, and 21, are devoted to elliptic func- 
tions, theta functions, and Jacobian elliptic functions, respectively. 


J.B. Diaz 
University of Maryland 


College Park, Maryland 


22[L].—L. J. Starer, Confluent Hypergeometric Functions, Cambridge University 
Press, New York, 1960, x + 247 p., 29 cm. Price $12.50. 


This is a valuable treatise on the subject, and the first of its kind in English. 
Tricomi [1] and Buchholz [2] have previously written books on the subject in Italian 
and German, respectively. Tricomi used the notation derived from the theory of 
hypergeometric functions, and this also is the principal notation employed in the 
Bateman Manuscript Project [3]. Buchholz uses the notation introduced by Whit- 
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taker. The notation of the present volume is a fusion of the two, and many results 
of the same character are given in all of the notations. This should prove most useful 
to research workers, as the confluent function has many applications, and notation 
is not uniform. 

Chapter I studies the confluent hypergeometric functions as solutions of dif- 
ferential equations. Power series expansions are developed and relations between 
the various functions are carefully detailed. 

Chapter II deals with differential properties including contiguous relations, 
Wronskians, addition theorems and multiplication theorems. 

Chapter III is principally concerned with definite integrals involving the con- 
fluent functions. These include integrals of Barnes, Euler and Pochhammer types, 
and Laplace, Mellin and Hankel transforms. 

Chapter IV takes up asymptotic expansions. It is a very useful compendium. 
In particular, as the author remarks, the theory is not complete, and it is in this 
area where important new results are anticipated. 

Related functions such as Coulomb wave functions, Bessel functions, incomplete 
gamma functions, etc., are briefly considered in Chapter V. 

Descriptive properties are the subject of Chapter V1. Zeros in x of ,F;(a; b; x) 
and W;,_(x), and zeros in a and b of the former are studied. Formulas and expan- 
sions for the zeros are provided. Numerical evaluation of ,/’; is discussed and nu- 
merous figures are provided to illustrate its behavior when a, b, and z are real. 

Three appendices give tables. These are as follows: 


Appendix I. Smallest positive zeros of ,Fi(a; 6; x) 
a = —4.0(0.1) — 0.1, b = 0.1(0.1)2.5, 7D 

Appendix II. F(a; b; x) 
a = —1.0(0.1)1.0, b = 0.1(0.1)1.0, x = 0.1(0.1)10.0, 
78 — 98 

Appendix III. ,F (a; b; x) 


a = —11.0(0.2)2.0, b = —4.0(0.2)1.0, z=1, 7D 


Appendix I was calculated to 8D on EDSAC, and the error does not exceed two units 
in the seventh decimal. Appendix II was calculated on the same machine, but no 
mention is made of accuracy. Some spot checks indicate that the error does not 
exceed four units in the last figure recorded. Appendices I and II are extensions of 
tables previously published by the author. Other short tables of the material in 
Appendix II have also appeared previously, but the present table is the most com- 
plete available. Appendix III was evolved using recurrence formulas and the data 
of Appendix II. It is quite surprising that nothing is said about interpolation. 

A detailed table of contents, general index, and symbolic index of notation en- 
hance the usefulness of the volume. The list of references, though not exhaustive, 
is fairly complete for work subsequent to that of Buchholz [2], whom the reader 
should consult for an extensive bibliography of work prior to about 1952. The 
reviewer notes that two important tables of Whittaker functions reviewed in [4] 
are not found in Miss Slater’s list of references. 


Y. L. L. 
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1. F. G. Tricom1, Lezioni sulle funzioni ipergeometriche confluenti, Gheroni, Torino, 1952. 
. 2, H. Bucnnouz, Die konfluente hypergeometrische Funktion mit besonderer Beriicksichtigung 
threr Anwendungen, Springer, Berlin, 1953. 
' 3. A. Erp&éiy1, W. Maenus, F. Opernettincer, & F. G. Tricomi, Higher Transcendental 
Functions, v. 1, Chapter VI, McGraw-Hill, New York, 1953. 

4. RMT 46, MTAC, v. 12, 1958, pp. 86-88. 


23[M, X].—G. Doetscu, Einfuhrung in Theorie und Anwendung der Laplace-Trans- 
formation, Birkhauser Verlag, Basel, Switzerland, 1958, 301 p., 24 cm. Price 
SFr 39.40. 


This book forms an excellent introduction to the subject of Laplace transforma- 
tions in one dimension, written by one of the leading experts in the field. From his 
wide knowledge of both the theoretical and applied aspects of the subject, the author 
has written a very readable, rigorous exposition which also indicates relations with 
appropriate physical concepts. 

After a general introduction, the basic properties of the Laplace transform are 
developed in ten short chapters. Included are discussions of half-planes of conver- 
gence, uniqueness of the inverse, analytic properties of the transform, the effect of a 
linear transformation of the independent variable on the transform, the effect of dif- 
ferentiation and integration, and the transformation of a convolution. Because each 
chapter is devoted to a separate topic, the book is very useful for reference purposes. 
Many examples of specific transforms are given. 

The next four chapters deal with the application of Laplace transforms to the 
following problems: the initial-value problem in ordinary differential equations with 
constant coefficients; the solution of differential equations for special input func- 
tions; homogeneous and non-homogeneous systems of differential equations; the 
initial-value problem for difference equations. 

The next group of chapters deals with further properties of the transform, some 
of which may not be familiar to the reader: the behavior of the transform at infinity; 
inversion formulas expressed as integrals along vertical lines, as integrals along de- 
formed paths in the complex plane, and as series of residues; conditions for the rep- 
resentation of a function as a transform; functions given as the sums of series of 
transforms; the analogue of Parseval’s formula and transforms of products; and the 
asymptotic behavior of the transform and of the “original” function. 

The book concludes with three chapters on further applications of the Laplace 
transform to differential equations with variable coefficients, to simple partial dif- 
ferential equations, and to certain integral equations. 

A useful feature of the book is the inclusion of necessary background material, 
particularly in the later chapters. 


Dororuy L. BERNSTEIN 
Goucher College 


Towson, Baltimore 4, Maryland 


24[S, X]—G. I. Marcuux, Numerical Methods for Nuclear Reactor Calculations, 
(An English translation of a work originally published in Russian as Supplement 
Nos. 3-4 of the Soviet Journal of Atomic Energy, Atomnaia Energiia. Atomic 
Press, Moscow, 1958.) Consultants Bureau, Inc., New York, 1959, 293 p. 28 cm. 
Price $60.00. 


A more accurate title for this book would have been Numerical Methods for 
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Nuclear Reactor Physics Calculations. As stated by the author, “the book is an at- 
tempt at a more or less systematic exposition of numerical methods for the calcula- 
tion of thermal, intermediate, and fast neutron reactors. Particular attention is 
devoted to the problems of critical mass, the space-energy neutron flux distribution, 
and the neutron importance.” This is the first book to become available in English 
that is entirely devoted to this area of calculation, and as such is a helpful addition 
to our literature. 

The book is clearly not addressed to mathematicians. A familiarity with nuclear 
reactor theory is assumed; terms such as “cross section” are used without definition. 
The mathematician will probably be disturbed by such statements as “the inte- 
grands in the resulting two terms are now continuous, and they can therefore be ex- 
panded in a Taylor’s series about r,” (p. 106), or, “primarily, this is due to the fact 
that a second-order difference equation has two eigenvalues, one of which is posi- 
tive” (p. 122). Emphasis throughout the book is on the derivation of equations and 
the mechanics of their solution. The stability of a number of the numerical methods 
against the exponential growth of round-off errors is considered. Problems related 
to the existence of solutions and the accuracy with which solutions of the various 
approximate equations represent the desired solutions are mentioned only occa- 
sionally. 

In the foreword the author mentions that the numerical-methods were ‘first 
tested on a large quantity of theoretical and experimental data and were later used 
in operation.” It is unfortunate that the results of some of these calculations were 
not included in the book. Results of specific calculations are not given. For example, 
in Chapter VIII, various sets of finite-difference equations are set up to approximate 
the diffusion equation. Some of these difference equations are more accurate than 
others. However, no discussion of the actual errors incurred when these various 
approximations are used is presented other than a statement, “with a sufficiently 
fine network good accuracy can be obtained.” This would have been an excellent 
place to include numerical results showing the kind of accuracy that could be at- 
tained when the various approximate equations are used in typical specific problems. 

A useful feature of this book is that corresponding to every set of equations de- 
rived, an adjoint set of equations is also obtained. 

In the first chapter a basic set of transport equations is derived: a “slowing- 
down” equation in which the assumption is made that the moderator nuclei are at 
rest while the neutrons are slowed down by elastic collisions, and a “thermal” 
equation in which the velocities of the moderator nuclei are assumed to have a 
Maxwellian distribution. In succeeding chapters, various sets of equations ap- 
proximating this basic set are obtained. In the second chapter the diffusion approxi- 
mation is derived following the method of R. E. Marshak, H. Brooks, and H. 
Hurwitz, Jr., published in Nucleonics, v. 4, 5 (1949), in which the angular de- 
pendence of the neutron flux and the collision function are both approximated by 
retaining the first two terms in a series expansion involving spherical functions. 
The corresponding boundary conditions are also derived. In the third chapter the 
diffusion-age approximation is obtained when the energy dependence of the flux is 
approximated by retaining the first few terms of a Taylor series expansion. Here 
again, the corresponding boundary conditions are derived. In the fourth chapter 
further approximations to the basic transport equation are derived which take 
more accurate account of energy dependence than the diffusion-age approximation, 
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and are therefore applicable in cases where there is stronger absorption of slowing- 
down neutrons. Equations are derived in which resonance absorption is accounted 
for based on a method suggested by Wigner. Also, the Greuling-Goertzel equations 
for the slowing-down of neutrons in an inhomogeneous medium are obtained. 

In Chapter V exact solutions are given of the basic and adjoint equations of the 
diffusion-age approximation for a homogeneous, unreflected reactor with an ex- 
trapolated boundary which is energy-independent. Also, for this special case the con- 
vergence of “the method of successive approximations” for finding the criticality 
factor is established. ‘“The method of successive approximations” means here that 
at each step, assuming an estimate of the fission source, the corresponding neutron 
flux is calculated, and then using this set of flux values, a new estimate of the fission 
source is obtained. The criticality factor is obtained using ratios of successive source 
estimates. 

In Chapter VI various sets of multigroup equations are derived which are appli- 
cable primarily to thermal reactors where weak absorption of slowing-down neu- 
trons is assumed; in Chapter VII various sets of multigroup equations are derived 
which are applicable primarily to intermediate reactors where strong absorption of 
slowing-down neutrons must be considered. 

Chapter VIII contains the derivation of various sets of three-point difference 
equations to approximate the one-dimensional diffusion equation in plane, cylin- 
drical, and spherical geometries. In addition, five-point difference equations are 
obtained for the two-dimensional diffusion equation in plane (2x, y) and cylindrical 
(r, 2) geometries. 

In Chapter [X various methods are discussed for solving sets of finite-difference 
diffusion equations of the types set up in the previous chapter. To solve the three- 
point difference equation a well-known factorization technique is described, and its 
stability against the growth of round-off errors is considered. For the solution of the 
five-point difference equations an analogous factorization technique, as well as a 
simple iteration procedure, and a third method suggested by N. I. Buleev are con- 
sidered. Error criteria which can be used to estimate the accuracy of results are 
also mentioned. 

Chapter X contains a discussion of the “net”? method for obtaining approximate 
solutions of the differential slowing-down equation in one spatial dimension. Dif- 
ference equations are set up to approximate the differential equation and solved 
for the unknown function in one energy interval in terms of values computed for 
the preceding interval. Both explicit and implicit schemes are discussed. Their 
stability is investigated by applying the von Neumann stability criterion. 

In Chapter XI perturbation theory is applied to obtain approximate formulas 
giving the change in the criticality factor when small changes are made in values of 
the various physical parameters. This is done starting with the basic equations of 
Chapter I as well as with the various equations obtained in later chapters to ap- 
proximate these basic equations. Some further applications of perturbation theory 
are discussed such as the calculation of the reactivity equivalent of a control rod 
in the center of the core of a cylindrical reactor, as well as for a whole system of con- 
trol rods in such a reactor. 

Chapter XII on heterogeneous effects in nuclear reactors contains descriptions of 
the Gurevich-Pomeranchuk theory of resonance absorption applicable in calculat- 
ing the effective resonance integral for “thin” U™ blocks, the Wigner theory for 
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calculating the resonance integral for “thick” blocks, and Orlov’s extension of Wig- 
ner’s method giving a formula for the effective resonance integral which can be used 
for U™ blocks of any thickness. Further, a method is given which was developed by 
Marchuk and Orlov for calculating the resonance capture of neutrons in a plane 
lattice of uranium blocks. Sedel’nikov’s extension of the Gurevich-Pomeranchuk 
theory to take account of self-shielding in intermediate neutron reactors is de- 
scribed. Effective boundary conditions for “black” and “grey” bodies are con- 
sidered. Finally, the Wigner-Seitz method is applied to a cell of a heterogeneous 
reactor to obtain effective constants for an equivalent homogeneous reactor. 

Chapter XIII contains a discussion of the calculation of the neutron flux in fast 
nuclear reactors. The transport equation including terms for inelastic scattering is 
given. Corresponding multigroup equations are obtained. For solving these equa- 
tions both the method of spherical harmonics and Carlson’s S, method are de- 
scribed. 

A final chapter on calculations for intermediate and thermal reactors with hy- 
drogenous moderators is included. This chapter was added after the original manu- 
script had been prepared for publication. 

It is unfortunate that this book was not more carefully edited. For example, 
there is confusion in the use of the terms ‘flux’ and ‘current’; on page 104 the follow- 
ing statement appears, “‘we wish to find a solution of (30.1) which has a continuous 
flux 


ap a6 ,, 
I = r°D 

Throughout the text, the equations which are usually referred to as “transport” 
equations are called “kinetic” equations. Too frequently proper names are spelled 
phonetically as a result of the transliteration process, rather than with their usual 
English spelling; for example, R. E. Marshak appears as R. E. Marchak (p. 6), 
R. Ehrlich appears as R. Erlich (p. 6), Neumann appears as Neiman (p. 142), 
etc. Also, it is unfortunate that the publishers did not take more care in their repre- 
sentation of symbols. The equations were apparently reproduced photographically. 
A parameter represented by a script letter in an equation is frequently represented 
in another form in the text. On pages 132 and 134 a parameter which appears in 
the equations as a “‘chi’’ is represented as a “kappa”’ in the text. 

This book is a very useful addition to our literature. It is hoped that its presence 
and imperfections will act as a stimulus for the publication of another book in the 
area of numerical methods for nuclear reactor calculations which will give a more 
satisfying mathematical treatment of this subject, and which will be made available 
at a more reasonable price. 

E. CuTuinn 
Applied Mathematics Laboratory 


David Taylor Model Basin 
Washington 7, District of Columbia 


25(W].—E. Burcer, Einfuhrung in die Theorie der Spiele, Walter de Gruyter & Co., 
Berlin, 1959, 169 p., 23 em. Price DM 28. 


This is the first book on the theory of games to appear in the German language. 
The author, professor of mathematics at the University of Frankfurt, has previously 
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given several courses on the subject for economists, but the emphasis in this book 
is to bring out the mathematical structure of the theory. He has succeeded in this 
in an admirable way. The comparatively few applications are meant as illustrations 
rather than as developments of the empirical fields for which the theory of games is 
particularly suited. There is an extensive discussion of zero-sum two-person games, 
and of equilibrium points for non-cooperative games. The presentation of linear pro- 
gramming is also clear and rather exhaustive. The theory of cooperative games, that 
is, those in which the formation of coalitions is advantageous and allowable, and in 
which side-payments by the players are freely admitted, is developed in fair detail, 
even including a discussion of Shapley’s value of an n-person game. At this point the 
author confesses that he is less sure of the intuitive background against which this 
theory has been placed, a position that is common to many mathematicians who 
have studied the problem of n-person games. However, this is a difficult issue and 
the author is wise not to have taken too definite a stand, rather withdrawing to the 
strictly mathematical aspects involved. These problems can only be solved by re- 
course to an improved description of the socio-economic world. If it turns out that 
the real problem involves a high degree of cooperation—and I have no doubts what- 
soever that this will be the case—then the mathematical theory will have to accom- 
modate itself to these facts, even if the mathematical structure is uncommon and 
cumbersome, until fundamentally new concepts are established. 

Dr. Burger possesses a very high didactic skill; the great clarity which pervades 
his whole book should make it a welcome tool for the novice in game theory who 
commands the mathematical knowledge expected of first or second year graduate 
students. A translation of the work should be seriously considered, since there is no 
similar book in English which accomplishes as much in such small compass as Dr. 
Burger’s does. 


OsKAR MORGENSTERN 
Princeton University 


Princeton, N. J. 


26[X].—R. P. Boas, Jr. & R. C. Buck, Polynomial Expansions of Analytic Func- 

tions, Springer-Verlag, Berlin, 1958, viii + 77 p., 23 cm. Price DM 19.80. 

A great part of numerical analysis is concerned with polynomial approximation 
to analytic functions, and so this booklet appears of immediate interest to the nu- 
merical analyst. However, with numerical analysis in its present state, it is more 
relevant for studies in the general theory of functions of a complex variable or in 
the theory of special functions. 


Given a set of polynomials {p,(z)} and a function f(z), it is reasonable to ask 
whether we can find coefficients {c,} such that 


(1) f(z) = Lo enpn(z) 


in some sense. This is the “expansion problem’’. 
It is also reasonable to ask whether, given linear functionals { L,,(f)}, for example, 


(2) L,(f) = f(0) 


or 


(3) Lon(f) = f°"(0), — Lensa(f) = f°" (1), 
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we can obtain a representation of f in the form (1). This is the “interpolation prob- 
lem”. The solution in case (2) is given by the Taylor expansion and that in (3) 
by the two-point Lidstone expansion 


f(z) = DL An(l — z)f°" (0) + DY An(z)f°"(1), 
where the polynomials A,(z) are defined by 


sinh zw 

sinh w 
This is valid, for instance, when f(z) is an entire function of exponential type less 
than z. 

These two problems are studied by a general method—‘“kernel expansion” — 
for wide classes of polynomials (defined by generating functions) and for the cases 
when f(z) is entire, or regular at the origin. 

The material is accessible to those familiar with the classical methods of complex 
variable theory, and its study by numerical analysts is recommended. It will, 
for instance, encourage us to get off the real axis, reveal some thought-provoking 
“bad examples” (e.g., non-uniqueness in (1)), and show us how mathematics 
should be written. 





= > A,(z)w™. 


Joun Topp 
California Institute of Technology 
Pasadena, California 


27(X|—RvupcuprH E. Lancer, Editor, Symposium on Numerical Approximation, 
The University of Wisconsin Press, Madison, 1959, x + 462 p., 24 cm. Price 
$4.50. 


The present volume auspiciously launches the publication efforts of Professor 
Langer’s group at the Mathematics Research Center. It contains the proceedings of 
a symposium conducted at the University of Wisconsin, April 21-23, 1958. The 
papers (original and expository) present such a well-rounded survey of numerical 
approximation because of the careful selection of the invited contributors. A listing 
of the authors and the titles of their papers will serve to indicate the level and scope 
of the work. 

A. M. Ostrowski On Trends and Problems in Numerical Approximation 


R. C. Buck Linear Spaces and Approximation Theory 

Z. Kopal Operational Methods in Numerical Analysis Based on Ra- 
tional Approximations 

P. J. Davis On the Numerical Integration of Periodic Analytic Func- 
tions 

H. E. Salzer Some New Divided Difference Algorithms for Two Vari- 
ables 

P. C. Hammer Numerical Evaluation of Multiple Integrals 

M. Golomb Optimal Approximation and Error Bounds 

A. Sard The Rationale of Approximation 

J. L. Walsh On Extremal Approximations 

E. L. Stiefel Numerical Methods of Tchebycheff Approximation 


L. Fox Minimax Methods in Table Construction 
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T. S. Motzkin Existence of Essentially Nonlinear Families Suitable for 
Oscillatory Approximation 


I. J. Schoenberg On Variation Diminishing Approximation Methods 

M. Golomb Approximation by Functions of Fewer Variables 

J.C. P. Miller Extremal Approximations—A Summary 

R. C. Buck Survey of Recent Russian Literature on Approximation 

F. L. Bauer The Quotient-Difference and Epsilon Algorithms 

J. B. Rosser Some Sufficient Conditions for the Existence of an Asymp- 
totic Formula or an Asymptotic Expansion 

J. W. Tukey The Estimation of (Power) Spectra and Related Quantities 

L. Collatz Approximation in Partial Differential Equations 

J. Todd Special Polynomials in Numerical Analysis 


E. I. 


28[X].—Rvupoitex E. Langer, Editor, Frontiers of Numerical Mathematics, The 
University of Wisconsin Press, Madison, 1960, xi + 132 p., 24 cm. Price $3.50. 


This book contains eight papers and a discussion of these papers presented at a 
Symposium held at Madison, Wisconsin on October 30-31, 1959. This symposium 
was conducted jointly by the Mathematics Research Center, the United States 
Army, and the National Bureau of Standards. Its purpose “was not intended to be 
an occasion for the presentation of research results, but one for a survey of the fu- 
ture; for the identification of some mathematical problems that will have to be faced 
in the lines of scientific advance.” 

The authors and the titles of their papers are: 


William Prager Stress Analysis in the Plastic Range 

Garrett Birkhoff Some Mathematical Problems of Nuclear Reactor 
Theory 

Zdenek Kopal Numerical Problems of Contemporary Celestial Me- 
chanics 

Lee Arnold Aeroelasticity 

Phillip M. Morse Operations Research 

Joseph O. Hirschfelder Mathematical Bottlenecks in Theoretical Chemistry 

S. Chandrasekhar Magnetohydrodynamics 

J. Smagorinsky On the Application of Numerical Methods to the 


Solution of Systems of Partial Differential Equa- 
tions Arising in Meteorology 

It is unfortunate that the paper by Lee Arnold was made from a tape recording 
and was not finally reviewed by the author. It is very difficult for a reader to follow 
this paper, for many statements in it seem to refer to illustrations which are not in- 
cluded. 

The reader of this book should not expect to find a detailed discussion of numer- 
ically formulated problems arising in various branches of science. He will find, in the 
main, discussions and reviews of various open mathematical problems in different 
scientific areas. Some of the papers refer to numerical treatment of these problems. 
A surprising number indicate a strong preference on the author’s part for analytical 
methods for dealing with their problems. 


One can but agree with J. Smagorinsky when he states “that computing ma- 
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chines cannot be considered a substitute for the ingenious mathematical and labora- 
tory techniques of analysis which have been devised. . . .”” However, this reviewer 
is convinced that such ingenuity when properly coupled with the power of modern 
computers will provide greater insight than analytical methods alone. 

A. H. T. 


/ 29[X|.—(a) D. S. Mrrrinovié, “Sur les nombres de Stirling de premiére espéce et 


les polynémes de Stirling,” Publ. de la Fac. d’ Electrotechnique de I’ Univ. a Bel- 
grade (Série: Math. et Phys.), No. 23, 1959, 20 p. (Serbian with French summary. 
Tables by Miss RuZica S. Mirrinovicé.) 

(b) D. S. Mrrrinovié & R. 8. Mrrrinovié, “Sur les polynémes de Stirling,” 
Bull. Soc. Math. Phys. Serbie, v. 10 (for 1958), p. 43-49, Belgrade. (Summary in 
Russian.) 

(c) D. 8S. Mirrinovié & R. 8. Mirrinovié, “Tableaux qui fournissent des 
polynémes de Stirling,” Publ. Fac. Elect. Univ. Belgrade (Série: Math. et Phys.), 
No. 34, 1960, 24 p. (Summary in Serbian.) 

These three papers are concerned with the Stirling numbers of the first kind, S,’, 
which may be defined for positive integral n by 


> S,'2'. 


r=0 


a(x — 1)(a — 2) --- (a@ —n +1) 


Altogether the numbers S,~” are tabulated for m = 1(1)32, n = m + 1(1)N, 
where N = 200 for m = 1(1)5, N = 100 form = 6, and N = 50 form = 7(1)32. 
The values for m = 1(1)7 are given in (a), for m = 8(1)13 partly in (a) and 
partly in (c), for m = 14(1)20 partly in (b) and partly in (c), and for m = 
21(1)32 in (c). The authors found no discrepancy as a result of some checking 
against unpublished tables by F. L. Miksa (see MTAC, v. 10, 1956, p. 37). 
Algebraic expressions for Si” in the form of binomial coefficients a - 1 
multiplied by polynomials (with factors n(n — 1) separated out if m is odd and 
not jess than 3) are given for m = 1(1)13 in (a) and for m = 1(1)9 in (ce). 
Si” may also be expressed as a sum of multiples of binomial coefficients in the 


form 
m—1 
n—m y & n 
ss" = d re a a ‘): 


Altogether the values of the coefficients C,,." are given for k = 0(1)31, m = 
k + 1(1)32, the values for k = 0(1)19, m = k + 1(1)20 being found in (b) 
and the remaining values in (c). 

A. F. 


30|Z].—_ Wayne C. Irwin, Digital Computer Principles, Van Nostrand Co., Inc., 
Princeton, 1960, vi + 321 p., 24 em., $8.00. 


This book contains material presented at a training course in the Electronics 
Division of the National Cash Register Company. It is an extremely elementary 
“book for the beginner. No previous acquaintance with computers, electronics or 
mathematics is necessary.” 
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The book is divided into nine sections and forty-two chapters. The section titles 
are: Methods of Computation, Symbolic Logic, Mechanization of Logic, Mechani- 
zation of Storage, Timing, Mechanization of Arithmetic, Control, Communication 
with the Computer, Preparation of Instructions, Reduction of Errors, and Present 
Trends. 

The author has placed more emphasis on symbolic logic than on any other sub- 
ject. This section of the book begins with Boolean algebra and terminates with a 
discussion of the Harvard Minimizing Chart. However, even in this section, the best 
written one in the book, the author does not develop the mathematical ideas in- 
volved adequately or prove statements. He describes various facts and illustrates 
the use of various techniques. 

The author does not go very far in the discussion of any topic he selects, and he 
omits many significant points in his discussion. Some omissions are noted in foot- 
notes but some are never mentioned. For example, “asynchronous machines’ are 
dismissed with a footnote on page 145. The words “round-off” and ‘“‘rounded multi- 
plication” are never mentioned in the text and neither word can be found in the 
index. The existence of a computer using a number representation other than one 
with a signed absolute value is never mentioned. 

It is unfortunate that the author did not see fit to include in the bibliography 
references to the reports by von Neumann and his co-workers at the Institute for 
Advanced Study. The discussion given in these reports of arithmetic performed in a 
computer with a two’s complement representation of numbers is clear and complete, 
and would supplement very nicely the author’s limited discussion of arithmetic in 
a binary machine. 


A. H. T. 


31(Z].—T. E. Ivau1, Electronic Computers, Principles and Applications, Second Edi- 
tion, Philosophical Library, New York, 1960, viii + 263 p., 22 em. Price $15.00. 


This is a revised version of the first edition published in 1956 under the same title 
[1]. In the first edition the author acts as an editor, publishing in book form the sep- 
arate contributions of a number of English writers working in the field of computers. 
In the present edition the author attempts to regroup and rewrite the material so 
as to present a more coherent picture. He also adds three short chapters, as follows: 

(1) In place of one chapter on Analogue Computing Circuits the new edition 
now contains two chapters, Analogue Computing Circuits—1, and Analogue Com- 
puting Circuits—2. 

(2) A chapter on Programming Digital Computers is added. 

(3) The chapter, Computers of the Future, somewhat rewritten, now becomes 
Recent Developments, and a new chapter, Computers of the Future, is added. 

The book is intended as an introduction to the computer field for the non-tech- 
nical reader. It can adequately fulfill this purpose. However, the book has a number 
of deficiencies. It also appears that the author has not really kept abreast with mod- 
ern developments in this field, nor with the rapid advances which are taking place in 
this country. In his brief chapter, Computers of the Future, he discusses mainly 
some future potential applications of computers rather than the exciting develop- 
ments which are taking place in the field of computer design. He sometimes betrays 
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a rather limited appreciation of the flexibility and power of modern digital com- 
puters. For instance, on page 214 he states: “It now appears that digital computers, 
too, may have some applications in the field of simulation, particularly when the 
system to be simulated has a digital nature.”’ Digital computers have for years been 
used to simulate very complex physical systems, for example, the temperature dis- 
tribution within the core of a nuclear reactor. 

This book can serve as a simple and not too technical introduction to computers. 
In the opinion of the reviewer, however, the interested reader can find a number of 
other books recently published that will more adequately serve this purpose, and 
probably at a lower price. 


H. P. 


1. T. E. Ivaux, Editor, Electronic Computers, Philosophical Library, New York, RMT 
118, MTAC, v. 12, 1958, p. 255. 


32[Z|—_Watrer J. Karpius, Analog Simulation, Solution of Field Problems, Me- 
Graw-Hill, New York, 1958, xv + 434 p., 23 em. Price $10.00. 


Analog Simulation by Karplus is essentially unique. There are neither competi- 
tive contemporary treatises nor older works with which to compare it in review. It 
almost automatically, then, belongs in any complete library of machine computa- 
tion. The real question, on the other hand, is whether a real need exists for either a 
text or a reference book on its subject matter. As a research engineer, the reviewer 
has a strong affirmative answer to the second half of this question. He can now throw 
away the collection of reprirts and notes which he has been carrying around for 
years. Its merit as a text is more debatable. 

The subject matter of the book should be made clear first. The title, Analog 
Simulation, conveys to most American engineers devices like aircraft and missile 
simulators and electronic amplifier differential analyzers. Even the subtitle, “‘Solu- 
tion of Field Problems,’’ evokes a momentary picture of a fire control computer in 
the “‘field.”’ The full title is correctly descriptive of the category of physical prob- 
lems considered, but a title like “Analog Techniques for Partial Differential Equa- 
tions’ would have made the material covered more readily apparent. 

The book is divided into three major parts. Part 1 (Chapters 2-4) gives 98 pages 
of mathematical background for analog study of field problems. In Part 2 (Chapters 
5-10), the actual “hardware” is discussed in 165 pages. Part 3 returns to a mathe- 
matically organized discussion of the applications of the previously discussed analog 
techniques to different classes of partial differential equations (Chapters 11-14, 115 
pages). It may be seen that most of the book is problem-oriented. 

In spite of this problem orientation, the book is highly recommended as a refer- 
ence on analog techniques. This is in large part a consequence of the superb list of 
references at the end of each chapter and the excellent 22-page general bibliography. 
The only serious omission is Dynamical Analogies by H. F . Olson, the closest equiva- 
lent to Karplus, although almost solely concerned with acoustics and vibration. 
If for no other reason, it would be recommended for a place on the engineer’s book- 
shelf alongside good source books on digital computers and electronic differential 
analyzers. 
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The reviewer feels that information on direct analog simulators should be more 
readily available to engineers. The practitioners of the computing art tend to develop 
rather strong prejudices for the techniques with which they are more familiar. In 
fact, the three major classes of computers (digital data processors, electronic differ- 
ential analyzers, and direct analog devices) all figure with equal importance in cur- 
rent engineering research work for which the reviewer is responsible. Not only do we 
have the continual rather pointless argument about the relative merits of digital 
versus analog computation, but computer engineers working with operational am- 
plifier machines tend to disparage passive-network and other direct analog devices. 
The Introduction (Chapter 1) and Chapter 9 on ‘Electronic Analog Computers” do 
an excellent job of putting this controversy into proper perspective. The analog 
systems discussed in Part 2 are not diminishing in importance to engineering. These 
are: Chapter 5, “Conductive-Solid Analogs”; Chapter 6, ‘Conductive Liquids— 
The Electrolytic Tank”; Chapter 7, ‘Resistance Networks’; Chapter 8, ‘‘Resist- 
ance-Reactance Network Analogs”; Chapter 10, ‘‘Nonelectric-analog Simulation 
Systems.” The value of Analog Simulation as a unique reference on distributed 
system analogs is indisputable. 

Some consideration must be given to the merits of the book as a student text, 
since it is intended for such by the author. The first proposed use for the second half 
of a one-semester upper-undergraduate or first-graduate introduction to analog 
computation does not appeal to the reviewer. At this level, I have strong prejudices 
in favor of orientation along the line of “‘mathematical and physical principles of 
engineering analysis,”’ with all techniques of computation introduced for illustration 
and some familiarization, plus some other topics such as scale modeling and dimen- 
sional analysis. My principal reason for objecting to a course based on Karplus is 
that I feel mathematics beyond an engineering course in ordinary differential equa- 
tions is required; probably at least an introduction to complex variables or the usual 
survey course in applied mathematics. As a text for a full-semester course for grad- 
uate students, however, it appears excellent. Part 1, ‘“The Mathematical Model,”’ 
would be an excellent “refresher” on the mathematical and physical principles in- 
volved. Graduate engineering students working in mechanical engineering areas 
should benefit by the familiarization with electrical engineering techniques. 

The reviewer’s major criticisms of the book are very few. The weakest part is 
Chapter 3, ‘Transformations.’ Although it is a reasonably adequate condensed 
summary with useful reference tables of transformations, as a text for someone with 
no previous introduction to the subject it would be useless. The caution on page 66 
about applicability of conformal transformations in space, although quite correct, 
would only be mystifying after a four-page discussion of the subject. The only mis- 
print of any consequence noted is also in this chapter, namely F(¢) for f(é) in for- 
mula (3.23). There are also occasional implications that physical intuition can be 
substituted for mathematical understanding, as on pages 11, 12, where it is stated 
that “It is entirely possible to develop analogs without taking a mathematical route 
—namely, by applying direct physical insight.” Again, on page 79, Karplus warns 
that use of finite-difference expansions requires an understanding of errors and ap- 
proximations inherent in the process, and conditions for a valid solution then im- 
mediately suggest that physical reasoning and insight may be substituted for mathe- 
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matical understanding. To the reviewer, a contradiction is implied, which, it is 
doubtful, was intended. It is considered dangerous to suggest to students that phys- 
ical intuition may be used in solving problems, without parallel mathematical 
analysis of the properties of the analogous systems involved. Otherwise, there is 
nothing that can be cited to detract from the value of Analog Simulation as a unique 
contribution to its field. 


W. H. Keen 
Bureau of Naval Weapons 


Navy Department 
Washington, D. C. 


33(Z].—W. Orre.t, Compiler, Anwendung von Rechenmaschinen bei der Berechnung 
von Regelvorgingen (Application of Computing Machines to the Computation of 
Servomechanisms), R. Oldenbourg, Verlag, Munich, 1958, 128 p., 24 em. Price 
DM 16.80. 


This book consists of ten papers presented at a meeting in Diisseldorf on No- 
vember 8, 1957. The authors and titles are: 


1. E. Bucovics On the Use of Small Digital Computers for the Solu- 
tion of Fundamental Problems of Servomechanics 

2. J. B. Reswick A Simple Graphical Method for Deconvolution 

3. A. Ieonhard A Special Computer for Polynomials 

4. R. Herschel On the Design of Analog Circuits for Problems of Ser- 
vomechanies 

5. O. Foellinger & Comparison of Computations for Servomechanisms 

G. Schneider using Computing Machines of Different Types 
6. E. Buehler On the Mechanical System with Friction and its Elec- 


tronic Equivalent 

7. H. Witsenhausen Application of Analog Computers for Optimizing Dis- 
continuous Control Circuits with Randomly Vary- 
ing Input Parameters 


8. D. Ernst Practical Work with Analog Computers 
9. Th. Stein On the Usefulness of Analog Experiments in Practical 
Applications 
10. W. Roth Investigation of the Control Characteristics of Electric 


Generator Sets by the Use of Analog Computers 
(All papers are in German, except the second.) 

Only two papers consider digital computers: No. 1 describes in much detail the 
use of an IBM 604 calculating punch for some basic problems; No. 5 compares the 
accuracy and speed of two analog computers (one mechanical, one electronic) with 
that obtained on the IBM 650 in the BELL interpretative mode. The comparison is 
unfair, however, since the Runge-Kutta method is used for solving a system of linear 
differential equations with constant coefficients (time: 27 sec./step for a system of 
six equations!). . 

There is no space here for reviewing each paper separately, so that the following 
statements may do some injustice to individual papers. Most authors report mainly 
on their practical work and experiences, not on new theoretical results. There are a 











112 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


fair number of references, but some are geographically biased and/or obsolete. The 
book gives an interesting picture of the local situation (Germany and Austria) in 
the mid-fifties; but as far as applications of computing machines, and in particular 
digital computers, are concerned, it was almost obsolete when it appeared in 1958 
and is, of course, even more nearly obsolete today. 


Hans J. MAEHLY 
Syracuse University 


Syracuse 10, N. Y. 
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TABLE ERRATA 


296.—K. MircuE.., “Tables of the function [ —log|1 — yI dy, with an account 


of some properties of this and related functions,” Phil. Mag., v. 40, 1949, p. 
351-368. 
On page 359, Table Il: in column 4, first entry, for 0.05 6039 292 read 0.05 
0639 292. 
R. P. Brumpacu, Jr. 
University of Notre Dame 
Notre Dame, Indiana 


297.—G. W. Spencetey, R. M. Spencetey, & E. R. Epperson, Smithsonian 
Logarithmic Tables to Base e and Base 10, The Smithsonian Institution, Wash- 
ington, D. C., 1952. 

The following errata supplement those recently reported by me (Math. Comp., 

v. 14, p. 308, MTE 283). 


For Read 
p. 5 log, 245 ... 98481 147 .. 98481 146 
p. 10 log. 472 ... 85553 ... .. 85555 ... 
p. 176 log, 8772 9.07832 ... 9.07932 ... 
p. 185 log, 9215 9.12856 ... 9.12858 ... 
p. 287 logio 4213 63459 ... 62459 ... 
p. 330 logie 6381 .. 87448 ... S060 


CHARLES R. Sexton 
2947 Elmwood Court 


Berkeley 5, California 


298.—I. M. Vinocrapov, Elements of Number Theory, Translated from the fifth 

edition by Saul Kravetz. Dover Publications, Inc., New York, 1954. 

The following typographical error has been found, supplementing errata previ- 
ously reported [1]: 

p. 227 for 5373, read 3373. 
Sipney Kravitz 

592 Herrick Drive 
Dover, New Jersey 


1. MTAC, v. 9, 1955, p. 125, RMT 59; v. 12, 1958, p. 108, MTE 260. 


CORRIGENDA 


Hersert FE. Sauzer, “A note on the. solution of quartic equations,” Math. 
Comp. v. 14, 1960, p. 279-281. 
p. 280, line 8, . . depending upon whether m is real... . 
should read .. depending upon whether n is real... . 


H. E. Sauzer 
Convair Astronautics 


San Diego, California 


113 











114 CORRIGENDA 


Dierricu Suscuowk, “Explicit formulae for 25 of the associated Legendre 
functions of the second kind,” MTAC, v. 13, 1959, p. 303-305. 
On p. 305, the first recurrence relation should read: 
zx 


Qei*(z) = (m — k)(n + k + L)Qs(2) — 2k + 1) Fo O(a) 


JoHN HERNDON 
Stanford Research Institute 
Menlo Park, Califcrnia 
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